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I .   INTRODUCTION 

A.    BACKGROUND 

The  end  of  the  cold  war  started  a  transition  from  the 
Major  Regional  Conflicts  to  the  Lesser  Intensity  Conflicts 
(Hughes,  1993, pp.1)  as  the  more  likely  scenario  to  be  devel- 
oped in  the  future  for  US  Forces.  Such  a  transition  renewed 
the  interest  and  efforts  to  study  what  is  called  "Littoral 
Warfare".  Other  allied  countries  have  had  the  Littoral  Warfare 
as  their  natural  scope  and  permanent  frame  for  the  design  and 
development  of  the  capabilities  of  their  naval  forces. 

Above  all,  Littoral  Warfighting  implies  sensible  changes 
in  the  physical  nature  of  the  properties  of  the  threat  and  the 
environment  which  together  with  the  mission  and  the  current 
capabilities  determines  the  imbalance  that  drives  research, 
development  and  acquisition. 

1.    The  environment  in  Littoral  Warfare 

The  changes  in  the  physical  properties  of  the  environment 
due  to  the  transition  from  the  Open  Sea  to  Littoral  or 
Coastal  waters  have  already  been  thoroughly  studied  and  are 
characterized  by  the  distinction  between  deep  and  shallow 
waters,  often  termed  blue  and  green  waters  respectively.  A 
particular  maritime  area  is  considered  a  shallow  water  en- 
vironment in  terms  of  the  influence  its  boundaries  exert  in 
its  physical  properties  rather  than  in  terms  of  a  particular 
predetermined  depth. 


From  the  standpoint  of  the  Search  and  Detection  problem, 
which  remains  crucial  to  Littoral  Warfighting  (Hughes, 1986, pp- 
126) ,  the  most  remarkable  differences  between  shallow  or  deep 
water  environments  are  the  intensity  of  noise  and  clutter  and 
the  spatial  and  temporal  variability  in  the  properties  of  the 
ocean  which  become  distinctive  and  determinant  factors 
affecting  the  performance  of  the  sensors  and  the  decision  aids 
associated  with  them.  Several  factors  have  a  direct  impact  in 
increasing  the  intensities  of  clutter  and  noise.  The  most 
important  factors  are  the  wind,  the  biological  content  of  the 
water,  the  proximity  of  the  bottom,  and  human  activity. 

The  wind  increases  both  ambient  noise  and  clutter.  It  has 
been  observed  that  for  the  typical  frequencies  used  in  passive 
sonar,  the  ambient  noise  in  coastal  locations  is  in  general 
more  intense  in  shallow  waters.  The  simplest  reason  is  the 
larger  mass  of  absorption  in  the  deep  ocean.  In  shallow  waters 
the  noise  produced  by  the  wind  (white  caps,  bubbles,  waves, 
etc.)  is  partially  reflected  by  the  bottom  and  the  surface, 
hence  the  dissipation  process  is  slower.  The  increase  in 
clutter  is  due  to  the  roughness  of  the  surface  which  in  turn 
increases  as  the  wave  length  decreases  due  to  the  reduction  of 
depth  (Urick,1975,  pp.  236).  Additionally,  particular  condi- 
tions of  short  waves  typical  of  shallow  waters  like  tide  rips, 
choppy  seas,  etc.,  make  detection  harder  due  to  the  movement 
transmitted  to  the  platforms  where  the  sensors  are  carried. 

Biological  activity  increases  as  depth  decreases  thus 


leading  to  an  increase  in  both  noise  (crabs,  shrimps, etc. )  and 
clutter  (whales,  schools  of  fish,  zooplankton,  etc. ) . (Urick, 
1975, pp. 193) 

Human  activity  produces  both  noise  and  clutter.  Shipping 
activity  increases  close  to  shore.  A  fishing  vessel,  for 
instance,  can  be  considered  a  source  of  both  noise  and 
clutter. 

Finally,  the  closer  bottom  of  the  ocean  not  only  increas- 
es the  backscattering  due  to  reflection  in  the  case  of  an 
active  sensor,  but  also  prevents  the  formation  of  caustics  or 
convergence  zones  typical  of  deep  waters . 

Meteorological  and  Oceanographical  changes  also  intensify 
when  transitioning  from  the  macroscale  of  the  open  ocean  to 
the  mesoscale  of  coastal  regions.  For  instance,  daily  varia- 
tions in  the  temperature  profile  in  the  water  and  in  the  air 
lead  to  drastic  changes  in  the  propagation  of  electromagnetic 
and  acoustic  energy. 

The  prediction  of  the  performance  of  a  sensor  for  a 
particular  target  and  environment  is  of  the  essence  in  the 
searching  problem.  The  more  intense  the  variability  in  the 
conditions  found  in  shallow  or  coastal  waters  the  harder  the 
prediction.  Consequently,  the  uncertainty  in  the  evaluation  of 
the  effectiveness  of  a  given  detection  effort  increases. 

2 .    Changes  in  the  threat 

With  the  end  of  the  cold  war,  the  most  likely  threat 


became  smaller  but  not  weaker,  at  least  from  a  conventional 
warfare  point  of  view.  Small  fast  patrol  boats  can  strike  fast 
with  remarkable  lethality  and  are  not  easy  to  detect  in 
choppy  seas .  Modern  conventional  submarines  have  the  same  or 
even  better  information  systems  than  older  nuclear  submarines, 
taking  advantage  of  multiple  sensors  and  massive  signal 
processing  based  on  inexpensive  and  commercially  available 
technology.  Except  for  the  snorkeling  cycle  where  they  become 
noisy  and  visible,  conventional  submarines  on  batteries  are 
almost  undetectable  by  current  passive  sensors.  Their  smaller 
size  makes  it  more  difficult  to  detect  them  with  active 
sensors  due  to  a  reduction  in  target  strength. 

3 .    Sensors  and  Tactical  Decision  Aids (TDA) 

The  threat  and  the  environment  are  determinant  factors  in 
the  decision  to  use  a  particular  type  of  sensor  and  whether  it 
will  be  used  in  a  passive  or  active  fashion.  Regardless  of  its 
type  or  mode,  a  sensor  normally  requires  a  Tactical  Decision 
Aid  for  a  more  intelligent  and  efficient  use. 

In  the  realm  of  the  ASW  Search  Problem,  the  fact  that  a 
diesel  submarine  on  batteries  is  almost  undetectable  with 
passive  acoustic  means,  especially  in  the  noisy  shallow  water 
environment,  drives  the  need  to  use  active  sensors.  This 
conclusion  can  be  extended  to  other  threat -environment 
combinations  in  shallow  waters,  for  instance  small  fast 
surface  combatants  because  they  normally  operate  in  silence 


unless  they  are  about  to  strike . 

A  sensor  typically  feeds  a  Tactical  Data  System  where 
software  tools  (TDA)  help  decide  the  best  course  of  action 
given  a  particular  profile  of  mission,  environment  and  threat. 
The  essence  of  a  Tactical  Decision  Aid  is  providing  qualita- 
tive and  quantitative  processing  of  the  data,  so  it  becomes 
information.  The  factors  of  the  situation  should  be  ranked, 
simplified  and  presented  in  a  way  such  that  the  determinant 
ones  are  emphasized  in  terms  of  the  mission. 

The  smaller  signal  to  noise  and  target  to  clutter  ratios 
due  to  the  more  noisy  and  cluttered  environment  will  not  only 
make  detection  much  more  difficult  but  ,also,  together  with 
strong  environmental  variability,  will  generate  more  uncer- 
tainty about  the  effectiveness  of  the  search.  As  a  conse- 
quence, the  need  for  a  stochastic  TDA  to  keep  adequate  track 
of  the  uncertainties  is  evident.  Figure  1.1  summarizes  the 
problem  structure  leading  to  the  development  of  a  TDA.  The 
chains  of  causes  and  effects  are  displayed  as  a  flow  pattern. 

B.    THE  PROBLEM 

For  any  Navy  that  does  not  have  a  TDA  capable  of  handling 
the  uncertainty,  the  problem  is  apparent  and  the  need  to 
develop  one  is  clear. 

A  large  array  of  TDA' s  in  the  U.S.  NAVY  like  VPCAS, 
ASWTDA  and  PACSEARCH  (Wagner,  1989,pp . II . 2)  responded  to  the 
needs  of  the  cold  war  where  the  most  important  means  of 


Figure  1.1.  The  structure  of  the  problem  of  Active  Search  in  Shallow 
water  is  remarkable  different  than  Passive  Search  in  Deep  Water. 


detection  was  the  passive  sonar.  In  a  passive  sonar  search, 
specially  from  airborne  platforms,  the  target  will  usually  not 
be  aware  of  the  presence  of  the  searchers .  As  a  consequence 
those  TDA's  do  not  consider  any  reaction  of  the  target  caused 
by  the  presence  or  activities  of  the  searcher. 

In  an  active  search,  the  target  will  almost  invariably 
know  of  the  presence,  location  and  activities  of  the  searcher 
before  detection,  A  simple  signal  intensity  budget  shows  that 
there  exists  a  distance  where  the  energy  suffices  to  produce 
a  detection  in  the  passive  receptor  of  the  target,  but  it  is 
not  enough  to  produce  the  same  effect  at  the  receptor  of  the 


active  device  in  the  searcher.  Geometrical  spreading  and 
absorption  in  the  medium  account  for  the  difference.  The 
active  searcher  is  thus  a  beacon  with  the  choice  of  avoidance 
or  approach  up  to  the  target .  Bearings  from  the  target  to  the 
searchers  change  noticeably  due  to  the  shorter  distances  and 
higher  speed  of  the  actions,  thus  enabling  a  relatively 
accurate  Motion  Analysis  (Wagner,  1989,pp. III-2)  by  the 
target . 

As  a  consequence,  the  target  is  not  only  aware  of  the 
presence  of  the  searcher  but  also  can  have  a  very  clear 
picture  of  the  situation.  It  will  invariably  have  a  mission 
driven  reaction,  either  to  attack  or  evade.  Information 
becomes  crucial  and  the  target  has  the  advantage.  This  feature 
of  the  target  is  the  essence  of  the  problem,  hence  the  name 
RTCAS,  Reactive  Target  Computer  Assisted  Search. 

Time  and  distance  scales  change  as  well  when  transi- 
tioning from  a  passive  to  an  active  search  in  the  case  of 
sonar.  A  passive  search  is  a  resource  employed  in  a  larger 
time  frame  and  over  a  more  extended  area  than  the  active 
search.  Hence  the  name  Area  ASW  for  the  actions  typically 
carried  out  in  the  open  ocean  either  to  "sanitize"  a  corridor 
or  to  survey  a  given  wide  area  for  strategic  purposes. 
Conversely,  active  search  in  ASW  is  a  focal,  more  concentrated 
effort.  The  detection  ranges  are  shorter  hence  the  space  scale 
of  the  search  is  smaller.  The  cost  in  terms  of  allocation  of 
resources  is  higher  in  an  active  search,  hence  its  sustainabi- 


lity  decreases  and  the  time  scales  shrink.  Minutes  count. 

C.    APPROACHING  A  SOLUTION 

The  objective  of  this  thesis  is  to  develop  a  Reactive 
Target  Tactical  Decision  Aid  for  Computer  Assisted  Search. 

TDA'  s  are  aimed  to  assimilate  and  conveniently  present 
the  data  (target,  own  assets  and  environment)  to  the  decision 
maker,  i.e.,  transform  the  data  into  information  and  "...ana- 
lyze the  tactical  problem  beyond  what  is  possible  by  humans  in 
a  timely  fashion..."  (Wagner,  1989,  pp.I.l) .  That  analysis  may 
either  yield  an  evaluative  picture  from  where  the  decision 
maker  will  construct  a  course  of  action  more  easily  or  produce 
a  recommended  course  of  action.  The  first  step  is  always 
needed  and  crucial  to  the  second,  so  this  thesis  is  con- 
strained to  an  evaluative  TDA  leaving  the  assessment  capabili- 
ties as  a  future  enhancement. 

Given  the  uncertainties  involved  in  the  problem  and 
following  approximately  the  same  approach  as  previous  TDA' s 
did,  RTCAS  uses  a  Monte  Carlo  Simulation  to  represent  a 
Reactive  Target,  a  Fixed  Scan  Model  to  measure  the  detection 
effort  and  a  Bayesian  Filtering  Process  to  fuse  the  simulation 
with  the  real  world  search  activity.  This  section  describes 
those  models. 

1.    The  Information  Fusion  Paradigm 

Figure  1.2  is  a  pictorial  description  of  the  information 
fusion  paradigm  applied  to  a  search  TDA.  Two  main  sources  of 


information  can  be  distinguished,  reality  and  a  Monte  Carlo 
simulation.  Reality  provides  the  information  from  the  tactical 
situation,  the  datum,  the  time,  and  searcher  activity.  The 
Monte  Carlo  Simulation  provides  a  large  number  of  possible 
locations  and  movement  of  the  target  based  on  the  datum,  each 
one  equally  likely,  at  least  at  the  beginning,  to  be  the  real 
world  target:  the  "prior." 

Fusing  the  information  consists  of  updating  the  "prior" 
or  probability  map  of  the  location  of  the  target  using  all 
available  information.  There  are  three  main  causes  of  updat- 
ing, namely  the  passage  of  time,  negative  information  and 
positive  information.  In  the  first  case  the  prior  will  move 
with  the  variety  of  courses  and  speeds  provided  by  the  Monte 
Carlo  Simulated  Targets.  A  large  number  of  simulated  tracks 
(500  in  RTCAS)  resemble  many  of  the  possible  positions  and 
movements  of  the  real  target.  If  no  searching  were  conducted 
the  "furthest  on  circle"  would  just  expand  continuously  until 
none  of  the  500  simulated  targets  were  left  in  the  area. 

Updating  for  negative  information  takes  into  account  the 
information  embedded  in  an  unsuccessful  search  over  a  period 
of  time.  During  a  certain  time  interval,  some  simulated 
targets  are  within  the  range  of  detection  of  some  searcher.  If 
the  real  target  is  not  detected,  then  those  simulated  tracks 
become  less  likely  and  those  further  apart  from  the  searcher 
become  more  likely. 

Updating  for  positive  information  takes  place  when  a 


The  Information  Fusion  Paradigm 
Environment 


Figure  1.2.  The  information  Fusion  Paradigm 

searcher  gains  a  contact .  Those  simulated  targets  within  a 
certain  distance  from  where  the  contact  is  produced  become 
more  likely  than  others  further  away.  The  probability  map,  the 
posterior  or  "new  prior, "  consists  of  a  partition  of  the 
Local  Area  into  a  number  of  cells.  Each  cell  has  a  probability 
value  which  results  from  adding  up  all  the  individual  likeli- 
hoods or  weights  of  all  simulated  tracks  within  the  cell. 
RTCAS  produces  a  visual  representation  of  the  prior  by  mapping 
cell  probabilities  to  a  color  code. 

This  successive  updating  of  the  individual  weights  of  the 
simulated  targets  (Bayes  Theorem)  and  the  cell  probabilities 
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(Law  of  Total  Probability)  at  a  given  time  interval  is  called 
Bayesian  Filtering  (Stone  et  al . , 1995,pp. 6) . 

2.    Detection  Effort 

In  a  deterministic  detection  model  like  the  "cookie 
cutter"  sensor  (Washburn, 1989) ,  a  target  within  a  certain 
range  is  detected  with  probability  one,  while  those  outside 
are  not  detected  at  all .  The  uncertainties  involved  in  the 
detection  process,  intensified  in  the  case  of  a  shallow  water 
environment  and  a  small  target,  prevent  the  use  of  a  determin- 
istic model  to  evaluate  the  detection  effort.  Instead, the 
Fixed  Scan  Model  (Washburn,  1989, pp. 4)  ,  is  used.  For  each  pair 
simulated  target -searcher,  the  Sonar  Equation  has  to  be  solved 
yielding  a  value  of  Signal  Excess  (Kinsler  et  al . ,  1976, pp. 
411)  that  represents  how  much  above  the  Detection  Threshold 
the  signal  return  from  the  target  is.  However,  the  Signal 
Excess  is  considered  a  random  variable  due  to  the  uncertain- 
ties in  the  problem,  so  a  detection  probability  is  calculated 
representing  the  detection  effort.  For  example,  a  simulated 
target  very  close  to  the  searcher  will  have  a  large  value  of 
Signal  Excess  chiefly  because  of  the  small  transmission 
losses,  thus  leading  to  a  high  detection  probability. 

The  same  detection  model  serves  to  quantitatively 
appraise  the  effort  in  an  unsuccessful  search  over  an  area  as 
well  as  the  likelihood  of  a  contact  report.  The  estimate  of 
such  an  effort  is  given  by  the  nondetection  probability  in  the 
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former  and  the  detection  probability  in  the  latter. 

3 .    The  Reactive  Target 

Active  search  turns  the  one  sided  activity  of  passive 
search  into  a  game  of  two  opposed  goals.  Parameterizing  the 
reaction  of  the  target  is  a  complicated  task  because  it 
involves  human  decisions.  Selecting  a  course,  depth,  and  speed 
or  releasing  a  decoy  or  a  weapon  are  the  result  of  a  sophisti- 
cated mental  process  that  takes  place  in  the  target .  The 
information  is  integrated  in  time  to  provide  a  picture  of  the 
threat  the  target  faces  versus  its  own  resources  to  accomplish 
a  goal.  Modeling  such  a  problem  is  a  task  out  of  the  scope  of 
this  thesis.  However,  there  are  some  basic  parameters  such  as 
course  and  speed  that  can  be  considered  the  drivers  of  the 
situation.  RTCAS  uses  those  basic  parameters  to  construct  a 
simple  model  of  the  behavior  of  a  real  target . 

D.    THESIS  STRUCTURE 

Chapter  II  develops  the  core  models  for  the  Reactive 
Target,  the  Bayesian  Filtering  and  the  Detection  Effort. 

Chapter  III  explains  the  implementation  of  the  models. 

Chapter  IV  verifies  the  correctness  and  robustness  in  the 
implementation  and  the  significance  of  a  Reactive  Target  Model 
as  opposed  to  a  Non  Reactive  one. 

Chapter  V  summarizes  future  enhancements  of  the  model  and 
the  implementation. 

Chapter  VI  states  the  conclusions  and  recommendations. 
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II.  THE  ANALYTICAL  MODEL 

A.    THE  TARGET  MODEL 

Three  basic  tasks  have  to  be  accomplished  in  the  target 
model:  generation,  updating  and  reaction. 

1.    Track  Generation 

The  initial  prior  in  RTCAS  is  provided  by  a  single 
scenario  given  by  a  lost  contact  or  a  flaming  datum,  i.e.,  the 
cue  about  the  position  of  the  target  given  by  a  ship  in  flames 
as  a  consequence  of  a  torpedo  attack. 

Other  TDA  models, (Wagner, 1989, pp. II .3) ,  allow  the 
generation  of  multiscenario  priors  but  those  are  associated 
with  large  scale  models  that  deal  mainly  with  passive  Area 
Search  where  several  initial  positions  and  estimated  movements 
are  likely  and  should  be  taken  into  account.  In  the  case  of  a 
flaming  datum,  the  expected  movement  is  to  flee  the  area. 

This  initial  prior  is  represented  by  a  bundle  of  500 
simulated  tracks  that  encompass  the  many  possibilities  for 
target  initial  position  and  movement. 

A  Datum  can  be  defined  by  the  following  elements 

•  Datum  identity, 

•  time, 

•  position, 

•  error, 

•  course  and  speed. 
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RTCAS  will  deal  with  a  single  datum,  hence  the  identifi- 
cation of  the  datum  is  irrelevant.  In  a  Local  Area  Search,  in 
general,  all  time  and  space  coordinates  are  considered  rela- 
tives to  that  of  the  datum,  i.e.,  the  geographical  location  of 
the  area  and  the  real  time  are  not  part  of  problem. 

The  time  scale  starts  with  the  searchers'  arrival  in  the 
Local  Area  after  a  delay  time   t^.  (Figure  2.1) 


^ 

DATUM  TIME 

zr 

0 

L,^ 

J > 

delay 

^                            (/me 

\\     SEARCHERS' 
^    ARRIVAL 

Figure  2.1.  The  time  scale. 

The  positions  of  the  tracks  representing  the  datum  are 
referred  to  the  center  of  a  grid  that  covers  the  Local  Area 
where  the  search  is  conducted.  The  coordinates  of  any  track  i 
are  [xTi,yTi]  (t)  at  time  t. 

The  error  is  taken  into  account  when  the  tracks  are 
generated  and  is  assumed  to  be  bivariate  normal, 

XTi(-T<i)  ~  N  (/ix  ,  o^)  (2.1) 

where  the  mean  \Lx  =  \Iy  are  assumed  to  be  equal  to  the  center 
of  the  Local  Area  and  the  standard  deviations  0^=0^  are  input 
by  the  user.  The  notation  X  ~  N()  means  that  X  is  a  random 
variable  normally  distributed. 
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The  magnitude  of  each  target's  velocity  vector  is  given 
by  a  random  variable  uniformly  distributed  between  a  lower  and 
an  upper  bound  given  by  the  user  and  does  not  change  unless 
there  is  a  reaction  because  of  the  presence  of  the  searcher. 
(See  Section  3) . 

The  course  is  also  a  random  variable  uniformly  distribut- 
ed between  user  input  limits.  The  purpose  of  this  model  is  to 
resemble  a  situation  where  a  preferred  course  of  evasion 
exists, e.g.,  windward,  towards  its  mission  essential  target, 
navigational  hazard  avoidance,  etc.  The  situation  where  no 
preferences  may  be  estimated  about  the  target's  intentions 
uses  the  same  model  including  all  possible  courses  between  0 
and  3  60  degrees,  thus  leading  to  an  omnidirectional  fleeing 
datum . 

The  delay  time  r  is  taken  into  account  by  updating  the 
position  of  all  tracks  using  each  one's  velocity  vector  to 
obtain  the  initial  target  distribution. 

[Xj.^,yj.i]  (0)  =[Xj,^,Xj.J  (-T^) +Tv^(-t^)         (2.2) 

A  typical  target  distribution  at  time  zero  can  be  seen  in 
Figure  2.2  representing  the  scenario  the  searchers  find  when 
they  arrive  at  the  Local  Area.  To  the  left  of  the  figure,  the 
target  is  assumed  to  flee  omnidirectionally.  To  the  right  is 
the  case  when  the  target's  course  is  between  10  and  100 
degrees . 
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Figure    2.2.    The   initial  target  distribution. 

2 .         Track  Updating 

At  every  time  increment  r  the  position  of  all  the  targets 
as  a  function  of  time  is  updated  using  each  track's  velocity 
vector. 


[Xj,j,yj,j]  ( t+x)  =  [Xj,j,Xj.^] +TV^  ( t) 


(2.3) 


The  speed  of  the  target  remains  unchanged  unless  there  is 
a  reaction. 

3 .    Track  Reaction 

The  reactive  updating  of  the  tracks  consists  of  changes 
in  the  properties  of  each  track  as  a  consequence  of  the 
presence  of  the  searchers.  All  reactions  of  a  real  target  may 
be  reflected  in  a  variation  of  speed,  course  or  depth,  RTCAS 
updates  the  course  only.  The  variations  in  depth  are  more 
likely  to  occur  in  deep  water  environments.  Speed  alterations, 
in  general,  are  made  within  ranges  that  do  not  have  a  great 
impact  in  the  overall  geometry  and  kinetics  of  the  problem. 
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The  target  has  two  goals,  evade  the  datum  and  evade  the 
searchers.  In  RTCAS  the  simulated  targets  evade  the  datum 
unless  they  have  to  evade  a  searcher  because  of  its  proximity. 
The  assumption  is  made  that  the  target  knows  the  searcher's 
location.  This  is  not  too  strong  an  assumption  in  active 
search.  In  general,  the  most  common  information  the  target 
knows  about  the  searchers  is  the  bearing.  In  the  close 
distances  normally  implied  in  Local  Area  actions,  the  bearing 
varies  rapidly  thus  allowing  the  target  to  perform  an  adequate 
passive  motion  analysis  (Wagner,  1989, pp.  III-l)  .  It  is  assumed 
that  the  target  reacts  when  a  searcher  at  a  distance  smaller 
than  some  threshold. 

The  simplest  evasion  maneuver  is  to  adopt  a  course  that 
generates  a  direction  of  relative  motion  (DRM)  opposed  to  that 
of  collision.  Figure  2.3  ,  left  depicts  such  a  situation.  Most 
of  the  time,  this  will  be  impossible  for  the  target  because 
its  speed  will  be  smaller  than  that  of  the  searcher  (Figure 
2.3,  right)  .  A  discussion  of  all  possible  courses  of  action  a- 
vailable  to  the  target  at  this  point  is  out  of  the  scope  of 
this  thesis.  However,  there  are  some  constraints  on  what  the 
target  should  do.  It  should  not  continue  at  the  same  course  if 
the  distance  is  decreasing,  nor  should  it  turn  and  end  up  with 
a  relative  bearing  of  ±90  degrees  to  the  searcher,  since  this 
would  increase  dramatically  the  target  strength  in  the  sonar 
equation.  A  180  degrees  relative  bearing  should  also  be 
avoided  since  this  may  imply  a  reduction  in  the  quality  of 
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Figure    2.3.    The  evasion  alternatives  as  a  f iinction  of  the  velocities 
of   the   searcher  Vg  and  the   target  "7^. 

information  obtained  by  the  target  because  of  the  blind  sector 
in  the  stern.  For  simplicity,  it  is  assumed  that  target  i 
turns  a  fixed  angle  dd  e  (90°, 180°)  off  the  bearing  B^^  to  the 
searcher.    This  new  course   for  target   i   is. 


dB 


(2.4) 


where  Db^g/dt  is  the  bearing  rate.  If  the  bearing  increases 
(decreases) ,  the  target  turns  to  a  course  i^0  degrees  measured 
counterclockwisely  (clockwisely)  from  B^^.  It  is  the  tangential 
component  of  the  searcher's  speed  relative  to  the  target's 
velocity  vector  that  determines  the  sign  of  the  last  term  in 
(2.4)  and  it  is  calculated  in  RTCAS  with  the  following 
formula. 


sign 


dB, 


JS  N  _ 


dt 


)  =sign {sin  (Q^-B^^)  ) 


(2.5) 


where  9,  is  the  searcher's  course . (Figure  2.4) 
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Figure  2.4.  The  evasive  course  parameters 

When  facing  several  searchers  within  the  reaction  range, 
the  target  decides  to  evade  the  closest. 


B. 


THE  BAYESIAN  FILTERING  MODEL 


1.    The  Probability  Map 

The  final  purpose  of  the  model  is  to  provide  a  "probabil- 
ity map"  representing  the  location  of  the  target  as  an  aid  in 
conducting  a  Local  Area  Antisubmarine  Active  Search. 

The  search  starts  with  a  "flaming  datum"  (Cheong,  1988, 
pp.1) ,  hence  a  relatively  small  time  and  space  frame  contains 
the  action,  e.g.,  1  to  4  hours  and  an  area  with  a  radius  from 
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10  to  20  nautical  miles. 

The  probability  map  is  a  partition  of  the  Local  Area  into 
a  grid  of  KxL  cells  each  of  which  has  a  probability  value 
associated  with  it.  We  define  the  following  events: 

I   :{  "The   target   is   within   the  Local  Area"} 

Gki  :{  "The   target   is   in   the  grid  cell   k,l"} 

Then  the  probability  map  is  the  bivariate  probability 
mass  function  of  P(Gj^2\I).  The  idea  of  a  flaming  datum  as  the 
starting  point  of  the  search  makes  the  condition  on  the  fact 
that  the  target  in  within  the  Local  Area  rather  irrelevant  at 
the  start.  However,  as  times  goes  by,  the  influence  of  I 
increases  as  the  target  may  have  already  left  the  box  where 
the  search  is  conducted. 

In  order  to  build  the  probability  map  it  is  necessary  to 
find  a  way  to  assign  a  value  to  P(G),i\I). 

This  Bayesian  analysis  begins  with  a  single  scenario 
prior  distribution  of  the  "state  of  nature"  (Washburn, 
1995, pp.1)  represented  hy  N  =  500  possible  tracks.  Intuitive- 
ly, when  conducting  a  Local  Area  search,  one  make  some 
assumptions  and  evaluate  a  series  of  "what  if"  questions, 
reacting  to  the  most  evident  and  easy  to  maintain  in  one's 
mind  as  a  chain  of  causes  and  effects.  In  this  case  N  simulat- 
ed tracks  try  to  keep  that  chain  alive  for  as  long  as  possi- 
ble. 

The  prior  is  generated  by  the  Target  Model  as  described 
in  the  previous  section.  The  output  of  this  model  is,  at  this 
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stage,  the  stochastic  location  and  movement  of  each  one  of  the 
N  simulated  tracks  according  to  the  assumed  prior.  We  define 
the  event, 

T^:  {"The  real    target  follows   track   i",  i=l..N} 
then,  if  all  tracks  are  equally  likely  to  be  the  real  target, 
assumption  embedded  in  the  idea  of  a  single  scenario,  then  the 
initial  "weight"  of  each  track  is 

P[r,](0)=%  (2.6) 

N 

and  the   initial  probability  map  can  be   depicted  as 

P[G,,]  it)=     X;      P[T,]{t)  (2.7) 

ieKjci  it) 

where   K,,i(t)    =    {i    :    [xTi,yTiJ  (t)    =    G;,^,    t=0}. 

Ideally,  if  there  were  no  searchers  and  (2.7)  was 
calculated  as  a  function  of  time,  it  would  convey  almost  the 
same  information  as  the  expanding  furthest-on-circle  (Cheong, 
1988,  pp.1)  .  However,  the  purpose  of  RTCAS  is  to  integrate  the 
information  from  the  real  world  provided  by  the  assets 
conducting  the  search  together  with  the  information  provided 
by  the  prior  (Stone  et  al . , 1995,pp. 1) .  The  search  done  up  to 
a  certain  time  t  may  be  unsuccessful  or  result  in  one  or  more 
searchers  obtaining  a  contact,  which  requires  updating  the 
prior  for  negative  or  positive  information  respectively  (Wag- 
ner, 1989, pp. II .5) . 
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2.    Updating  for  Negative  Information. 

The  valuable  information  of  an  unsuccessful  search  during 
the  interval  r  is  that  the  real  target  is  less  likely  to  be  in 
the  cells  that  had  been  searched  as  a  result  of  the  detection 
effort  imposed  to  the  tracks  located  in  or  near  those  cells. 
The  value  of  each  track  weight  will  be  decreased  in  accordance 
with  that  effort.  The  track  weights  will  change  from  the  prior 
values  at  time  t  to  the  posterior  values  at  time  t+r,  hence 
changing  the  probability  map  of  the  location  of  the  target. 

Define  the  event, 

Dg^  :{      "The     target     is     detected    by    the     searcher    s 
during      time  interval    r"  } 
then 

cdp^^^{t)=P[D^^\T^]   (t)  (2.8) 

is  the  cumulative  probability  that  the  target  was  detected  by 
the  searcher  s  during  the  time  interval  t  starting  at  time  t 
conditioned  on  that  the  real  target  was  located  at  the 
position  of  simulated  target  i  during  that  interval.  This 
value  is  provided  by  the  Detection  Model  to  be  developed  in 
section  C. 

If  D's^  is  the  complement  of  D^^    ,    then 

P[D^,\T^]   (t)=l-ccfp^i,(t)  (2.9) 

is  the  nondetection  probability  of  searcher  s  over  target  i 
during  interval  r  starting  at  time  t.  Consequently,  the 
probability  that  the  target  is  following  track  i    and  is  not 
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detected  by  searcher  s   during  the  interval  r   is, 

P[TfUD^,]{t)=P[Di,\T^]   (t)P[rJ  it)  (2.10) 

Using  Bayes  Theorem,  the  Negative  Information  Updating  of  the 
prior  is  , 

ieN 

Note  the  convention  that  P  [Ti]  (t+r)  is  conditioned  on  all 
information  available  up  to  time  t+r. 

The  denominator  in  (2.11)  represents  the  unconditional 
probability  P[D' ^^1  (t)  obtained  using  the  Total  Probability 
Theorem.  Furthermore  it  is  the  normalizing  factor  that  will 
allow  the  new  prior  to  be  a  "proper"  probability  mass  function 
i.e.  E,,„  P[TJ  (t+r)    =  I. 

An  iterative  application  of  this  process  permits  the 
calculation  of  the  new  prior  for  any  time.  The  update  2.11 
must  be  applied  once  per  time  interval  for  each  searcher  that 
does  not  detect  the  target . 

3.    Updating  for  Positive  Information 

The  probability  map  should  be  updated  for  Positive 
Information  when  a  searcher  reports  a  contact .  It  may  be 
contended  that  once  a  contact  is  obtained,  the  relative  value 
of  the  infonnation  provided  by  the  simulation  of  tracks  is 
somehow  negligible  with  respect  to  the  real  world  information 
implied  in  the  contact  report.  However,  in  a  shallow  water 
acoustic  environment  where  many  elements  of  the  "clutter" 
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space, e.g.,  a  whale,  a  rock,  an  eddies,  may  be  considered  a 
contact,  the  fusion  of  information  of  both  sources  remains  in 
place  as  a  means  to  smoothing  the  reactions  and  providing 
continuity  whenever  the  contact  is  confirmed  not  to  be  the 
desired  target . 

Define  the  following  events 

CgcT  '  {  "The     searcher     s     reports     a     contact     at     the 

position    c    during    the    interval     t",  s=1..S, 
c=l.  .C} 
where  S   is  the  total  number  of  searchers  and  C  is  the  total 
number  of  reported  contacts  among  all  searchers,  and, 

Eg^  :{      "The  real    target  is  located  at   the  position  of 
the   contact  report   c  during  the   time  interval 

T"} 

then, 

C,^=D,,nE,  (2.12) 

where,  as  before,  D^^  is  the  event  that  searcher  s  detects  the 
target.  Thus, 

P[C^„\Ti]  {t)=P[D^^\T^]  it)P[E^\D^^,Ti]  it)  (2.13) 

The  first  factor  in  2.13  is  the  detection  effort  alloted 
by  searcher  s  at  the  position  of  simulated  track  i  during 
period  r, 

P[D^,\E^]  it)  =cdp^,^  it)  (2.14) 
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Assuming  that  location  errors  are  circularly  normally  distrib- 
uted the  second  factor  is 


P[E^\D,,,T^]=ke 


2a: 


(2.15) 


where  r^^^  is  the  distance  between  the  position  of  the  contact 
report  c  and  track  i  during  interval  t  (Figure  2.5),  Jc  is  a 
constant  independent  of  c  and  i,  and  o^^  is  the  standard 
deviation  of  reporting  errors  made  by  searcher  s  for  a  target 
following  track  i  during  interval  t.  If  S2g  is  the  beamwidth  of 
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Figure    2.5.    The   Error   in   the    contact   report . 
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the  sensor  and  r^^    is  the  distance  between  searcher  s    and 
tracki,  then  RTCAS  takes  a^^   to  be  Q  ^sJ2.    Thus  2.15  becomes, 


P\E^\D^,Tr\=e 


^^  (2.16) 


■S-^Si 


where  the  constant  Jc  was  dropped  because  it  cancels  when  Bayes 
Theorem  is  applied.  Substituting  2.14  and  2.16  into  2.13, 


PiC^„\T,\  {t)=cdp,,,it)e 


Qlr^At)^  (2.17) 


represents  the  probability  that  searcher  s  reports  a  contact 
at  position  c  given  that  the  target  follows  track  i  during  the 
interval  t.  Equation  2.17  embodies  the  mathematical  core  of 
fusing  information  between  the  real  world  contact  report  and 
the  simulation. 

The  probability  that  the  real  target  follows  track  i 
given  that  the  searcher  s  reported  a  contact  at  the  position 
of  c  during  period  t  applying  Bayes  Theorem  is  (Washburn, 
1995, pp. 6), 


i=l 


J2p[Cs„'t,]p[t,]  it) 

which  using  2 . 17  becomes 


(2.18) 


cdD  ■  e    "s-^s-i  prr.!  (t) 

2  2 


(2.19) 


N 


Y,cdp,,,e    ^^^-Plr,]  (t) 
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where  the  +  superscript  implies  that  the  contact  is  actually 
a  true  contact  on  the  target.  In  practice,  updates  for 
positive  information  must  allow  for  the  possibility  of  false 
alarms,  so  the  formula  actually  used  by  RTCAS  is, 

p[rj  (t+T)=p[rj*(t+T)C)+P[rj  it)  (i-o)  (2. 20) 

where  Q  is  the  "credibility"  (Washburn,  1995, pp.  2)  of  the 
contact  report.  When  Q=l,  2.20  is  the  same  as  2.19,  but  when 
Q=0,    2.20  is  just  the  prior  P[T^](t). 

It  might  be  contended  that  the  posterior  alternate  to 
P [Tj^]*  (t+T )  should  include  the  negative  information  implied  in 
the  unsuccesful  search  that  a  failed  contact  report  conveys. 
However,  it  is  very  plausible  and  commonly  seen  in  the 
practice  that  a  searcher  investigating  a  contact  will  focus 
onto  it,  not  only  because  of  a  natural  tendency  in  the 
distribution  of  the  attention  of  the  sonarman,  but  also 
because  investigating  a  contact  requires  more  often  than  not 
using  modes  or  techniques  in  the  sensors  that  prevent  the 
operator  from  an  adequate  scanning  of  the  rest  of  the  picture 
even  without  intention  of  doing  so.  Examples  of  those  modes  or 
techniques  are  doppler  and  image  analysis,  noise  reduction  by 
windowing  the  signal  processing  in  the  target,  angular 
scanning  reduction  spotting  to  concentrate  acoustic  intensity 
on  the  target,  etc.  As  a  consequence  and  for  computational 
convenience,  RTCAS  considers  an  absolute  focusing  on  the 
contact,  i.e., the  individual  track  probability  alternative  to 
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the  one  updated  for  positive  information  is  the  value  at  the 
beginning  of  the  scan  period  as  in  equation  2.20. 

4.    The  undetected  probability  mass 

Although  each  probability  map  is  the  best  estimate  of  the 
location  of  the  target,  it  says  nothing  about  how  much  of  the 
original  probability  mass  remains  undetected.  A  measure  of  the 
effect  of  the  search  is  necessary  since  the  searcher  needs  to 
know  how  much  probability  mass  is  still  undetected  in  order  to 
be  able  to  decide  when  to  terminate  the  search. 

The  denominator  in  (2.11)  is  the  probability  that  the 
real  target  is  following  one  of  the  500  simulated  tracks  and 
was  not  detected,  hence  is  the  undetected  probability  mass 
during  the  interval  t  at  time  t+r: 

5,(t)=^P[r,nDi,]  (2.21) 

for  instance,  at  time  t=0,  before  any  search  is  conducted 
S(t)=l.  At  each  updating  step  the  searching  will  reduce  the 
original  probability  mass  such  that  at  any  time  t,  the 
Cumulative  Undetected  Probability  Mass  is 

U{t)=U{t-T)S,{t)  (2.22) 

This  recursive  definition  allows  the  calculation  of  its  value 
at  every  time  step  in  order  to  produce  the  desired  assessment . 

C.    THE  DETECTION  MODEL 

Several  possibilities  were  taken  into  account  in  select- 
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ing  the  detection  model.  The  idea  is  to  model  active  search. 
The  Fixed  Scan  Model  seems  to  be  naturally  applicable  since 
(Washburn, 1994,  pp.  3): 

•  Each  pulse  can  be  thought  of  as  an  independent 
look, 

•  Detection  in  the  absence  of  a  pulse  is  impossible, 

•  There  exists  a  pulse  rate  which  can  be  associated 
with  the  pulse  repetition  frequency  of  a  sonar. 

The  first  statement  is  seemingly  too  strong  an  assumption 
but  in  RTCAS  only  the  conditional  independence  of  the  pulses 
given  the  target  track,  is  required.  The  validity  of  the 
assumption  depends  on  the  rapidity  of  fluctuations  in  the 
sonar,  equation  terms  relative  to  the  time  between  active 
pulses . 

The  pulse  repetition  frequency  (prf )  is  normally  associ- 
ated with  the  maximum  range  at  which  it  is  desired  to  attempt 
a  detection.  Given  a  typical  value  of  maximum  scale  range  in 
active  sonar  of  about  18000m,  the  round  trip  of  a  pulse  takes 
24  seconds  at  a  sound  speed  of  1500m/s,  thus  yielding  a  prf  of 
approximately  2  per  minute.  In  conventional  active  sonar  the 
frequency  is  normally  greater  than  4  kilohertz.  The  strongest 
time  variability  in  the  propagation  is  given  by  the  time 
change  in  the  refraction  patterns  which  in  turn  depend  on  the 
depth- velocity  profile.  In  shallow  waters,  the  fastest  change 
in  the  upper  layer  of  the  velocity  profile  can  be  regarded  as 
semi-diurnal  as  a  consequence  of   the  heat   interchanges 
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provoked  by  the  daily  cycles  of  the  sun.  Nonetheless,  the 
effect  of  cloud  cover,  rain,  wind  or  other  many  factors  may 
produce  noticeable  changes  in  the  propagation  conditions. 
Furthermore,  there  are  some  changes  in  the  geometry  of  the 
detection  problem  that  may  cause  variations  in  the  time  frame 
of  interest.  For  instance,  given  a  frequency  of  4  . 5  kilohertz, 
the  wavelength  will  be  0.3  meters.  In  shallow  waters,  it  is 
very  likely  that  the  receiver  will  sense  a  multipath  scatter- 
ing from  the  target.  If  the  difference  between  two  paths  is  a 
multiple  of  0.15  meters,  there  will  be  destructive  interfer- 
ence thus  leading  to  an  eventual  lack  of  detection.  Although 
a  deterministic  geometric  model  may  account  for  each  situation 
regarding  the  relative  positions  of  target  and  searcher  with 
respect  to  the  bottom,  the  surface  and  the  layer,  the  outcome 
of  such  a  model  may  be  regarded  as  randomly  distributed  with 
respect  to  space  and  time.  Other  phenomena  may  concur  in  order 
to  increase  the  randomness  in  the  signal  at  the  receiver 
between  different  pulses,  for  instance,  the  presence  of  the 
bubbles  due  to  the  wake  of  other  ships  or  the  target,  changes 
in  the  angle  of  incidence  of  the  pulse  at  the  target,  the 
destructive  interference  with  the  signals  from  other  plat- 
forms, the  macroscopic  biological  activity,  etc.  All  of  these 
sources  of  randomness  are  typical  of  the  shallow  water  active 
signal  propagation  problem. 

The  assumption  in  RTCAS  is  that  each  pulse  is  condition- 
ally independent  of  the  others.  Of  course  the  computational 
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difficulty  in  considering  nonindependence  had  to  be  taken  into 
account  and  favored  the  chosen  option. 

1.    The  Fixed  Scan  Model 

This  thesis  applies  the  Fix  Scan  Model  as  a  paradigm  of 
the  use  of  hull  mounted  and  dipping  active  sonars  and  sono- 
bouys  in  searching  for  a  conventional  submarine.  However,  it 
is  applicable  to  other  types  of  sensors  like  radar,  lidar, 
photo  electronic  devices,  etc.  as  well  as  other  types  of 
targets  like  surface  ships,  life  rafts,  land  vehicles,  etc. 

The  sensor  is  assumed  to  emit  pulses  at  a  repetition 
frequency,  prf ,    given  by 

prf=—  (2.23) 

M 

where  C  is  the  speed  of  the  sound  in  the  water  and  M  is  the 
maximum  unambiguous  range,  assumed  always  to  be  equal  to  the 
maximum  range  scale  selected  in  the  device. 

The  probability  of  detection  will  be  calculated  only  for 
those  targets  within  the  maximum  range  scale  of  any  searcher 
at  the  end  of  every  updating  time  interval  r.  The  value  of  t 
must  satisfy  a  compromise  between  the  need  of  an  adequate 
"sampling"  of  the  search  process  and  the  computational 
difficulty  of  too  many  calculations  in  a  given  time  frame. 

As  an  example,  for  a  maximum  range  scale  of  18000m,  the 
interval  between  pulses  will  be  24  seconds,  which  may  be 
approximated  as  3  0  seconds.  Hence,  if  the  updating  interval  is 


31 


T  =  5  minutes,  10  independent  scans  are  going  to  be  made  with 
the  same  sonar  operating  characteristics,  i.e.,  same  pulse 
type,  prf ,  frequency,  etc.  Selecting  5  minutes  as  the  updating 
interval  is  convenient  from  the  point  of  view  of  the  computa- 
tional effort  as  well  as  tactically  reasonable;  e.g.,  the  rate 
of  turn  of  a  conventional  submarine  at  6  knots  is  such  that 
it  would  take  about  4  minutes  to  complete  a  180°  turn. 

During  one  updating  interval  t,  both  the  searcher  and  the 
target  are  assumed  to  remain  stationary  at  a  point  midway 
between  their  respective  locations  at  the  beginning  and  at  the 
end  of  T  (Figure  2.6),  for  purposes  of  detection  probability 
calculations. 

Assuming  the  above  mentioned  independence  between  each 
scan,  (Washburn, 1995,pp. 3) ,  the  probability  that  at  least  one 
of  the  pulses  k  from  searcher  s  detects  the  target  given  that 
it  follows  track  i,  at  time  t,  is. 


cdp,,,  ( t)  =l-nL  (l-P^si(t))  (2.24) 


where  T,     the  total  number  of  pulses  during  the  updating 
interval  t,  may  be  calculated  as, 

T=x*prf  (2.25) 

The  detection  probability  of  each  pulse  k  is  given  by, 

p^^.(t)=$(^?5i!^)  (2.26) 

where  <l>  is  the  Standard  Normal  Cumulative  Density  Function  and 
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Figure  2.6.  The  kinematics  of  the  detection  process 

SEg^  (t) ,  calculated  by  the  Propagation  Model  to  be  described 
in  the  next  section,  is  the  Signal  Excess  from  target  i  at  the 
receiver  in  searcher  s  at  time  t,  and  a^  is  the  standard 
deviation  in  the  calculation  of  SE^^   due  to 

•  the  randomness  in  the  detection  process,  and 

•  the  cumulative  error  or  uncertainties  in  the  calcu- 
lation of  each  one  of  the  terms  in  the  sonar  equa- 
tion. 

The  value  of  o^j  can  be  estimated  either  by  experience  or 
by  using  a  property  of  the  variance  of  the  sum  of  random 
variables  (Devore, 1991,pp. 102) , 
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"-^ 


.2  (2.27) 
k 

k=l 


E'i 


where  a,^  is  the  error  in  the  calculation  or  estimation  of  each 
k  of  the  S  terms  in  the  sonar  equation.  Wagner  (1989, pp.  II- 
28)  suggests  values  of  6-9  dB  for  a^',  RTCAS  uses  15  dB  since 
it  is  intended  for  use  with  relatively  poor  sonar  forecasts. 

2 .    The  Propagation  Model 

This  model  corresponds  to  the  propagation  of  active  sonar 
signals  in  a  shallow  water  environment .  The  values  of  all 
terms  in  the  sonar  equation  will  be  either  assumed  to  be 
constants  or  approximated  by  simple  models.  It  is  the  user 
that  will  finally  input  the  values  corresponding  to  more 
accurate  calculations. 

For  each  pair  searcher  s  and  target  i  at  time  t  the 
signal  excess  is  be  calculated  using  the  active  sonar  equa- 
tion, (Kinsler  et  al .  1982)   limited  either  by  noise, 

SEN^^=SL^-DNL^  ( t)  -DT^+DI^-2  TL^^  ( t)  +TS^^  (t)  (2.28) 

or  by  reverberation, 

SER^^  ( t)  =TS^^  ( t)  -DT^-SS^+TLG^^  ( t)  -lOlog.o  (-^)  (2  .29) 

All  levels  in  the  formulae  above  are  in  decibels  with 
respect  to  one  microPascal .  The  Source  Level  5L,  the  Detection 
Threshold  DT    and  the  Directivity  Index  Dl    are  considered 
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constant  for  each  searcher  s.  The  rest  of  the  terms  are 
variables  calculated  by  RTCAS  and  are  developed  in  the 
following  sections.  DNL^  is  the  Detected  Noise  Level  for 
searcher  s,  TL^^  (t)  is  the  Transmission  Loss  for  each  pair 
(s,i)  at  time  t,  TS^iit)  is  the  Target  Strength  of  target  i 
with  respect  to  searcher  s,  SS  is  the  Scattering  Strength  for 
the  sonar  in  searcher  s ,  TLG^^  (t)  is  the  transmission  loss  due 
to  geometrical  spreading  and  the  last  term  in  (2,10)  is  the 
area  of  the  reverberating  surface  within  the  beam  of  the 
transmitted  pulse  corresponding  to  the  speed  of  the  sound  c, 
the  pulse  length  t  and  the  horizontal  angle  of  the  sonar  beam 
Q  (Kinsler  et  al.,1892).  The  smaller  of  (2.28)  and  (2.29)  is 
to  be  selected  and  used  in  calculating  (2.26) . 

In  active  search  in  shallow  waters, the  sonar  equation  is, 
in  general,  more  likely  to  be  limited  by  reverberation  rather 
than  by  noise.  However,  in  some  cases  the  effect  of  the  noise 
will  prevail  at  distances  within  the  maximum  range  scale  since 
the  echo  level  and  the  reverberation  level  decrease  with  range 
but  the  noise  level  remains  the  same .  A  typical  behavior  can 
be  seen  in  Figure  2.7,  where  up  to  7500  meters  the  noise 
restricted  Signal  Excess  is  higher  than  the  reverberation 
restricted  one. 

The  following  sections  discuss  the  terms  in  the  sonar 
equations  (2.9)  and  (2.10)  that  are  not  considered  a  constant. 
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Figure  2.7.  Noise  and  Reverberation  Limited  detection. 


3.    Target  Strength 

Target  strength  is  defined  as  the  echo  returned  by  the 
target  and  hence  is  a  factor  when  dealing  with  active  sonar. 
The  corresponding  term  in  the  sonar  equation  is: 


r5^,=101og,o(-i),r=l 

-'■zo 


(2.30) 


where  J^  is  the  intensity  of  return  at  1  yd  from  the  target 
and  I^  is  the  incident  intensity. 

The  intensity  ratio  Ir/^i     is  a  function  of  distance 

between  the  target  and  the  source  (Urick,  1975, pp. 284) ,  the 

shape  of  the  target  (Urick, 1975, pp. 274) ,  the  pulse  length 

(Urick, 1975, pp. 285) ,  the  frequency  of  the  incident  signal 
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(Urick, 1975,pp.283) ,  and  the  geometry  of  the  reflection  (U- 
rick, 1975,pp. 281) .  The  influence  of  all  those  factors  except 
the  last  can  be  considered  constant. 

In  the  case  of  a  monostatic  sonar  the  direction  of 
incidence  is  the  same  as  the  direction  in  which  the  signal  is 
received  at  the  transducer  except  for  small  variations  due  to 
the  relative  motion  target -searcher  during  the  time  the  pulse 
travels.  The  target  strength  is  a  function  of  the  angle 
between  the  heading  of  the  submarine  and  the  angle  of  inci- 
dence of  the  signal,  which  is  called  aspect  angle.  If  lateral 
symmetry  is  assumed,  then  the  aspect  angle  can  be  considered 
from  0  to  ±  TT,  with  0  being  at  the  bow.  The  great  variation  of 
TS  as  a  function  of  the  aspect  angle  requires  the  aspect  angle 
to  be  considered.  The  Target  Strength  will,  in  general,  vary 
from  its  maximum  value  at  broadside  to  a  minimum  at  the  bow 
and  stern. 

Let  r  be  the  distance  target -searcher  and  v  the  target 
speed,  then, 

cos(0)=f-v-  (2.31) 

indicates  the  dot  product  between  the  corresponding  unit 
vectors  and  0  is  the  smallest  angle  between  them,  in  this  case 
the  aspect  angle  (Figure  2.8). 

The  typical  values  of  TS  vary  from  about  25  dB  at  beam 
aspect  to  about  15  dB  at  bow  and  stern  aspects  (Urick,  1975 
p. 282)  .  The  following  equation  is  an  approximation  adopted  as 
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a  simplified  model  of  the  variations  between  those  extremes, 

TS{Q)=TS^+ATS*{l-cos^{d))  (2.32) 

where  TS^,  is  the  value  of  Target  Strength  at  bow  and  stern, 
ATS  is  the  typical  increment  to  reach  the  value  at  beam  aspect 
and  (1-cos^  (B) )  models  the  variation  as  a  function  of  the 
aspect  angle  G.  Substituting  (2.12)  into  (2.13)  yields, 

TS^^  ( t)  =TS^^ATSil-{f^^  ( t)  'V^  ( t)  )  2)        (2.33) 

where  r^^  (t)  is  the  instantaneous  distance  between  searcher  s 
and  target  i,  and  Vi(t)  is  the  speed  of  target  i  at  time  t. 
RTCAS  uses  a  square  of  the  cosine  as  a  model  because  of 


Figure  2.8.  Vector  representation  of  the  target  strength 
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computational  convenience.  Note  in  2.33  that  the  trigonometric 
function  in  2.32  is  replaced  by  a  dot  product,  which  can  be 
calculated  using  the  components  of  the  corresponding  unit  vec- 
tors instead  of  the  trigonometric  function. 

4.    Transmission  Loss 

Transmission  loss  involves  the  combined  effect  of  geomet- 
rical spreading  and  absorption. 

For  simplicity,  assuming  that  the  distance  r^^  (t)  is  much 
larger  than  the  depth  H,  geometrical  spreading  is  calculated 
with  the  following  isospeed  shallow-water  channel  perfectly 
rigid  bottom  expression  (Kinsler, 1984,pp.429) , 

rLG^^(t)=101og[r^^(t)+^)]  (2.34) 

where  the  last  term  is  the  correction  for  the  transition 
range . 

The  total  Transmission  Loss  is 

TL^^  ( t)  =TLG^^  ( t)  +ar^^  ( t)  (2  .  35) 

where  a,  the  absorption  coefficient,  is  calculated  with 

^   8x10-5  ^    0.04   ^.      Tn-7^2 

a  = + +4x10  'f'^  ,^    ^^. 

0.7  ^j^   6000  ^^  (2.36) 

f2  f2 

from  Kinsler  et  al . (1976, pp. 398)  where  f  is  the  frequency  in 
Hertz  at  which  each  sonar  operates. 
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5.    Noise  Level 

The  Noise  Level  accounts  for  the  background  ambient  noise 
(an)    and  the  platform  self  noise  (sn) , 


NL=2Qlo<3 


10 


an+sn 


(2.37) 


IjiPa 

Instead  of  calculating  the  intensities  an  and  sn  their 
corresponding  levels  ANL  and  SNL  in  decibels  are  interpolated 
from  observational  data. 

The  Ambient  Noise  Level  is  calculated  by  means  of  a 
third  order  polynomial  fit  on  windspeed  and  a  linear  fit  on 
the  logarithm  of  frequency  of  observational  data  obtained  from 
curves  representing  the  noise  spectra  for  some  coastal 
locations  in  Figure  7.7.  in  Urick  (1975, pp. 191), 

AZVL=46  .12+2  .22W-0  .OAlw^-ll  .01[logio(f)  -3]     (2.38) 

where  w  is  the  wind  speed  in  m/s  and  f  is  the  frequency  in 
hertz.  This  approximation  is  centered  around  f  =  1000  Hz. 
Figures  2.9  and  2.10  compare  Equation  2.38  to  the  data  from 
the  above  mentioned  reference . 

The  Self  Noise  Level  is  calculated  as  a  function  of 
frequency  and  the  searcher's  speed  using  figure  11.11  in 
Urick(1975,  pp.340),  which  shows  equivalent  isotropic  self- 
noise  levels  at  25  khz  on  a  number  of  World  War  II  American 
and  British  destroyers.  A  polynomial  fit  of  this  data, 
including  an  extrapolation  term  for  other  frequencies  assuming 
a  slope  of  -6dB  per  octave,  yields. 
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Figure     2.9.    The    Ambient    Noise    as    a    function    of    frequency    and    the 
windspeed  as   a  parameter  between   1   and  16   knots . 


SNL=23 +2.  02v  At)  +20log(-^) 
^  f 


(2.39) 


where  v^  is  the  searcher's  speed  in  knots  and  f   is  the 
frequency  in  kilohertz.  Figure  2.11  compares  this  equation 
with  the  corresponding  data.  SNL  and  ANL  are  power- summed  to 
obtain  NL. 

6.    Scattering  Strength 

According  to  Urick  (1975,p.253)  ,  when  downward  refraction 
occurs,  i.e.,  negative  sound  speed  profile,  bottom  back 
scattering  is  dominant  and  the  reverberation  level  can  be 
calculated  using  the  bottom  scattering  strength.  Conversely, 
when  upward  refraction  occurs,  i.e.,  positive  sound  speed 
profile  as  in  isothermal  water,  the  surface  backscattering 
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Figure      2.10.     The    Ambient    Noise     as     a     function    of     frequency    and 
windspeed  as  a  parameter  between  17  and  35  knots. 

becomes  most   important.    This  model   follows  that   rule. 

Following     Urick      (1975, pp. 241) , the     Surface     Scattering 
strength   is   calculated  with, 

55L=10logio(fxi2xsin(e))°-55-45.3  (2.40) 

where  f  is  the  frequency  in  hertz,  6  is  the  grazing  angle  in 
degrees  and  h  =  0.0026  l>f^^  where  W  is  the  windspeed 
in  meters  per  second.  In  a  shallow  water  environment , the 
backscattering  of  the  sonar  signal  will  occur  at  all  possible 
grazing  angles  due  to  the  multiple  interactions  between  the 
surface  and  the  bottom.  Due  to  the  analytical  and  computation- 
al difficulty  in  accounting  for  those  angles,  9  was  chosen  to 
be  35°,  an  average  between  the  typical  limit  angles  in  Figure 
8.21  in  Urick(1975,  pp.  240). 
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Figure  2.11.  Self  Noise  data  for  25kHz,  its  polynomial  fit  and  the 
extrapolations  for  other  frequencies  using  a  slope  of  -6dB  per  octave. 

In  the  calculation  of  bottom  strength,  Urick  (1970, pp. 
396)  found  that  the  spectrum  slope  in  shallow  waters  was 
similar  to  measurements  made  in  deep  waters,  therefore  figure 
8.30  in  Urick  (1975,  pp.  248)  could  be  used  as  an  approxima- 
tion for  shallow  waters.  Depending  on  three  generic  types  of 
bottom,  the  backscattering  Bs  is  calculated  using  a  third 
order  polynomial  fit  of  the  data  in  figure  8.30  as  a  function 
of  frequency  as  follows 

•  TYPE  I  :  Plain  bottom,  little  or  no  roughness,  high 
angular  and  frequency  variations. 


BS=0 . 012f  2+0 . 283f 2+0 . 0062f-28 . 153 


(2.41) 


•  TYPE  II:  Intermediate  between  a  very  rough,  rocky 
bottom  and  a  Type  I  bottom 
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BS=0  .  076/2-0  .  233/2+2  .  8706/-40  .  41 


(2.42) 


•  TYPE  III:  Heavily  dissected  bottom,  with  underwater 
ridges,  very  rough.  Almost  no  angular  or  frequency  dependence. 


B^=-18 


(2.43) 


Figures  2  .  12  and  2  .  13  show  the  data  and  the  corresponding 
polynomial  fit  for  bottom  types  I  and  II.  The  user  must  decide 
which  type  of  bottom  to  use. 
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Figure  2.12.  The  bottom  backscattering  for  plain  surfaces . 
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Figure  2  .  13  .  The  bottom  backscattering  for  intermediate  surfaces 
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III.  THE  IMPLEMENTATION  OF  THE  MODEL 


A. 


PROGRAM  STRUCTURE 


The  analytical  model  described  in  Chapter  II  was  imple- 
mented in  a  structural  procedural  program  written  in  Turbo 
Pascal  7.0  for  DOS  (Borland,  1991).  The  structure  of  the 
program  can  be  observed  in  Figure  3.1. 

The  program  consists  of  a  main  program  named  RTCAS 
(Reactive  Target  Computer  Assisted  Search)  and  seven  units 
namely,  InterManager,  ScoutManager,  TargtManager, Detection 
Manager,  CellsManager  and  GraphManager . 


RTCAS  STRUCTURE 


Figure  3.1.  The  Program  Structure  of  RTCAS 
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The  main  program  controls  the  sequence  in  which  each  one 
of  the  "Managers"  will  execute  their  tasks. 

A  single  data  structure  contains  all  the  information 
about  the  simulated  tracks,  the  searchers,  the  environment  and 
the  simulation  control  required  to  run  the  program. 

All  units  share  and  perform  their  tasks  in  a  common  data 
structure.  The  units  could  be  changed  provided  their  inputs 
and  outputs  remain  the  same,  thus  allowing  future  expansions 
or  modifications  in  the  software. 

The  implementation  of  the  Data  Fusion  Model  by  the  main 
program  RTCAS  is  described  in  the  algorithm  shown  in  Figure 
3.2.  The  code  of  RTCAS  is  in  APPENDIX  Al . 


RTCAS  general  algorithm 

Input  initial   simulation  data 
Input   initial    searchers   data 
Loop  until   selected  time 

Get   information   from  Searchers 

Loop   for  all    searchers 

If  searcher  is   in   contact 

Update  for  Positive   Information 

else 

Update  for  Negative   Information 

end 

Normalize  Tracks'    Weights  with  Bayes  Theorem 
Update   Traks'   movements 
Update   cells  probabilities  with   Total   Probability  Law 
Output  Probability  Map 


Figure     3.2.     The    algorithm    implemented    in    the    main    program    of    RTCAS 
organizes   the   sequence   or   the   component   siibprograms . 
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B.    DATA  STRUCTURE 

1.  Tracks  Data 

The  number  of  tracks  was  fixed  at  500  for  computational 
convenience.  Processing  two  searchers  during  150  minutes  of 
activity  in  the  Local  Area  took  2  minutes  running  in  a 
Personal  Computer  with  a  4  86  DX2  66  MHZ  processor,  which 
implies  a  time  compression  ratio  of  1/75.  Any  future  increase 
in  computational  speed  may  permit  an  increase  in  the  number  of 
simulated  targets  in  order  to  augment  the  statistical  richness 
of  the  simulation. 

Each  track  is  represented  by  a  record  with  fields  for 
position,  course  and  speed  and  weight.  Information  common  to 
all  targets  is  stored  out  of  the  records  for  memory  space 
convenience,  e.g.,  standard  deviation  in  the  initial  posi- 
tions, speed  and  course  ranges,  reaction  distance  and  reaction 
turn  angle . 

All  the  information  concerning  the  tracks  is  input  by  the 
user  through  the  InterManager' s  corresponding  procedures. 

2.  Searchers'  Data 

The  program  can  process  up  to  15  searchers.  However  it  is 
very  unlikely  that  such  a  large  number  of  assets  will  be 
concentrated  for  a  Local  Area  or  Lost  Contact  search. 

Each  searcher  is  represented  by  a  record  with  fields  for 
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position,  course,  speed,  sonar  parameters  like  frequency, 
status,  source  level,  etc.  and  information  regarding  the 
contact . 

The  searchers'  data  structure  is  composed  of  static  data 
like  sonar  frequency,  detection  threshold,  etc.,  input  by  the 
user  and  controlled  by  the  InterManager  unit  and  by  dynamic 
data  like  position,  speed,  sonar  range  selected,  etc.  read 
from  a  file  by  a  procedure  in  the  ScoutManager  unit .  Reading 
from  a  file  resembles  receiving  the  dynamic  searchers' 
information  through  a  data  link  or  from  the  main  stream  of  a 
tactical  data  system. 

3 .  Environmental  data 

Input  by  the  user  through  the  User  Interface  includes 
bottom  type,  sound  velocity  profile  type,  wind  speed  and 
depth.  This  part  of  the  data  structure  is  used  by  DetecManager 
in   the    determination   of    the    transmission    losses. 

4.  Simulation  control  data 

This  part  of  the  data  structure  stores  information 
required  in  the  control  of  the  program  like  the  name  and 
location  of  input /output  devices,  (files  in  this  implementa- 
tion) ,  times  for  which  probability  maps  are  desired,  and  delay 
in  the  start  of  the  search  from  the  datum  time.  This  informa- 
tion is  input  by  the  user  using  Intermanager  procedures. 
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c. 


ALGORITHM  STRUCTURE 


This  section  outlines  the  algorithmic  structure  executed 
by  the  units  Interemanager,  Targe tManager,  ScoutManager, 
DetecManager  and  CellsManager  following  the  sequence  estab- 
lished by  the  main  program  RTCAS  to  follow  the  Data  Fusion 
Model . 

1.    Interface  Manager 

This  unit  declares  the  complete  common  data  structure, 
creates  the  user  input  or  static  part  of  that  structure  and 
provides  control  of  the  program  by  means  of  a  series  of 
interactive  menus  called  the  User  Interface. 

For  computational  convenience,  inside  the  program  the 
calculations  are  done  in  measurement  units  different  from 
those  input  by  the  user.  The  Interface  Manager  converts  the 
data  from  the  external  to  the  internal  unit  system,  i.e., 
discrete  time  unit  (dtu) ,  pixels,  and  discrete  speed  units 
(dsu)  ,  according  to  Table  3.1.  The  code  for  this  unit,  the 
file  INTERMAN.PAS  is  in  APPENDIX  B  and  the  sequence  of  menus 
is  included  in  APPENDIX  C. 


Internal  Units 

External  Units 

Time 

1  dtu 

5  minutes 

Distance 

1  pixel 

0.1  Nautical  Miles 
(185.2  meters) 

Speed 

1  dsu 

1.2  Knots 

Table  3.1.  The  measurement  units  conversion. 
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2 .    Target  Manager 

This  unit  implements  the  model  described  in  Chapter  II, 
Section  A  for  the  generation,  updating,  and  reaction  of  the 
simulated  targets.  Figure  3.3  depicts  the  process  to  decide  an 
evasion  course  at  the  target.  Note  that  an  evasion  course  is 
calculated  for  each  searcher  (evasion)  and  the  one  correspond- 
ing to  the  closest  searcher  {d(s) <d(s-l))  is  selected  as 
course.  The  code  for  this  unit,  the  file  TARGTMAN . PAS ,  is  in 
APPENDIX  D. 


s  [searcher)  d(s.-1)-  BIG  NUMBER 

i   (sim.  target 


dCs]-distance(i3] 


b-bearingG^] 


b-U(port;starboard] 


evasion 'b 


evasion*  b-'-tack 


evasion  ■  b^aclc 


— <:.^w(s-i)^^> 

course- evasion 

inc(s) 

Figure  3.3.  The  deteirmination  of  the  evasion  course  at  the  target 
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3 .    Scout  Manager 

Implements  the  core  of  the  Data  Fusion  Paradigm  as 
described  in  Chapter  II  section  B. 

The  public  procedures  InitScouts,  CloseScoutsDataLink, 
ConnectWithScouts  and  ReadScoutDataLink  have  the  purpose  of 
controlling  the  input  of  information  from  the  searchers.  In  a 
real  world  implementation  this  procedure  would  receive  the 
information  through  a  Data  Link  device  or  from  the  main  stream 
of  the  Tactical  Data  System  of  the  platform.  In  RTCAS  the 
information  is  read  from  a  file  with  a  format  shown  in  Table 
3.2. 


ID 

SEARCHER 

SONAR 

CONTACT 

LOCATIOM 

MOVEMENT 

ON 

SCALE 

HOT 

LOCATION 

MOVEMENT 

CREDIBILITY 

time 

SI 

214 

205 

336 

10 

1 

10 

0 

- 

- 

- 

25 

HI 

195 

202 

171 

0 

1 

7 

1 

200  250 

120  7 

0.75 

H2 

195 

280 

271 

90 

0 

- 

- 

- 

- 

- 

SI 

204 

225 

336 

10 

1 

10 

0 

- 

- 

. 

30 

HI 

195 

202 

171 

0 

1 

5 

1 

205  245 

120   7 

0.85 

H2 

195 

280 

271 

0 

1 

5 

1 

207  252 

150   4 

0.75 

•• 

... 

... 

... 

• 

• 

• 

• 



... 



Table  3.2.  The  Data  Link  information  simulated  by  the  content  of  a  file. 

The  file,  resembling  the  format  of  a  data  bus  in  a 
communications  link,  has  a  row  for  every  searcher  at  every 
time  updating  period,  fixed  in  the  program  to  5  minutes.  Each 
row  contains  information  regarding  the  identification  of  the 
searcher  as  a  ship,  helicopter  or  sonobuoy,  position,  course 
and  speed,   sonar  state,   sonar  range  scale  selected  and 
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contact  information  like  position,  course,  speed  and  credibil- 
ity. Course  and  speed  are  used  to  keep  track  of  searchers  as 
well  as  contacts  and  also  used  in  the  calculations  of  the 
transmission  losses  by  DetecManager .  For  example,  in  Table 
3.2,  one  surface  ship  and  two  helicopters  are  in  the  Local 
Area.  At  time  25  the  ship  SI  is  emitting  but  does  not  have  a 
contact,  whirlwind  HI  is  dipping  with  its  sonar,  holds  a 
contact  and  assigned  a  credibility  value  of  .75.  Helicopter  H2 
is  in  transit  at  90  knots.  The  inclusion  of  the  sonar  condi- 
tions in  the  "data  link  message"  makes  the  program  generic  for 
any  type  of  platform.  All  that  is  needed  to  fuse  the  informa- 
tion with  the  simulation  is  position  and  whether  or  not  the 
sonar  is  on  or  not.  The  identification  at  the  beginning  of  the 
row  is  required  to  display  an  appropriate  symbol .  The  value  of 
credibility  may  be  considered  as  analogous  to  the  standard 
confidence  level  in  the  Allied  doctrine  (NATO,  1990) . 

The  procedure  ReadScoutsDataLink  updates  all  fields  in 
the  searcher  data  base  and  when  a  contact  is  produced,  the 
corresponding  fields  for  contact  information. 

The  basic  structure  of  the  Data  Fusion  Paradigm  performed 
by  the  ScoutManager  is  shown  in  Figure  3.4. 

The  code  for  ScoutManager  is  in  files  SCOUTMAN.PAS  and  is 
shown  in  APPENDIX  E. 
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Figure  3.4.  The  basic  algorithm  for  updating  the  prior. 
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4 .  Detection  Manager 

The  sole  purpose  of  this  unit  is  the  calculation  of  the 
value  of  cdp5,  the  five  minutes  cumulative  detection  probabil- 
ity as  developed  in  Chapter  2,  Section  C.  The  code  for  this 
unit  is  in  the  file  DETECMAN . PAS  and  is  shown  in  APPENDIX  F. 

5.  Cells  Manager 

This  unit  renders  the  Probability  Map.  The  final  product 
of  the  procedure  UpdateProbinCell  is  a  bidimensional  array 
with  a  probability  value  associated  with  each  cell.  This 
product  is  the  probability  map.  It  could  be  used  in  a  differ- 
ent kind  of  graphics  interface  or  transmitted  by  the  procedure 
BroadcastProbabilityMap.  Figure  3.1  shows  this  output  from 
CellsMan.  In  RTCAS  the  procedure  outputs  its  broadcasting  to 
a  file  that,  in  a  real  world  application,  would  be  replaced  by 
an  Input/Output  device.  The  purpose  of  broadcasting  is  to 
provide  the  same  information  about  the  probability  map  to  all 
searchers  and  other  recipients. 

The  code  for  this  unit  is  in  the  file  CELLSMAN. PAS  and 
can  be  seen  in  APPENDIX  G. 

6.  Graph  Manager 

This  unit  produces  the  final  output  of  the  program  in  the 
form  of  a  color  coded  Probability  Map.  Additionally,  it 
produces  the  rest  of  the  Graphic  Information  Display  (GID) 
elements  like  grids,  color  scale,  static  and  dynamic  informa- 
tion display,  etc.  Figure  3.5  show  the  layout  of  the  GID. 
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The  colors  in  the  probability  map  are  calculated  by  the 
procedure  SetProbabilityColors  using  the  color  scale  shown  in 
the  lower  left  corner  of  the  GID.  The  colors  are  relative  to 
the  largest  probability  value  on  any  cell. 

Instead  of  using  one  of  the  natural  color  scales  like  RGB 
or  RYB  (Elliot  et  al .  ,  1995),  white  was  selected  as  the 
brightest  color  in  correspondence  with  the  maximum  probability 
value  because  of  the  black  background  RTCAS  uses . 

The  GID  shows  the  local  area  with  the  probability  map 
together  with  graphic  information  regarding  the  searchers 
activity.  The  symbology  found  in  most  Tactical  Data  Systems 
was  adopted  to  represent  the  different  types  and  conditions  of 
searchers  and  contacts. 

The  contact  information  is  displayed  with  a  hostile 
submarine  symbol  in  the  map  and  with  information  regarding  the 
time  of  initial  contact,  course,  speed,  searcher  holding  it, 
etc.  in  the  bottom  of  the  screen. 

The  code  for  this  unit  is  in  the  file  GRAPHMAN.PAS  and 
can  be  seen  in  APPENDIX  H. 
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Figure  3.5.  The  Graphic  Information  Display  (GID^ 
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IV.  DATA  ANALYSIS  METHODS:  VERIFICATION  AND  EVALUATION  OF 
THE  MODEL 

The  Data  Analysis  Methods  described  in  this  chapter 
numerically  verify  the  correctness  in  the  implementation  of 
the  Bayesian  Filtering  Model  and  evaluate  the  significance  of 
using  a  Reactive  Target  Model  as  compared  to  using  a  Non 
Reactive  one. 

The  final  product  of  the  model  is  a  probability  map, 
hence  the  correctness  of  the  model  implementation  as  well  as 
the  relative  efficiency  of  different  methodologies  must  be 
evaluated  on  the  model's  ability  to  accurately  construct  the 
probability  distribution  of  the  location  of  the  target. 

The  verification  consists  of: 

•  Selecting  a  "dummy,"  assumed  to  be  the  real  target, 
among  the  simulated  ones. 

•  Running  several  scenarios  with  different  sets  of 
variables. 

•  Extracting  the  test  statistics  of  interest. 

•  Calculating  the  Measures  of  Effectiveness   (MOE) 
based  on  those  statistics 

•  Executing  a  test  to  verify  a  given  hypothesis. 

A.    MEASURES  OF  EFFECTIVENESS 

Accuracy,  Area  of  Uncertainty,  and  Mean  Missed  Distance 
were  selected  as  MOE's  to  respectively  verify  the  correctness, 
dispersion  and  location  of  the  distribution  represented  by  the 
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probability  map  and  the  significance  of  a  Reactive  Target 
Model  with  respect  to  a  Non  Reactive  one. 

Mean  Detection  Probability  and  Mean  Time  to  Detection 
were  selected  to  evaluate  the  performance  of  the  search 
conducted  with  a  Reactive  Target  Model  as  compared  to  a  Non 
Reactive  one  when  the  information  provided  by  the  probability 
map  was  used  in  conducting  a  search.  The  rest  of  this  section 
describes  each  one  of  the  selected  MOE's. 

1 .    Accuracy 

Accuracy  is  a  statistical  measure  of  how  well  the 
probability  map  represents  the  model's  uncertainty  in  the 
target's  location  (Widdis, 1995,pp.  14).  Furthermore,  accuracy 
is  a  measure  of  correctness  in  the  implementation  of  a 
Bayesian  Updating  algorithm.  The  containment  statistic  X  is 
the  cumulative  probability  of  all  cells  G^^-  with  a  higher 
likelihood  than  that  of  the  cell  G  where  the  dummy  is. 

^=  E  ^ij^^  E   ^ij  (4.1) 

(i.j)eK  {i.j)€K' 

where 

K    :{    (i,j):    P(Gij\l)     >    P(G\I)    }, 
K':(    (i,j):    P(G,j\l)     =    P(G\I)     }, 
and  L7  is  a  uniformly  distributed  random  variate  over  the 
interval  [0,1] .  The  second  term  in  (4.1)  represents  a  uniform 
proportion  of  the  all  cells  where  there  is  a  tie, i.e.,  the 
cell  probability  is  equal  to  that  of  the  cell  where  the  dummy 
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is.  Assuming  a  discrete  extension  of  the  Probability  Integral 
Transformation  (Dudewicz  et  al . ,  1988,  pp674)  already  used  by 
Stone  et  al (1991, pp.  12)  and  Widdis  (1995, pp.  15),  the  random 
variable  X  should  be  uniformly  distributed  over  the  interval 
(0,1) . 

The  MOE  described  in  this  Section  is  a  discrete  version 
of  a  one  sample  Kolmogorof -Smirnoff  Test  Statistic  (Dudewicz 
et  al.,  1988, pp.  671).  The  K-S  statistic  represents  the 
maximum  deviation  of  the  empirical  distribution  of  X  from  the 
Uniform, 

D=supre7?7LZ/n^g[Q  -L]  \x-Fix)  \  (4.2) 

where  the  second  term  is  the  sample  cumulative  distribution 
function  of  X.  Accuracy  is  defined  as  (Stone,  1991, p.  12) , 

A=l-D  (4.3) 

It  varies  between  0  and  1,  the  larger  the  value  the  more 
accurate . 

Performing  the  K-S  test  not  only  provides  a  measure  of 
correctness  of  the  model  by  means  of  the  accuracy  but  also 
permits  a  pictorial  evaluation  of  the  behavior  of  the  model. 
If  the  percentiles  of  the  sample  CDF  of  Xj^  are  plotted  against 
the  percentiles  of  a  uniform  distribution,  the  discrepancy 
between  both  becomes  apparent  for  all  values  in  the  interval 
[0,1].  The  nature  of  those  discrepancies  characterize  a 
tracker  as  pessimistic  or  optimistic. 

An   optimistic   tracker   constructs   a   "concentrated" 
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probability  map, i.e.,  it  is  optimistic  about  its  own  ability 
to  determine  the  location  of  the  target.  As  a  consequence,  the 
position  of  the  dummy  target  is  more  often  than  not  far  from 
to  the  center  of  the  empirical  distribution,  thus  yielding 
much  more  frequently  larger  values  of  the  containment 
statistic  X    (Figure  4.1). 

A  pessimistic  tracker  yields  a  more  "dispersed" 
probability  map,  hence  the  frequency  of  lower  values  of  the 
containment  statistic  is  larger.  For  a  given  location  of  the 
target,  the  containment  statistic  X  will  be  stochastically 
larger  in  an  optimistic  tracker  than  in  a  pessimistic  one.  In 
Figure  4.1,  assume  that  both  the  pessimistic  and  the 
optimistic  trackers  represent  the  same  situation.  The  value  of 
the  cumulative  probability  or  "probability  mass"  inside  any 
given  equi-area  contour  in  the  optimistic  tracker  would  be 
greater  than  that  in  the  pessimistic  one.  As  a  consequence, 
if  a  contact  were  produced,  for  instance,  at  any  point  in  such 
a  contour, e .g. ,  the  point  indicated  with  a  chevron  in  both 
maps,  then  the  containment  statistic  X  would  tend  to  be 
greater  in  the  optimistic  tracker  than  in  the  pessimistic  one 
(0.9  and  0.2  respectively).  Such  a  tendency  translates  into 
the  shape  of  the  relative  cumulative  frequency  distribution  of 
X.  The  example  in  Figure  4.1.  shows  that  a  given  intermediate 
value  of  X  ,  for  instance  0.6,  is  much  more  frequent  in  a 
pessimistic  tracker,  i  .e .  ,  82%  of  the  time,  than  in  an 
optimistic  one, i.e.,  18%  of  the  time. 
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OPTIMISIIC  MODEL 


X-0.9 


FESIMISTIC  MODEL 


X-0.2 


"Pmbftbililr  Miq)* 


*ProlMbmiy  Map* 


Figure    4.1.    Examples  of  Pessimistic  and  Optimistic  Trackers 


2 .        Area   of   Uncertainty    (AOU) 

The  Area  of  Uncertainty    (AOU)    can  be  defined  as  the  area 
of   those   cells    (i,j)    where      P(Gij\I)    >   P(G\I), 


AOU=    -£     a,j 


(4.4) 


{i,j)eK 


where  a^^  is  the  area  of  an  individual  cell  and  K  is  defined 
as  in  4.1.  The  value  of  AOU  provides  a  quantitative  idea  of 
how  spread  the  probability  map  is.  The  smaller  the  AOU  the 
better  since  the  probability  mass  is  more  concentrated,  thus 
leading  to  a  better  localization. 
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It  may  be  argued  that  the  AOU  and  the  accuracy  signify 
the  same  MOE,  but  it  is  clear  that,  for  example,  a  constant 
value  of  accuracy  may  coexist  with  different  values  of  AOU. 

3.    Mean  Missed  Distance (MMD) 

The  MOE  selected  to  verify  the  location  features  of  the 

probability  map  with  respect  to  the  dummy  is  the  Root  Mean 

Square  Missed  Distance,   which  can  be  expressed  using  a 

straightforward  extension  from  the  one  used  in  Stone  et  al . 

(1988,  pp.  14) , 

(4.5) 

N 

RMSMD=^ 


N 


where  Pi=  P(T^\I)  is  the  weight  of  each  simulated  target,  N  is 
the  total  number  of  simulated  targets  and  r^i  is  the  distance 
between  the  dummy  d  and  each  simulated  target  i.  This  MOE 
represents  the  standard  error  in  the  location  of  the  target. 

4 .  Mean  Detection  Probability  (MDP) 

The  Mean  Detection  Probability  is 

DR=-^  (4.6) 

N 

where  C^.  is  the  number  of  times  out  of  N  identical  replica- 
tions of  an  identical  experiment  in  which  a  searcher  obtains 
at   least   one   contact   within   a   period   of   time  T. 

5.  Mean  Time  to  Detection  (MTD) 

The  Time  to  Detection  TD  is  the  period  of  time  from  the 
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beginning  of  a  search,  i.e.,  time  zero  of  the  simulation  when 
the  searcher  enters  the  Local  Area,  to  the  time  at  which  the 
first  detection  is  produced  or,  if  no  contact  is  gained,  the 
end  of  the  experimental  scenario.  Then  the  Mean  Time  to 
Detection  is  the  average  among  N  replications  of  the  experi- 
ment. 

-1  ^ 
MTD=-Y  TD.  (4.7) 

B.    VERIFICATION  OF  THE  MODEL  IMPLEMENTATION  CORRECTNESS 

1.    Test  Structure  and  Experiment  Design 

The  purpose  of  this  experiment  was  to  verify  whether  the 
implementation  correctly  models  a  Bayesian  filtering  process. 
The  MOE  or  rather  the  measure  of  correctness  was  the  Accuracy. 

For  this  experiment  a  dummy  target  was  selected  at  random 
among  the  500  simulated  targets.  The  Accuracy  was  measured 
with  respect  to  the  location  of  the  dummy  target.  The 
detections  were  sampled  based  on  the  five  minute  cumulative 
detection  probability  (cdpg)  ,  corresponding  to  each  pair 
searcher -dummy . 

The  model  was  tested  using  Updating  for  Negative  informa- 
tion only  and  for  Negative  and  Positive  information.  In  the 
former  case  the  simulation  ceased  at  the  first  contact.  The 
Reactive  Target  model  was  tested  in  combination  with  each 
updating  type. 
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six  different  scenarios  or  searcher  patterns  were  played 
as  portrayed  in  Figure  4.2.  In  the  first  two  the  searcher  is 
stationary  resembling  sonobuoys  or  dipping  sonars,  in  the 
second  pair  the  searcher  moves  through  the  datum  area  in  a 
straight  line,  and  in  the  fifth  and  sixth,  one  and  two 
searchers  respectively  proceed  with  a  random  walk  type  of 
search  around  the  datum. 

The  number  of  replications  was  50  and  for  each  one  20 
time  samples  were  taken.  The  K-S  test  requires  a  random 
sample, i.e.,  the  observations  should  be  independent  and 
identically  distributed  (Dudewicz  et  al.,1988,  pp.280).  In 
this  case,  the  1000  time  samples  are  not  independent.  Though 
only  one  sample  could  have  been  taken  from  each  replication, 
all  time  samples  were  retained  for  they  provided  smoothing 
within  a  replication  "batch" . 

The  tests  were  identified  with  a  3  digit  label.  The  first 
digit  in  the  test  label  indicates  whether  Updating  for 
Positive  information  took  place,  1  if  yes,  2  if  not.  The 
second  digit  denotes  the  Target  Speed,  1  for  slow  (2-5  kts) 
and  2  for  fast  (6-12  kts)  the  last  digit  corresponds  to  the 
scenario. 

The  test  structure  is  shown  in  Table  4 . 1  and  the  set  of 
default  values  for  the  variables  used  in  the  experiment  are 
shown  in  Table  4.2. 
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Search  Patterns 


1]  On  Datum  Top 


2|  On  Datum 
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5]  One  Random  Searcher 


I  '5f 
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U 

6]  Two  random  searchers 


Figure  4.2.  The  search  patterns  utilized  in  the  tests. 
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SLOW  TARGET 

(2-5  KTS) 

NEGATIVE  INFORMATION 

11 

UPDATING  ONLY 

12 

13 
14 
15 
16 

POSITIVE  AND  NEGATIVE 

11 

INFORMATION  UPDATING 

12 
13 
14 
15 
16 

xai3J.e  t.J.,  'me  experiment  btruccure  ror  testing  tne  correctness  in  the 
implementation  of  the  Bayesian  Updating.  In  all  Tests  the  Target  Model  and 
the  Dummy  are  reactive . 
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VARIABLES 

VALUE 

TARGET  VARIABLES 

Datum  error 

2  nm 

Course  Range 

0-360° 

Minimum  Estimated  Target  Speed 

3  kts 

Maximum  Estimated  Target  Speed 

7  kts 

Reaction  Range 

5  nm 

Reaction  Tack 

150° 

SEARCHER  AND  SONAR  VARIABLES 

Source  Level 

200  Db  re  /^Pa 

Ambient  Noise  Level 

80  Db  re  /xPa 

Detection  Threshold 

20  Db  re  /xPa 

Directivity  Index 

30  Db  re  /iPa 

Frequency- 

5  Khz 

Pulse  Length 

0.05  sec 

Beam  Width 

12° 

Circular  Error  in  Contact  Reports 

1  nm 

Overall  Standard  Error  in  Sonar  Terms. 

6  Db  re  /iPa 

Credibility 

1 

ENVIRONMENT  VARIABLES 

Type  of  Bottom 

Flat , Sandy 

Sound  Speed  Profile 

Isothermal 

Depth 

100  m 

Wind 

10  kts 

TIME  VARIABLES 

Total  Simulated  Search  Time 

100  min 

Scouts  Delay 

15  min 

Table  4.2.   The  default  values  of  variables  used  in  the  Experiment  to 
test  the  correctness  of  the  model . 
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The  correctness  in  the  implementation  of  the  model  for 
each  scenario  was  decided  using  a  Kolmorgorv-Smirnov  One 
Sample  Test.  The  null  hypothesis  regarding  the  sample  CDF  of 
the  containment  test  statistic  Xj^, 

H^iF{x)=x  (4.8) 

tested  for  all  x   against, 

H^iF{x)*x  ^^'^^ 

is  rejected  at  level  a   if 

Dn>d,,a  (4.10) 

where  D„  is  calculated  with  (4.2)  and  d^Q,  is  the  asymptotic 
critical  point  corresponding  to  the  sample  size  n,  obtained 
from  Dudewicz  et  al .  (1988, pp. 670) . 

The  accuracy  for  each  test  was  calculated  and  plotted  for 
a  visualization  of  the  pessimistic  or  optimistic  nature  of  the 
resulting  probability  map  and  then  logged  in  a  table. 

2 .    Results 

Table  4.3  displays  the  results  of  this  test.  The  first 
two  columns  describe  the  tests  regarding  the  type  of  updating 
that  was  allowed,  the  target  model  in  use  and  the  type  of 
scenario  played. 

The  third  column  is  the  sample  size  n  used  in  performing 
the  tests.  The  variations  are  due  to  the  fact  that  in  the  case 
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of  updating  for  negative  information  only,  the  simulation 
stopped  at  the  first  contact .  This  happened  in  a  random 
fashion  in  time  hence  the  number  of  collected  samples  for  each 
replication  was  random  too. 

The  fourth  column  shows  the  average  nature  of  Optimistic, 
Pessimistic  or  Variable  in  the  probability  map  produced  at 
every  time  step  in  accordance  with  Figure  4.1. 


ID 
# 

UPDATING 
NEG   POS 

SAMPLE 
SIZE  n 

OPT 
PES 

KS 
TEST 

ACCURACY 

11 

Y     N 

173 

VAR 

arr 

0.8919 

12 

Y     N 

879 

OPT 

rrr 

0.9343 

13 

Y     N 

671 

VAR 

aaa 

0.9625 

14 

Y     N 

918 

VAR 

arr 

0.9538 

15 

Y     N 

189 

VAR 

arr 

0.9052 

16 

Y     N 

98 

PES 

aar 

0.8767 

21 

Y     Y 

1000 

VAR 

rrr 

0.9343 

22 

Y     Y 

1000 

OPT 

arr 

0.9269 

23 

Y     Y 

1000 

VAR 

rrr 

0.9684 

24 

Y     Y 

1000 

VAR 

aaa 

0.9842 

25 

Y     Y 

1000 

VAR 

rrr 

0.9243 

26 

Y     Y 

1000 

VAR 

aar 

0.9123 

Table  4.3.   The  implementation  correctness  test  results 


The  fifth  column  in  Table  4.3  shows  the  result  of  the 
Kolmorgorof -Smirnoff  test  of  hypothesis  defined  in  the 
previous  section.  The  three  letters  indicate  the  acceptance  of 
rejection  of  the  null  hypothesis  for  three  different  levels  : 
0.01,  0.05  and  0.10.  The  null  hypothesis  is  accepted  most  of 


71 


the  time  for  the  first  level  and  rejected  by  a  slim  margin, 
0.0152  in  the  average,  most  of  the  time  for  the  other  two. 

The  last  column  shows  the  final  values  of  Accuracy- 
obtained  for  each  test.  The  results  range  from  0.8214  to 
0.9842  with  an  average  of  0.8945. 

C.    VERIFICATION  OF  THE  IMPLEMENTATION  ROBUSTNESS 

1.  Test  Structure  and  Experiment  Design 

A  sensitivity  analysis  was  perfoinned  on  Accuracy  as  a 
function  of  Target  Speed,  Searcher  Speed,  Datum  Age,  Source 
Level,  Sample  Size  and  Credibility.  The  purpose  was  to  have  an 
idea  of  the  range  within  which  the  model  remained  correct .  For 
this  experiment,  with  the  exception  of  the  variable  of 
interest,  the  conditions  in  test  25  (Table  4.2)  were  used  for 
all  tests. 

2 .  Results 

Figure  4.3.  shows  the  plots  of  Accuracy  as  a  function  of 
Target  Speed,  Searcher  Speed,  Datum  Age  and  Source  Level.  No 
conspicuous  dependency  can  be  observed  within  the  explored 
ranges  of  each  one  of  those  variables,  and  accuracy  remains 
high. 

An  increase  in  Accuracy  as  a  function  of  Sample  Size  can 
be  observed  in  Figure  4.4.  This  relationship  conveys  the  idea 
that  the  quality  of  the  tests  is  noticeably  affected  by  the 
sample  size. 
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Figure  4.3.  Accuracy  as  a  function  of  Target  Speed,  Searcher  Speed, 
Datum  Age,  and  Source  Level. 
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Figure  4.4.  The  quality  of  the  test  is  sensible  to  the  sample  size 
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D.    THE  INFLUENCE  OF  CREDIBILITY 

1.  Test  Structure  and  Experiment  Design 

A  sensitivity  analysis  was  also  performed  on  Accuracy  as 
a  function  of  Credibility  since  the  latter  is  one  of  the 
hardest  parameters  to  estimate.  The  purpose  was  to  have  an 
idea  of  the  range  within  which  the  model  remained  correct .  In 
this  experiment  the  conditions  in  test  25  (Table  4.2)  were 
used,  except  for  the  credibility  Q.  In  the  tests  performed  in 
Section  B,  credibility  was  set  to  be  equal  to  1 .  In  RTCAS,  the 
user  is  allowed  to  input  other  values.  This  experiment  tested 
values  of  Q  between  0.6  and  1.  Since  in  this  test,  as  in 
Sections  B  and  C,  the  contacts  were  generated  only  at  the 
position  of  the  dummy,  i.e.,  there  were  not  false  alarms,  the 
true  value  of  credibility  in  the  simulation  was  always  1. 

2 .  Results 

As  can  be  observed  in  Figure  4.5,  the  Accuracy  curves 
were  plotted  for  several  values  of  Credibility  larger  than  0.5 
since  smaller  values  are  very  unlikely  to  be  used.  When  Q 
decreases  the  Accuracy  decreases  as  well.  A  tendency  towards 
a  pesimistic  probability  map  output  can  be  seen  as  decreasing 
values  of  Q  causes  the  model  to  ignore  good  information  and 
produce  a  less  "concentrated"  probability  map. 
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Figure  4.5.  Decreasing  the  value  of  Credibility  makes  the  probability- 
map  more  pessimistic . 

E.    THE  SIGNIFICANCE  OF  A  REACTIVE  TARGET  MODEL 


1.    Test  Structure  and  Experiment  Design 

The  purpose  of  this  experiment  was  to  test  whether  or  not 
a  Reactive  Target  Model  is  operationally  better  than  a  Non 
Reactive  one.  All  of  the  MOE's  defined  in  section  A  were  used 
in  order  to  compare  the  results  from  both  Target  Models. 

The  dummy  was  selected  among  the  five  hundred  simulated 
targets  by  means  of  a  unifoirm  variate  drawn  and  was  set  to  be 
always  reactive.  This  implies  the  assumption  that  a  real 
target  is  reactive.  All  MOE's  were  measured  with  respect  to 
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the  location  of  the  dummy  target.  The  detections  were  decided 
as  in  Section  B.  The  model  always  updated  for  negative 
information  only,  stopping  at  the  first  contact.  The  purpose 
of  not  updating  for  positive  information  is  that  a  contact 
report  will  noticeably  improve  most  of  the  MOE's  due  to  the 
"concentration"  of  the  probability  mass  in  the  surroundings  of 
the  target  thus  obscuring  the  difference  between  a  Reactive 
and  a  Non  Reactive  Target  Model.  Furthermore,  the  differences 
between  both  Target  models  are  more  pivotal  to  the  problem 
before  the  first  detection. 

The  Reactive  Target  Model  was  tested  in  the  same 
conditions  as  the  Non  Reactive  Target  Model.  Two  different 
ranges  of  target  speed  and  three  different  values  of  datum 
delay  were  tested  as  can  be  observed  in  the  test  structure 
shown  in  Table  4.4.  The  main  purpose  of  the  test  was  to 
compare  how  good  each  model  (reactive  or  non  reactive)  deals 
with  the  same  scenario  when  constructing  the  probability  map. 
As  a  consequence  variables  with  the  greatest  impact  in  the 
kinematics  of  the  problem  were  selected. 

The  number  of  replications  was  the  same  as  in  the 
experiment  in  Section  D.  The  tests  were  labeled  with  a  3  digit 
code.  The  first  digit  in  this  case  indicates  the  speed  range, 
the  second  the  delay  and  the  last  one  the  scenario  or  search 
pattern,  which  are  the  same  as  in  the  test  in  Sections  D  and 
E  with  the  addition  of  a  seventh  pattern  representing  the 
first  270  degrees  in  an  exhaustive  spiral  search.  The  set  of 
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default  values  for  the  variables  used  in  the  experiment 
remained  the  same  as  in  Table  4.2  with  the  exception  of  the 
target  speed  and  delay  time  which  were  varied  from  test  to 
test.  For  each  running  all  MOE's  were  recorded  in  separate 
files  for  the  reactive  and  non  reactive  case  and  then  compared 
by  means  of  a  Wilcoxon  Test  for  paired  data  (Stone  et  al .  , 
1991,pp.28) .  The  Wilcoxon  Rank-Sum  Test  (Devore, 1991, pp . 612) 
is  a  non  parametric  test  used  to  compare  paired  observations 
without   regard  of   the  distributions  they  come   from. 


Target  Speed/ 
Delay 

2-5  kts 

5-11  kts 

0  min 

111 
112 
113 
114 
115 
116 

211 
212 
213 
214 
215 
216 

15  min 

121 
122 
123 
124 
125 
126 
127 

221 
222 
223 
224 
225 
226 
227 

3  0  min 

131 
132 
133 
134 
135 
136 
137 

231 
232 
233 
234 
235 
236 
237 

Table    4.4.     Test   structure   for   comparing   the  Reactive 
Target  Model   with  the  Non  Reactive   one.    Each  one   of 
the   tests   in  the   table  was  performed  for  each  model , 


77 


Let  i?i  and  N^^  be  the  corresponding  MOE's  for  the  case  of 
a  Reactive  and  a  Non  Reactive  Target  Model  respectively  for 
each  run  i=l.  .n.  We  want  to  compare  whether  or  not  one  of  them 
has  a  tendency  to  be  larger  than  the  other.  R^  and  Nj^  are 
correlated  for  each  i,  but  the  differences  Dj^=  R^-Nj^  are 
independent  for  all  i;  the  Wilcoxon  test  makes  their  original 
distributions  irrelevant  to  the  comparison.  A  statistically 
valid  statement  about  such  a  comparison  can  be  done  by  means 
of  a  test  of  hypothesis  of  the  mean  ii^  and  /ij,  of  each  r.v. 
where 

H^:{n^-^l^=o}  (4.11) 

and, 

H^:{\ij^-\Xj^>0}  (4.12) 

are  the  null  and  alternate  hypotheses.  The  order  in  the 
differences  between  both  mean  values  can  be  varied  depending 
on  the  nature  of  the  MOE.  For  instance  if  R^  and  N^ 
represented  the  Accuracy,  we  would  want  to  know  which  is 
larger  since  the  larger  the  better  in  such  a  case.  Conversely 
if  they  were  the  Area  of  Uncertainty,  we  would  want  to  find 
out  which  is  smaller. 

In  this  test,  the  differences  D^^  =  (R^  -  N^)  are  obtained 
and  sorted  by  ascending  absolute  value .  The  rank  of  such  a 
sorting  is  r^.  Adding  up  the  rank  of  all  negative  D^  yields  the 
test  statistic 
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^=Yl  ^i  (4.13) 

Di>0 

If  there  is  not  a  tendency  of  the  random  variables  R  to 
be  smaller  than  N,  then  W  is  expected  to  be  roughly  close  to 
the  value  of  its  mean  jUj^  =  n  (n+1)  /4,  i.e.,  half  of  the  sum  of 
the  ranks  from  1  to  n.  Conversely,  if  R^  tends  to  be  smaller, 
the  test  statistic  W  will  be  smaller  than  n(n+l)/4. 

For  n>30  W  is  a  normal  random  variable  with  mean  fi^  and 
variance  Oy/  =  (n+1)  (2n+l) /24 .  As  a  consequence  the  test  of 
hypothesis  can  compare  the  value  of  W  with  the  critical  value 
of  a  normal  distribution  C(a,n)  corresponding  to  the  sample 
size  22  and  the  desired  significance  level  a.  If  W  >  C(oL,n) 
then  we  reject  the  null  hypothesis. 

2.    Results 

Table  4.5.  sumarizes  the  results  of  the  Wilcoxon  Tests  of 
Hypothesis  conducted  to  verify  whether  or  not  using  a  Reactive 
Target  Model  improves  the  performance  of  the  tracker  from  the 
standpoint  of  Accuracy,  Area  of  Uncertainty  and  Mean  Missed 
Distance.  The  tests  were  performed  on  the  data  shown  in  Tables 
4.6  and  4.7  ,  where  for  each  MOE  and  each  run,  an  asterisk 
indicates  when  the  Non  Reactive  Target  Model  wins  for  a  visual 
appraisal  of  the  tendencies. 

On  the  average,  as  can  be  observed  in  Table  4.5,  second 
column.  Accuracy,  Area  of  Uncertainty  and  Mean  Missed  Distance 
show  a  better  result  for  a  Reactive  Target  Model.  However, the 
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Wilcoxon  Test  does  not  reject  the  null  hypothesis  that  there 
is  no  difference  between  the  RT  and  the  NRT  models. 
Furthermore,  the  null  hypothesis  should  be  accepted  for 
Accuracy  and  Area  of  Uncertainty.  This  result  is  somewhat 
disappointing  for  we  expected  the  former  to  be  the  MOE  where 
the  improvement  would  be  the  greatest.  Nonetheless,  at  a  10% 
significance,  H^   would  be  rejected. 

The  RT  model  proved  to  be  definitely  better  than  the  NRT 
one  with  respect  to  the  Mean  Missed  Distance,  yet  by  a  slim 
margin . 


MEAN 
NRT 

MEAN 
RT 

W 

Mw 

^w 

p- val- 
ue 

Ho 

TEST 
RESULT 
Oi=    1% 

ACC 

0.8292 

0.8567 

285 

451 

152 

0.2743 

/^N-/^R=0 

Accep- 
ted 

AOU 

9.8121 

8.3223 

187 

451 

152 

0.0834 

/^r-Mn=0 

Accep- 
ted 

MMD 

7.0542 

5.8874 

50 

451 

152 

0.0084 

/iR-MN=0 

Rejec- 
ted 

Table  4.5.  Significance  of  a  Reactive  Target  Model  in  the  average  and  in 
a  Wilcoxon  Test  of  Signed  Ranks.  W  is  the  Wilcoxon  Test  Statistic. 

F.    FXJRTHER  EVALUATION  OF  THE  SIGNIFICANCE  OF  A  REACTIVE 
TARGET  MODEL:  A  NON  OPTIMAL  MYOPIC  SEARCH  EXPERIMENT 


1.    Test  Structure  and  Experiment  Design 

Although  the  experiment  described  in  Section  E  ascertains 
that  a  Reactive  Target  Model  performed  better  than  a  Non 
Reactive  one  in  most  of  the  cases,  the  difference  in 
effectivness  reflected  by  the  tests  is  much  less  significant 
than  the  expectations.  In  the  previous  set  of  tests  the 
searchers  proceeded  without  regard  to  the  probability  map 
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output  by  the  model.  They  followed  typical  lost  contact 
exhaustive  initial  search  patterns. 

In  order  to  further  evaluate  the  significance  of  a 
Reactive  Target  Model  as  compared  to  a  Non  Reactive  Target,  an 
experiment  was  conducted  in  which  the  searchers  were  set  to  be 
fed  back  by  the  output  of  the  model,  conditioning  their 
behavior  by  the  probability  map  rendered  at  every  time  step. 
At  all  times,  the  single  searcher  simply  headed  for  the 
brightest  cell. 

The  visualization  of  the  search  pattern  conducted  with 
the  above  mentioned  rules  was  often  irrational  from  a  tactical 
point  of  view  in  the  sense  that  the  searcher  would,  for 
instance,  move  from  one  brightest  cell  to  the  other  leaving 
unexplored  cells  that  were  close  to  his  way.  Again,  the 
purpose  was  to  have  a  standard  set  of  rules  of  behavior  of  the 
searcher  in  order  to  obtain  more  insight  in  the  comparison 
between  a  Non  Reactive  Target  Model  and  a  Reactive  one.  The 
dummy,  uniformly  drawn  from  the  total  number  of  simulated 
targets  at  each  replication,  was  set  to  be  always  reactive  and 
the  actions  were  set  to  last  until  a  contact  was  obtained  or 
150  minutes,  whichever  occurred  first.  The  detections  were 
decided  as  in  section  B.  The  set  of  tests  was  again  identified 
with  a  3  digit  code  as  shown  in  Table  4.8.  The  first  digit 
indicates  the  target  speed  range,  the  second  one  the  datum  age 
and  the  third  one  the  searcher  speed  and  initial  position. 
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ID  # 

ACC 

AOU 

MMD 

NRT 

RT 

NRT 

RT 

NRT 

RT 

111 

0.7523 

0.8500 

6.6238 

4.4575 

4.2770 

3.4273 

112 

0.9361 

0.9636 

10.0275 

5.9462 

3.5910 

2.5723 

113 

0.9295 

0.9353 

9.3438 

6.5120 

4.4159 

3.9566 

114 

0.9518 

0.9551 

12.2125 

9.5123 

4.6987 

3.6002 

115 

0.7795 

0.8541 

5.4262 

2.5634 

3.4049 

2.7164 

116 

0.7401 

0.8456 

2.8725 

2.1560 

2.8419 

1.2563 

117 

0.9089* 

0.8687 

10.3038 

3.4521 

3.8001 

2.5800 

121 

0.7134 

0.9635 

7.9650 

4.7425 

4.9179 

4.1629 

122 

0.9441* 

0.8242 

13.2612 

10.2040 

4.6058 

4.2763 

123 

0.9230 

0.9261 

10.6025 

9.5623 

6.1347 

5.4798 

124 

0.8035 

0.9012 

17.2188 

12.9650 

6.0323 

5.1757 

125 

0.8872 

0.9162 

10.1212 

8.5612 

4.5867 

5.3936 

126 

0.7731* 

0.6647 

3.8525 

2.3350 

3.6402 

2.8612 

127 

0.8810 

0.8623 

11.618* 

12.2564 

5.3115 

4.7462 

131 

0.6532 

0.7823 

5.1925* 

5.5478 

5.5896 

4.5200 

132 

0.8452* 

0.8124 

9.5412 

8.4150 

5.9528 

3.5600 

133 

0.8687* 

0.8536 

11.2462 

9.9871 

6.6959 

6.4589 

134 

0.8049 

0.9304 

17.3588 

15.2412 

7.2086 

6.0514 

135 

0.7689 

0.8120 

7.4200 

5.8650 

5.6523 

4.0254 

136 

0.8897 

0.7025 

4.6812 

2.7475 

4.6841 

3.8912 

137 

0.8873 

0.9151 

9.430* 

10.2547 

5.9379 

4.2356 

avg 

0.8400 

0.8637 

9.3485 

7.2992 

4.95149 

4.0451 

Table  4.6.  The  results  of  the  Comparison  between  a  Reactive  Target 
Model  (RT)  with  a  Non  Reactive  one  (NRT)  for  slow  target  speeds.  Asterisks 
indicate  when  NRT  wins. 
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ID  # 

ACC 

AOU 

MMD 

NRT 

RT 

NRT 

RT 

NRT 

RT 

211 

0.6947 

0.7942 

7.7325 

5.7474 

6.5314 

4.1247 

212 

0.896* 

0.7123 

13.0422 

5.4373 

6.5159 

2.8195 

213 

0.8624 

0.9237 

9.6524 

7.4512 

6.9158 

6.1266 

214 

0.893* 

0.8222 

12.2575 

8.5647 

7.9634 

4.9579 

215 

0.7242 

0.8520 

7.6672 

7.6542 

5.2720 

3.2549 

216 

0.733* 

0.7047 

3.3073 

2.0707 

4.0580 

3.1092 

217 

0.8962 

0.9014 

13.1071 

8.4475 

5.9064 

3.6117 

221 

0.7530 

0.8562 

12.5350 

9.4358 

10.0125 

8.6542 

222 

0.8444 

0.9115 

8.445* 

8.4512 

8.4972 

5.7633 

223 

0.9134 

0.9139 

13.2425 

9.5623 

9.8595 

9.1776 

224 

0.901* 

0.8425 

13.4175 

12.6275 

9.7686 

8.3196 

225 

0.8223 

0.8366 

11.8777 

7.5158 

8.8233 

7.1634 

226 

0.7219 

0.8675 

4.050* 

5.6475 

7.7192 

6.1558 

227 

0.926* 

0.9027 

12.025* 

12.9775 

9.8413 

6.2208 

231 

0.8125 

0.8457 

19.6851 

10.2125 

13.4794 

11.8224 

232 

0.7846 

0.9315 

6.267* 

8.6958 

10.9819 

10.2591 

233 

0.903* 

0.8728 

8.965* 

9.8561 

12.9581 

12.8511 

234 

0.8044 

0.8841 

9.795* 

15.6056 

12.9458 

12.7729 

235 

0.7337 

0.8412 

6.325* 

12.7345 

9.817* 

11.4833 

236 

0.7951 

0.8324 

15.1725 

7.4125 

12.1102 

12.0041 

237 

0.7711 

0.7958 

6.962* 

20.1525 

12.3220 

11.6747 

AVG 

0.81849 

0.84975 

10.26347 

9.345557 

9.157086 

7.729829 

Table  4.7.  The  results  of  the  Comparison  between  a  Reactive  Target 
Model  (RT)  with  a  Non  Reactive  one  (NRT)  for  fast  target  speeds.  Asterisks 
indicate  when  NRT  wins 
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TEST 
ID 

TARGET 
SPEED 

DATUM 
AGE 

SEARCHER 
SPEED 

111 

2-4  Kts 

10  min 

12  kts 

112 

2-4  Kts 

10  min 

18  kts 

121 

2-4  Kts 

60  min 

12  kts 

122 

2-4  Kts 

60  min 

18  kts 

211 

4-8  Kts 

10  min 

12  kts 

212 

4-8  Kts 

60  min 

18  kts 

221 

4-8  Kts 

10  min 

12  kts 

222 

4-8  Kts 

60  min 

18  kts 

311 

8-12  Kts 

10  min 

12  kts 

312 

8-12  Kts 

60  min 

18  kts 

321 

8-12  Kts 

10  min 

12  kts 

322 

8-12  Kts 

6  0  min 

18  kts 

Table    4.8.    The  non  optimal  myopic  search  experiment   structure. 


For  each  test  3  0  replications  were  run  and  the  average 
values  of  Detection  Rate,  Time  to  Detection  and  Search  Speed 
were   extracted, 

2 .         Results 

Table  4.9  shows  the  results  of  this  test.  The  Mean 
Detection  Probability  showed  the  strongest  difference  between 
a  Non  Reactive  model  and  a  Reactive  one,  favorable  to  the 
latter.  The  Wilcoxon  Test  of  Signed  Ranks  Sum  performed  on 
Mean  Detection  Probability  and  Mean  Detection  Time  yielded  the 
results  shown  in  Table  4.10.  At  a  5  %  significance  level,  the 
null  hypothesis  that  there  is  no  difference  between  a  Reactive 
Target  Model  and  a  Non  Reactive  one  can  be  Rejected  with 


84 


respect  to  the  Mean  Detection  Probabilities  and  should  be 
accepted  regarding  the  Mean  detection  Time. 

The  obseirvation  of  the  evolution  of  the  search  step  by- 
step  on  the  screen  output  of  the  model,  provides  more  insight 
about  the  real  gist  of  Mean  Time  to  Detection  in  the  context 
of  the  test.  In  the  Non  Reactive  Target  Model,  the  searcher 
promptly  starts  "consuming"  probability  mass  since  the 
expanding  ring  of  the  probability  map  does  not  change  because 


MDP 

MDT 

NRT 

RT 

NRT 

RT 

111 

0.700 

0.900* 

73.8 

54.3* 

112 

0.767 

0.933* 

28.1* 

28.3 

121 

0.267 

0.700* 

48.8 

39.6* 

122 

0.867 

1.000* 

55.6 

32.1* 

211 

0.067 

0.167* 

75.5 

67.5* 

212 

0.667 

0.667* 

45.1 

39.2* 

221 

0.167 

0.267* 

45.2* 

46.1 

222 

0.833 

0.900* 

61.3* 

63.3 

311 

0.000 

0.000 

150.0 

150.0 

312 

0.567 

0.633* 

54.2* 

56.7 

321 

0.000 

0.000 

150.0 

150.0 

322 

0.267 

0.367* 

58.6 

39.6* 

AVG 

0.431 

0.486* 

45.5 

38.9* 

Table  4.9.  The  non  optimal  myopic  search  experiment  structure 
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MOE 

W 

Mw 

Ow 

p-value 

Ho 

Hx 

RESULT 

MDP 

0 

39 

15.7 

0 

Mr=Mn 

Mr>Mn 

Rejected 
(5%) 

MDT 

18 

39 

15.7 

0.059 

Mr=/^n 

Mr>Mn 

Accepted 
(5%) 

Table  4.10.  The  Reactive  Target  Model  performed  better  than  the  Non 
Reactive  one  with  respect  to  Mean  Detection  Time  but  no  difference  was 
appreciated  with  respect  to  Mean  Detection  Probabilities. 


of  the  searching  activity.  Therefore,  when  it  proceeds  towards 
the  brightest  cell,  it  tends  to  have  a  lot  of  simulated 
targets  within  its  detection  range.  This  is  reasonable  since 
the  relative  velocity  between  each  simulated  target  and  a 
searcher  will  be  high  because  of  the  average  outbound  nature 
in  the  course  of  the  latter  and  the  inbound  course  in  the 
former  (See  Figure  4.6).  Conversely,  in  the  Reactive  Target 
Model,  the  presence  of  the  searcher  will  "deform"  the 
expanding  ring  of  the  probability  map  in  a  way  such  that  the 
relative  velocity  between  searcher  and  targets  will  be  low  in 
general,  i.e.,  a  chasing  situation  will  be  the  norm  most  of 
the  time.  As  a  consequence,  the  rate  of  consumption  of 
probability  mass  will  be  lower  than  in  the  Non  Reactive  case 
and  the  time  to  detection  will  be  larger,  thus  leading  to  a 
slower  Search  Speed.  (See  Figure  4.7) 

Coalescing  the  ideas  of  the  three  MOE's  analyzed  in  this 
section,  it  can  be  stated  that  using  a  Reactive  Target  Model 
yielded  a  slower  yet  more  effective  search.  Nonetheless  the 
margins  by  which  this  sentence  can  be  stated  are  rather  slim. 
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Figure  4.6.  The  searcher  "consumes"  probability  mass  to  the  N.E.  of 
the  expanding  datum  of  a  Non  Reactive  Target  Model.  The  darker  colors 
indicate  those  cells  where  the  likelihood  is  reduced  by  the  search. 
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"consuming" 


Figure  4.7.  The  reaction  of  the  target  prevents  the  searcher  from 
too  much  probability  mass  as  it  proceeds  with  the  search  pattern.  The  expan- 
sion of  the  datum  is  somewhat  irregular. 
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V.  EXTENSIONS  OF  THE  CONCEPT  AND  IMPLEMENTATION  ENHANCEMENTS 
SUGGESTED  FOR  FUTURE  RESEARCH 


A.    EXTENSIONS  OF  THE  CONCEPT 

The  Information  Fusion  Paradigm  (IFF)  can  be  the  engine 
of  several  problems  other  than  the  one  solved  in  this  thesis. 

1.  Other  Sensors 

Extension  of  the  model  to  other  sensors  like  radar, 
lidar,  infrared,  etc.  is  straightforward.  Even  the  idea  of  a 
multisensor  search  is  plausible  provided  that  the  correspond- 
ing models  to  evaluate  the  effectiveness  of  the  search  are 
implemented.  The  structural  design  of  the  program  allows  the 
interchange  of  the  sonar  detection  model  almost  without 
penalties. 

2.  RTCAS  cuid  Java 

The  appearance  of  the  programming  language  Java  (Sun 
Microsystems  Corporation,  1994)  offers  a  new  dimension  to  the 
Information  Fusion  Paradigm.  Java  allows  relatively  easy 
networking  between  different  types  of  computers.  This  new 
level  of  interchange  of  information  not  only  favors  the 
handling  of  the  real  part  in  the  IFP,  but  also  brings  the 
possibility  of  a  Distributed  Simulation  (Buss,  et  al.,1996), 
as  an  option  to  deal  with  the  "virtual"  part  of  the  IFP.  Such 
a  connectivity  is  a  perfect  solution  when  for  instance  the 
search  is  conducted  with  heterogeneous  assets,  as  in  the  realm 
of  the  Search  and  Rescue  (SAR)  problem.  According  to  interna- 
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tional  regulations,  any  vehicle,  e.g.,  fishing  or  merchant 
vessels,  airplanes,  lorries,  etc.  can  be  convoked  to  partici- 
pate in  a  search  by  a  SAR  agency.  Nowadays,  the  only  existing 
links  between  those  heterogeneous  assets  are  the  International 
Maritime  Radio  Circuits.  Given  that  most  of  those  platforms 
can  carry  a  computer,  Java  brings  the  possibility  of  linking 
them  together.  Then,  it  is  possible  to  implement  a  Real  Time 
Bayesian  Filtering  in  which  each  one  of  the  "convoked" 
searchers  will  not  only  input  its  searching  activity  via,  for 
instance,  an  Internet  browser,  but  also  will  run  its  own 
Stochastically  Synchronized  Simulation,  i.e.,  the  seed  of  the 
random  generation  is  broadcasted  in  the  network.  Hence  each 
one  of  the  searchers  will  have  its  own  picture  of  the  situa- 
tion thus  leading  to  a  more  intelligent  allocation  of  effort. 

B.    IMPLEMENTATION  ENHANCEMENTS 

1.    Model  enhancements 

We  envision  a  lot  of  potential  growth  without  the  penalty 
of  very  expensive  research  and  development .  Some  ideas 
susceptible  to  further  research  are  outlined  in  this  section. 

A  Multiprior  target  generation  scheme  would  enhance  the 
ability  of  RTCAS  to  capture  reality  by  means  of  several 
weighted  options. 

An  Intelligent  Reactive  Target  Model  would  allow  each 
target  to  evaluate  more  factors  before  choosing  a  Reaction 
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Pattern.  One  of  those  factors  may  be  the  counterdetection 
estimated  range.  Solving  the  sonar  equation  for  the  one  way- 
path  from  the  searcher  to  the  target  is  more  realistic  than 
using  a  fixed  distance  as  a  reaction  parameter.  Artificial 
Intelligence  tools  and  concepts  would  be  in  place  when 
developing  the  new  reactive  features. 

In  water  depths  in  the  transition  between  deep  and 
shallow  from  the  acoustic  point  of  view,  shadow  zones  can  be 
expected  within  several  cases  of  sound  velocity  profiles.  As 
a  consequence,  the  search  effort  is  not  uniform  over  the  water 
column.  A  possible  solution  to  this  problem  would  be  enhance 
RTCAS  with  the  addition  of  a  third  dimension  in  the  data 
structure.  The  probability  map  could  be  rendered  for  different 
layers  that  would  help  decide  how  to  allocate  the  search 
effort  not  only  in  time  and  xy  positions  but  in  depth  as 
well. 

2 .    Implementation  Enhancements 

Pascal  was  selected  as  a  programming  language  because  of 
its  robustness  and  maintenance  simplicity.  However,  RTCAS  may 
benefit  from  implementation  in  other  languages  like  C  or  C++ 
mainly  because  of  the  speed  provided  by  the  pointer  arithmetic 
so  useful  in  the  simulation,  or  Java,  because  of  the 
connectivity  that  comes  with  the  language.  The  decision 
should  be  done  only  after  evaluating  advantages  and  disadvan- 
tages of  such  a  change  or  expansion. 
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Regardless  of  the  language,  RTCAS  must  migrate  to  a 
graphic  interactive  user  interface  in  order  to  allow  easier 
manipulation  of  its  possibilities.  The  more  the  interface  of 
the  program  resembles  other  popular  programs,  the  faster  the 
learning  curve  and,  as  a  consequence,  the  efficiency  in  the 
use  of  its  resources.  In  the  realm  of  the  user  interface,  the 
need  can  be  also  mentioned  for  zooming  and  scaling  on  the  fly. 
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VI.  CONCLUSIONS  AND  RECOMMENDATIONS 

A.  CONCLUSIONS 

RTCAS  correctly  implements  a  Bayesian  Filtering  Process 
for  Fusing  Information  from  a  real  world  active  search  with 
a  simulation  encompassing  many  possible  positions  and 
movements  of  a  reactive  target  to  produce  a  probability  map  of 
target  location. 

In  an  active  search,  the  target  will  invariably  react. 
Using  a  very  rudimentary  Reactive  Target  Model  proved  to  be 
better  than  a  non  reactive  one  in  most  of  the  cases. 

B.  RECOMMENDATIONS 

Intensive  testing  of  RTCAS  should  continue  to  further 
verify  its  behavior  and  to  validate  it  for  use  in  real 
environments . 

Extensions  of  the  Information  Fusion  Paradigm  should  be 
explored  and  enhancements  to  RTCAS  should  be  implemented  to 
increase  its  power. 

The  possibility  of  implementing  a  reactive  target  model 
in  currently  used  TDA's  should  be  explored. 
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APPENDIX  A.  RTCAS  MAIN  PROGRAM 


program  rtcas; 

uses 

GRAPH,  CRT,GraphManager,TargtManager,DetecManager, 
CellsManager,Scoutmanager,lnterManager; 

begin 

cirscr; 

Graphmanager. Start; 
Scoutmanager.lnitscouts; 
while  WannaContinue  do  begin 
WannaSeeData:  =  FALSE; 

Scoutmanager.ConnectWithScoutsThruDataLink; 
Intermanager.AskOptions; 
Detecmanager.UpdateSonarParameters; 
if  WannaContinue  then  begin 
Targtmanager.GenTargets; 
Graphmanager. Pleasewait; 
ScoutManager.ResetUndetectedMass; 
if  WannaSeeData  then  begin 

Graphmanager. DisplayData; 
end  else  begin 
for  time:  =  1  to  totalsimtime  do  begin 
Scoutmanager.ReadScoutsDataLink; 
Scoutmanager.UpdatePriorWithScoutslnfo; 
Targtmanager.UpdateWeightsWithBayes; 
Scoutmanager.UpdateUndetectedMass; 
Cellsmanager.UpdatelnCeilProbabilities; 
Graphmanager. OutputGraph; 
end; 
end; 
end; 

Scoutmanager.CloseScoutsDataLink; 
end; 
end. 
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APPENDIX  B.  UNIT  INTERMANAGER 


unit  InterManager; 


interface 
<: 


uses 
CRT; 

{TARGETS  CONSTANTS.  TYPES  AND  VARIABLES} 


const 

NUMBEROFTGTS 

=500; 

BIGNUMBER 

=1E30; 

DEFMOTIONMOO 

='OMNI 

DEFSIGMARINIT 

=20; 

DEFPORTMOSTCOURSE=0; 

DEFCOURSERANGE 

=360; 

DEFMINSUBSPEED 

=3; 

DEFMAXSUBSPEED 

=7; 

DEFREACTRANAGE 

=40; 

DEFREACTIONTACK 

=150; 

OMNI 

='OMNI 

type 

targetype 

=  r( 

X  :  integer; 

y  : integer; 

course      :real; 

speed       :real; 

cellx       : integer; 
celly       : integer; 
prob        :reaL; 
threatdis    :real; 
evas  i  oncourse : r ea  L ; 
datumx      : integer; 
datumy      : integer; 
end; 
tgtbundletype   =  array [1. .NUMBEROFTGTS]  of  targetype; 


{NUMBER  OF  TARGETS> 

CUT  I  LI TY  VARIABLE} 

{OMNIDIRECTIONAL  MOTION  MODEL} 

{DATUM  ERROR} 

{LEFT  LIMIT  COURSE} 

{COURSE  RANGE} 

{MIN  SUB  SPEED} 

{MAX  SUB  SPEED} 

{REACTION  RANGE} 

{REACTION  CHANGE  IN  COURSE} 

{MOTION  MODEL} 


{TARGET  STRUCTURE} 

{X  POSITION} 

{Y  POSITION} 

{TARGET  COURSE} 

{TARGET  SPEED} 

{TARGET  CELL  IN  X} 

{TARGET  CELL  IN  Y} 

{TARGET  INDIVIDUAL  WEIGHT} 

{THREAT  DISTANCE  ACCUMULATOR} 

{TACK  ANGLE} 

{TARGET  MIND  DATUM  POSITION  IN  X} 

{TARGET  MIND  DATUM  POSITION  IN  Y} 

{SET  OF  ALL  TARGETS} 


var 

mytarget 

sigmarinit 

mi  nsubspeed 

maxsubspeed 

courserange 

portmostcourse 

motionmodel 

react range 

react iontack 

totprob 


: tgtbundletype; 

:byte; 

:byte; 

:byte; 

: integer; 

: integer; 

: St ring; 

: integer; 

: integer; 

: rea I ; 


{STANDARD  DEVIATION  OF  THE  DATUM} 

{ESTIMATED  MIN  TARGET  SPEED} 

{ESTIMATED  MAX  TARGET  SPEED} 

{ESTIMATED  RANGE  OF  COURSES  FOR  THE  TARGET} 

{PORT  LIMIT  IN  THE  COURSE  RANGE} 

{TYPE  OF  MOTION  MODEL} 

{DISTANCE  OF  REACTION  FROM  SEARCHER} 

{ANGLE  OF  EVASION  FROM  THE  BEARING  TO  THE  SEARCHER} 

{ACCUMULATOR  OF  TARGETS  WEIGHTS  AFTER  UPDATING} 
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^SEARCHERS  CONSTANTS.  TYPES  AND  VARIABLES} 

r                                                                                                                                                                                                                                                                    -v 

const 

J 

MAXSCOUT    =15; 

CMAX  NUMBER  OF  SCOUTS} 

ANGLE       =360; 

<UTILITY  CONSTANT} 

NUMOFPULSES  =10; 

(NUMBER  OF  SONAR  PULSES  PER  SCAN  PERIOD} 

CDE FAULT  VARIABLE  VALUES} 

DEFSL       =200; 

(SOURCE  LEVEL} 

DEFNL       =80; 

(NOISE  LEVEL} 

DEFDT      =20; 

(DETECTION  THRESHOLD} 

DEFDI       =20; 

(DIRECTIVITY  INDEX} 

DEFABSCOEFF  =0.004; 

(ABSORPTION  COEFFICIENT} 

DEFREVRBAREA=8; 

(REVERBERATION  AREA} 

DEFSCATTSTRG=-40; 

(AREA  SCATTERING  STRENGTH} 

DEFCONT_A0U  =40; 

(AREA  OF  UNCERTAINTY  IN  CONTACT  REPORTS} 

DEFFREQUENCY=6; 

(FREQUENCY} 

DEFPULSELENG=0.05; 

(PULSE  LENGTH} 

DEFBEAMWIDTH=14; 

(BEAMWIDTH} 

DEFNOISESIG  =30; 

(UNCERTAINTY  IN  SONAR  EQN.  TERMS} 

type 

scout type      =  record 

(SEARCHERS  DATA  STRUCTURE} 

X 

: integer;                        (SEARCHER  X  POSITION} 

y 

: integer;                        (SEARCHER  Y  POSITION} 

course 

: integer;                          (SEARCHER  COURSE} 

speed 

:real;                              (SEARCHER  SPEED} 

SONAR_ON   : boo lean; 

(SONAR  CONDITION  ON  OFF} 

maxscale   : integer; 

(SONAR  SCALE  SELECTED  IN  THE  SEARCHER'S  SONAR} 

IS_HOT     : boo lean; 

(WHETHER  THE  SEARCHER  HOLDS  A  CONTACT  OR  NOT} 

abscoeff   rreal; 

(ABSORPTION  COEFFICIENT} 

frequency  :real; 

(SONAR  FREQUENCY} 

pulselength:real; 

(PULSE  LENGTH} 

reverbarea  : integer; 

(AREA  OF  REVERBERATION  FOR  GIVEN  LENGTH  AND  BEAMWIDTH} 

beannwidth  :real; 

(SONAR  BEAMWIDTH} 

SL        : integer; 

(SOURCE  LEVEL} 

NL        :byte; 

(NOISE  LEVEL  CORRESPONDING  TO  THE  SEARCHER  SPEED} 

DT        :byte; 

(DETECTION  THRESHOLD} 

DI        :byte; 

(DIRECTIVITY  INDEX} 

SS        : integer; 

(SCATTERING  STRENGTH  CORRESPONDING  TO  SONAR  CONDITIONS} 

subx      : integer; 

(X  POSITION  OF  A  CONTACT} 

suby      : integer; 

(Y  POSITION  OF  A  CONTACT} 

subcourse  : integer; 

(CONTACT  COURSE} 

subspeed   : integer; 

(CONTACT  SPEED} 

cred      :real; 

(CREDIBILITY  IN  THE  CONTACT  REPORT} 

id        :char; 

(PLATFORM  IDENTITY,  SURFACE  AIR  BUOY} 

initialhot  : integer; 

(TIME  OF  INITIAL  CONTACT} 

end; 

scout teamtype   =  array [1. .MAXSCOUT]  of  scouttype;        ^STRUCTURE  TYPE  OF  ALL  SEARCHERS} 

var 

myscout        :scoutteamtype; 

(INSTANCIATED  SEARCHERS} 

numberof scouts  :byte; 

(TOTAL  NUMBER  OF  INSTANCIATED  SCOUTS} 

noises igma     :byte; 

(STANDARD  DEVIATION  IN  THE  SOLUTION  OF  SONAR  EQQUATION} 

undetectedmass  :real; 

(UNDETECTED  PROBABILITY  MASS} 
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<ENVIRONMENTAL  VARIABLES  AND  CONSTANTS> 
< 


const 

DEFB0TT0MTYP=1; 
DEFSVTPR0FIL=1; 
DEFWIND     =10; 
DEFDEPTH    =200; 


(BOTTOM  TYPE> 

CSOUND  VELOCITY  PROFILE  TYPE> 

{WIND  INTENSITY  IN  KNOTS} 

(DEPTH  IN  M> 


var 

bottomtype  :byte; 

svtprofiletype  :byte; 

wind  : integer; 

depth  : integer; 

(SIMUUTION  AND  CONTROL  FLOW  PARAMETERS} 
< 


(TYPE  OF  BOTTOM  I  II  OR  III> 

(TYPE  OF  VERTICAL  SOUND  SPEED  PROFILE} 

(WIND  INTENSITY  IN  KNOTS> 

(DEPTH  IN  M> 


const 

DEFTOTALSIMTIME  =15; 
DEFSCOUTSDELAY  =  6; 


(TOTAL  SIMULATION  TIME> 
(SCOUTS  DATUM  DELAY> 


var 

UannaContinue  : boo lean; 

WannaSeeData  : boo lean; 

data  I  ink  :text; 
dlinktotaltime  : integer; 

totalsimtime  :word; 

time  :word; 

scoutsdelay  :real; 


(EXIT  BOOLEAN  CONDITION} 

(GRAPHIC  DISPLAY  BOOLEAN  CONDITON} 

(LOCATION  OF  THE  DATA  LINK  DEVICE  DRIVER} 

(TOTAL  SIMULATION  TIME  IN  DATA  LINK  DEVICE  DRIVER} 

(TIME  IN  MINUTES  FOR  WHICH  A  RENDER  IS  DESIRED} 

(DISCRETE  TIME  UNITS  COUNTER} 

(TIME  OFF  OF  THE  SEARCHERS} 


(PUBLIC  METHODS} 
( 


procedure  GetTargetData; 
procedure  GetScoutsData; 
procedure  GetEnvironmentalData; 
procedure  GetSimData; 
procedure  AskOptions; 
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implementation 

uses  detecmanager; 

procedure  PromptOptions; 

var  i  : integer; 

begin 

clrscr; 

textcolor(Cyan); 

clrscr; 


writelnC 

writelnC 

writelnC 

writelnC 

writeln; 

for  i:=1  to  16  do  writeln; 

textcolor(Red); 

writelnC 

writelnC 

writelnC 

writelnC 

writelnC 

writelnC 

writelnC 

writeln; 

write  C 

end; 


Please,  select  one  of  the  following  options:'); 


0.  QUIT  PROGRAM'); 

1.  UPDATE  TARGET  DATA'); 

2.  UPDATE  SEARCHER  DATA'); 

3.  UPDATE  ENVIRONMENTAL  DATA'  ); 
U.  SET  TIME  VARIABLES'); 

5.  START  SIMULATION  WITH  CURRENT  DATA' 

6.  DISPLAY  INITIAL  TARGET  POSITIONS'  ) 

Enter  your  option:  '); 


); 


procedure  AskOptions; 
var 

answer         :char; 
origmode       : integer; 
begin 

answer:='9'; 
or  i  gmode: =LastMode; 
TextMode(C80+Font8x8); 

while  ((answero'5')  and  <answer<>'6'))  and(answer<>'0')  do 
isegin 

PromptOptions; 
answer : =readkey; 
case  answer  of 

'0':  WannaContinue:=FALSE; 
GetTargetData; 
GetScoutsData; 
GetEnvi  ronmentalData; 
GetSimData; 


'1' 

'2' 

'3' 

'4' 

'5' 

'6' 
end; 
end; 

TextMode(ori gmode); 
end; 


WannaSeeData:=TRUE; 
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procedure  GetTargetReaction; 

var 

count  :byte; 

begin 

textcolor(Red); 

clrscr;  

wr  i  te  I  n( '    ^s0^ill^jllpls=ilH^I;l|||pj=p=jBpip|^  ' ) ; 

writelnC    IIIJiii-ks;5i51ilJ^I=g^lijliJl!^  ) ; 

writelnC    Please  input  the  following  data  for  500  targets:'); 

for  count:=1  to  12  do  writeln; 

writeln;writeln; 

writelnC    TARGET  REACTION  PARAMETERS  '); 

writeln; 

writeln; 

writeln; 

writelnC    Please  input  the  estimated  REACTION  RANGE  in  nautical  miles,    '); 

writeln; 

writelnC    The  current  value  is  :  ',reactrange/10:3:1); 

write  C    Enter  your  value  :     '); 

readln(reactrange); 

react range :=reactrange*10; 

writeln; 

writelnC    Please  input  the  estimated  REACTION  TACK  in  nautical  degrees,    '); 

writelnC    (Evasion  course  relative  to  the  bearing  to  the  closest  scout     '); 

writeln; 

writelnC    The  current  value  is  :  ',reactiontack:3); 

write  C    Enter  your  value  :     '); 

readln(reactiontack}; 

end; 


procedure  GetDatutnError: 

van 

count  :byte; 

begi  n 

textcolor(Red); 

clrscr; 

writelnC         l^gp;ipipSijipi|S TARGET  DATA       p|pii:ipii|iiiii ' ) ; 

writelnC    Iglli=jsll^^l=^iljj^il^^i^ing5g!lll^i^^§si^^  ) ; 

writelnC    Please  input  the  following  data  for  500  targets:'); 

for  count :=1  to  12  do  writeln; 

writeln;writeln; 

writelnC    DATUM  OMNIDIRECTIONAL  ERROR  (integer)  '); 

writelnC    (Standard  Deviation  in  Tgts  Position,  in  nm  '); 

writeln; 

writelnC    Current  value  is  :   ',sigmarinit/10:3:0); 

write  C    Enter  your  value  :     '); 

readln(sigmarinit); 

sigmarinit:=sigmarinit*10; 

writeln; 

writeln; 

writeln; 

writeln; 

end; 
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procedure  GetMotionModel ; 

var 

ans   :byte; 

count  :byte; 

begin 

textcolor(Red); 

clrscr; 

writelnC 

writelnC 

writelnC 

writelnC    Please  input  the  following  data  for  500  targets:'); 

for  count:=1  to  12  do  writeln; 

writeln;writeln; 


writelnC 
writeln; 
writelnC 
writelnC 
writeln; 
writelnC 
write  C 
readln(ans); 
case  ans  of 
1:  begin 


TARGET  MOTION  MODEL 

(1)  OMNI  :  Omnidirectional  fleeing 

(2)  FAN  :  Directional  movement'); 

Current  motion  model  is 

Enter  your  model  (1  OMNI,  2  FAN) 


'); 


'); 


',motionmodel); 

'); 


2: 
end; 
end; 


motionmodel  :=  'OMNI' 

portmostcourse:=0; 

courserange:=360; 

end; 

motionmodel  :=  'FAN'; 


procedure  GetlargetCourses; 

var 

count  :byte; 

begin 

textcolor(Red); 

clrscr; 

writelnC    Please  input  the  following  data  for  500  targets:'); 

for  count:=1  to  12  do  writeln; 

writeln;writeln; 

writelnC    TARGET  COURSE  LIMITS  '); 

writeln; 

writeln; 

writeln; 

writelnC    Please  input  the  courses  port  limit         "  '); 

writeln; 

writelnC    The  current  value  is  :  ',portmostcourse:3); 

write  C    Enter  your  value  :     '); 

readln(portmostcourse); 

writeln; 

writelnC    Please  input  the  course  range  '); 

writeln; 

writelnC    The  current  value  is  :   ',courserange:3); 

write  C    Enter  your  value  :     '); 

readln(courserange); 

writeln; 

end: 
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procedure  GetTargetSpeeds: 

var 

count  :byte; 

begin 

textcolor(Red); 

clrscr; 

writelnC    jy±!!i!!rS^='TiiJiiiS!gI^i^^==^^  ' ) ; 

writelnC    Please  input  the  following  data  for  500  targets:'); 

for  count:=1  to  12  do  writeln; 

writeln;writeln; 

writelnC    TARGET  SPEED  LIMITS  '); 

writeln; 

writeln; 

writeln; 

writelnC    Please  input  the  Minimum  estimated  Target  speeed  '); 

writeln; 

writelnC    The  current  value  is  :  ',round(minsubspeed*1.2):3); 

write  C    Enter  your  value  :     '); 

read  I n(mi  nsubspeed) ; 

minsubspeed:=  round(0.8333*minsubspeed);<converting  to  discrete  speed  units> 

writeln; 

writelnC    Please  input  the  Maximum  estimated  Target  speed  '); 

writeln; 

writelnC    The  current  value  is  :  ',round(maxsubspeed*1.2):3); 

write  C    Enter  your  value  :     '); 

readln(maxsubspeed); 

maxsubspeed:=  round(0.8333*maxsubspeed>;<:converting  to  discrete  speed  units> 

writeln; 

end; 


procedure  GetTargetData: 

var 

answering   :  char; 

count      :  byte; 

begin 

answering: ='P'; 

while  answering  <>  '0'  do  begin 

textcolor(Red); 

clrscr; 

writelnC  ^^^^^^^^^^^^^^-^M^^^^^^^^^^^^^^^M^^^^^^^^^M ' ) ; 

writelnC  Default  values  of  the  parameters  used  to  generate  and  update  '); 

writelnC  500  simulated  tracks:  '); 

for  count :=1  to  8  do  writeln; 

writeln; 

writelnC      (0)  Back  to  Main  '); 

writeln; 

write  C      (1)  DATUM  ERROR  IN  NM  :  ',sigmarinit/10:3:0); 

writeln;writeln; 

writelnC       (2)  TARGET  MOTION  MODEL  :   ',motionmodel); 

writelnC  (speed  &  course  uniformly  distributed)'); 

writeln; 

writelnC       (3)  COURSE  LIMITS  (FOR  FAN  MOTION  MODEL  ONLY)        '); 

writeln; 

writelnC  PORT  LIMIT  (0-360)  :   ',portmostcourse:3); 

writelnC  RANGE      (Not  starbord  limit)        :  ',courserange:3); 

writeln; 

writelnC       (4)  SPEED'); 
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wn'teln; 

writelnC  MINIMUM  ESTIMATED  TARGET  SPEED  IN  KTS   :   ',round(minsubspeed*1.2):3); 

writelnC  MAXIMUM  ESTIMATED  TARGET  SPEED  IN  KTS   :   ' ,round(maxsubspeecl*l .2):3); 

w^iteln; 
w^iteln('       (5)  TARGET  REACTION  MODEL  '); 
write In; 

writelnC  REACTION  RANGE  :   ',reactrange/10:3:0); 

writelnC  REACTION  TACK  :  ',reactiontack:3); 

writeln; 

writelnC      Enter  the  nuniber  of  the  parameters  you  want  to  change  or  0  '); 
write  C      to  go  back  to  main.  Return  is  not  necessary  :  '  ); 
answer  i  ng : =readkey ; 
case  answering  of 

'1'   :  GetDatumError; 

'2'   :  GetMotionModel; 

'3'  :  GetTargetCourses; 

'4'  :  GetTargetSpeeds; 

'5'  :  GetTargetReaction; 
end; 

end; <wh  i I e> 
end : {GetTargetData} 


procedure  DisplayDataLinklnfo; 

var 

k       :byte; 

begin 

clrscr; 

write  C    Connected  to  Data  Link,  #  of  scouts  present    :  '); 

textcolor(140); 

writelnC    — >  '.numberof scouts); 

textcolor(red>; 

write  C  Time  searching  so  far  :  '); 

textcolor(140); 

writelnC    — >  ',dlinktotaltime*5,'  min'); 

textcolor(red); 

writeln; 

writelnC    Receiving  the  following  initial  conditions:'); 

writeln; 

writeln; 

writelnC     SCOUT  #    XPOS   YPOS    COURSE   SPEED   SONAR_ON   SCALE        '); 

writeln; 

for  k:=1  to  numberofscouts  do  begin 

with  myscoutCk]  do  begin 

writelnC    ',k:8  ,'  ',  x:8  ,  y:8,  course:8,  speed:8:0,SONAR_ON:8,maxscale*185:8); 

end; 
end; 

writeln; 
writeln; 
writelnC      Hit  any  key  to  continue...'); 


readkey; 
end: 
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procedure  ChangeSonarDefaults(scout.param:byte) ; 

var 

f       rreal; 

k       :byte; 

begi  n 

clrscr; 

uriteinc   5~;.;  :::ii:i;^te^^?:^::^::^i=====^^!M^  ' ) : 

writelnC    Jj^^j^siJJSgj^giiPp  ""sCOUTS  "  t^       J^IS^^^i^^ig'); 

writelnC    Scouts  Sonar  Parameters:'); 
wnteln; 
write In; 

writelnC    )  H  '); 

writelnC  j  (1)  (2)  (3)  (4)  (5)  (6)  |  '); 
writelnC  |SCOUT  SL  DT  DI  Frq.  PulseL.  Beam  W.  f  '); 
writelnC  |  #  (dB)  (dB)  (dB)  (Khz)  (msec)  (deg)  |  '); 
writelnC    {  }  '); 

for  k:=1  to  niJii)erofscouts  do  begin 
with  myscout[k]  do  begin 

writelnC    |',k:5  ,'   ',SL:5,'   '  ,DT:5,'   ',DI:5,'  ',frequency:8:1,' 
,pulselength*1000:8:0,'  ',beamwidth:8:1,'    |'); 
end; 
end; 

writelnC    H  »'); 

writeln; 

writelnC    You  are  going  to  produce  modifications  to  scout  number:  '  .scout); 
with  myscout [scout]  do  begin 
if  param  =  1  then  begin 
writeln; 

writeC    Please  input  new  SOURCE  LEVEL  :  '); 

readln(SL); 
end; 

if  param  =  2  then  begin 
writeln; 

writeC    Please  input  new  DETECTION  THRESHOLD  :  '); 

readln(DT); 
end; 

if  param  =  3  then  begin 
writeln; 

writeC    Please  input  new  DIRECTIVITY  INDEX  :'); 

readln(DI); 
end; 

if  param  =  4  then  begin 
writeln; 

writeC    Please  input  new  FREQUENCY(updates  aborption  coeff)   :  '); 
read I n( frequency); 

abscoeff:=AbsorptionCoefficient (frequency); 
end; 

if  param  =  5  then  begin 
writeln; 

writeC    Please  input  new  PULSE  LENGTH(in  ms)  :  '); 

readln(pulseLength); 
pulselength:=pulselength/1000; 
end; 

if  param  =  6  then  begin 
writeln; 

writeC    Please  input  new  BEAM  WIDTH(in  degrees)  :  '); 

readln(beamwidth); 
end; 
end; 
end: 
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procedure  DisplayandChangeSonarDefaults; 
var 

k,item  :byte; 
scoutmod  :byte; 
option  :char; 
begin 

while  scoutmod  <>  0  do  begin 
clrscr; 


writelnC 
writelnC 
writelnC ' 

■■■:=•■=====•■: 

DATA       ^~|p=^|^=s 

If); 

:r=HHH::?=:Hr 

writelnC 

Scouts 

Sonar  Parameters:'); 

writeln; 

writeln; 

writelnC 
writelnC 

/  \ 

(1)     (2)     (3) 

(4)     (5)       (6) 

1 ') 

writelnC 

SCOUT 

SL       DT      DI 

Frq.    PulseL.    Beam  U. 

1 ') 

writelnC 

# 

(dB)    (dB)    (dB) 

(Khz)   (msec)    (deg) 

1 ') 

writelnC 

-^     M 

~\        ' 

for  k:=1  to  numberof scouts  do  begin 
with  myscout[k]  do  begin 
writelnC    |',k:5  ,'   ',SL:5,'   '  ,DT:5,'   ',DI:5,'  ',frequency:8:1,'  ' 
,pulselength*1000:8:0,'  ',beamwidth:8:1,'    |'); 
end; 
end; 

writelnC    t »'); 

writeln; 

writelnC    Please  input  the  number  of  the  searcher  you  want  to  modify  '); 

write  C    or  0  to  go  back  to  Scouts  Menu  :  '); 

read I n( scoutmod ) ; 

writeln; 

if  scoutmod  <>  0  then  begin 

if  scoutmod  >  numberof scouts  then  begin 

writelnCPlease  input  scout  #  smaller  than  '.numberof scouts); 
writelnC Hit  any  key  to  continue..'); 
readkey; 
end  else  begin 

write  ('     Please  input  the  coluim  of  the  item  you  want  to  modify    :  '); 
read I n(i tern); 
writeln; 

ChangeSonarDefaults(scoutmod,item); 
end; 
end;<:if> 
end; {while} 
end; 


procedure  Di  spl ayandChangeErrorDef aul ts ; 

var 

th  i  sanswer  : i  nteger; 


begin 
clrscr; 
writelnC 


wr  i  tein( '       Hiiiiiiipiiyjiji|piiiiiiiiiiiiijii_    scouts     data 

writelnC    ^^^^^^^^^^^^^^^^^^^^^^^^^M^ 


writelnC  ERROR  VARIABLE  DEFAULT  VALUE 

writeln; 

writeln; 

writeln; 

writelnC    Total  estimated  error  in  the  calculation  of  the  terms  in  the     '); 

writelnC    calculation  of  the  sonar  equation  and  in  the  uncertainties  in    '); 
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the  detection  process  in  dB  re  /vPascal, 


"NOISESIGMA"     current  value  is 


writelnC 

'); 

writeln; 
writelnC 
writeln; 

write  ('    Please  input  the  new  valaue 
readln(noisesigma>; 
end:  {procedure} 


',noisesignia:6  ,'  dB  re  /vPascal'); 
'); 


procedure  GetScoutsData: 
van 

response     :char; 
begin 

response: ='9'; 

while  response  <>  'O'do  begin 
clrscr; 


writeln 
writeln 
writeln 
writeln 
writeln 
writeln 
writeln 
writeln 
writeln 
writeln 
writeln 
writeln 
writeln 
writeln 
writeln 
wr  i  te 


Please  select  one  of  the  following  options: 


'); 


0.  BACK  TO  MAIN  '},' 

1.  DISPLAY  INFO  RECEIVED  FROM  SCOUTS      '); 

2.  DISPLAY  AND  OR  CHANGE  DEFAULT  SONAR  PARAMETERS'); 

3.  DISPLAY  AND  OR  CHANGE  DEFAULT  ERROR  VALUES  '); 

Enter  your  option:  '); 


response: =readkey; 
case  response  of 

'1':  DisplayOataLinklnfo; 

DisplayandChangeSonarDefaults; 

DisplayandChangeErrorDefaults; 


'2 
'3 

'0 
end; 
end; 
end: 


procedure  GetSimData; 

var 

answer: char; 

begin 

clrscr; 

writelnC         Izii^in 

writelnC 

writelnC 

writelnC 

writelnC 

writelnC 

writelnC 

writeln; 


SIMULATION 


TIME 


The  simulation  will  be  carried  out  from  time  0  to  TotaiSimt 
in  updating  steps  of  5  minutes.  Totalsimtime  should  be  smal 
than  the  last  updating  time  to  the  data- link  device  file. 


5-  =•:;'); 
,-  _  '); 
"  "  '); 

ime  '); 
ler  '); 

'); 

'); 
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write In; 

write In; 

writelnC    Total  Simulation  Time  (TotalSimtime,  a  multiple  of  5  min)'); 

write  ('    No  info  from  scouts  available  after         :  '); 

textcolor(140); 

writeln(dlinktotaltime*5,'min'); 

textcolor(red>; 

writeln; 

writelnC    Current  value  is  :  ',totalsimtime*5,'  min'); 

write  ('    Enter  your  value  :  '); 

readln(totalsimtime); 

totalsimtime:=totalsimtime  div  5; 

if  totalsimtime  >  dlinktotaltime  then  begin 

writelnC    Simulation  time  should  be  smaller  than  ',dlinktotaltime*5); 

writeln; writeln; writeln; writeln; writeln; 

writelnC    Hit  any  key  to  continue '); 

readkey; 
end; 

writeln; 
writeln; 

writelnC    Delay  Time  of  scouts  after  DATUM  time, (DATUM  "age"  in  min)  '); 
writelnC    Current  value  is  :  ',scoutsdelay*5:4:0); 

write  C    Enter  your  value  :  '); 

read I n( scoutsde I ay) ; 
scoutsdelay:=scoutsdelay/5; 
end; 


function  Profileinwords: string; 

var 

a:string; 

begin 

case  svtprof iletype  of 

1:  a:='POSITIVE  OR  ISOSPEED'; 
2:  a:='NEGATIVE'; 

end; 
Prof i leinwords:=a; 
end; 


function  Bottominwords: string; 
begin 

case  bottomtype  of 


end; 


1 
2 
3 
end; 


Bottominwords:='SANDY,  FLAT'; 

Bottominwords:=' INTERMEDIATE  BETWEEN  VERY  ROUGH  AND  FLAT' 

Bottominwords:=' ROCKY,  VERY  ROUGH'; 


procedure  Getwindspeed; 

begin 

clrscr; 

writelnC         HilSBIIiBIII 

writelnC         Pl=ni~=;i^^i 

writelnC        llJOJ^IIJJJJjJllHl 

writelnC 

writeln; 

writeln; 

writeln; 

writelnC 

writelnC 

writelnC 

writelnC 


ENVIRONMENTAL   DATA 


This  is  the  current  environmental  data: 


(1)  WINDSPEED 

(2)  SOUND  SPEED  PROOF  I LE 

(3)  BOTTOM  TYPE 

(4)  DEPTH 


',wind:4,'  Knots' 
', Profileinwords 
', Bottominwords 

',depth:4,'  m' 


111') 

B'> 

10') 

') 


); 
>; 
); 
); 
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writeln; 

write  ('    Please  input  wind  speed  (knots): 

readln(wind); 

end; 


); 


This  is  the  current  environmental  data: 


procedure  GetSoundProfile; 

var 

OUTTAHERE  : boo lean; 

begin 

CXJTTAHERE:=FALSE; 

while  not  OUTTAHERE  do  begin 

clrscr; 

writeln 

writeln 

writeln 

writeln 

writeln 

writeln 

writeln 

writeln 

writeln 

writeln 

writeln 

writeln 

writeln 

writeln 

writeln 

writeln 

writeln 

writeln 

write  ('    Please  enter  your  option  here  :   ' 

readln(svtprof i letype); 

if  (svtprof iletype=  1)  or  (svtprofi letype  =  2)  then  begin 

OUTTAHERE :=TRUE ; 
end; 
end: 


I'); 


(1)  UINDSPEED 

(2)  SOUND  SPEED  PROOFILE 

(3)  BOTTOM  TYPE 

(4)  DEPTH 


',wind:4,'  Knots' 
',Prof ileinwords 
',Bottominwords 

',depth:4,'  m' 


); 
); 
>; 
); 


You  can  choose  the  sound  speed  profile  according  to  the  following'); 
code:'); 


(1)  POSITIVE  OR  ISOSPEED 

(2)  NEGATIVE 


'); 


); 


procedure 

var 

OUTTAHERE 

begin 

OUTTAHERE; 

wh  i  I  e  not 

clrscr; 

writeln! 

writeln! 

writeln! 

writeln! 

'); 

writeln; 
writeln; 
writeln! 
writeln! 
writeln! 
writeln! 
writeln; 
writeln; 
writeln! 
writeln; 
writeln! 
writeln! 
writelni 


GetBottomtype: 

: boo lean; 

=FALSE; 

OUTTAHERE  do  begin 


ENVIRONMENTAL   DATA 


This  is  the  current  bottom  &  depth  environmental  data: 


(1)  WINDSPEED 

(2)  SOUND  SPEED  PROOFILE 

(3)  BOTTOM  TYPE 

(4)  DEPTH 


',wind:4,'  Knots' 
'.Prof ileinwords 
' .Bottominwords 

', depth :4,'  m' 


You  can  input  bottom  type  according  to  the  following  code'); 

(1)  SANDY, FLAT  BOTTOM, NEGLECTABLE  ROUGHNEESS 

(2)  INTERMEDIATE  BETWEEN  FLAT  AND  VERY  ROUGH 

(3)  ROCKY,  VERY  ROUGH 


); 
I'); 


); 
); 
); 
>; 


'); 
'); 
'); 


109 


writeln; 

write  ('    Please  enter  your  option  here  :   ' 

read I n( bot tomtype ) ; 

if  (bottomtype  =  1)  or  (bottotntype  =  2)  or  (bottomtype  =  3)  then  begin 

OUTTAHERE:=TRUE; 
end; 
end; 
end: 


); 


procedure  Getdepth: 

begin 

clrscr; 

writelnC 

writelnC 

writelnC 

writelnC 

writeln; 

writeln; 

writeln; 

writelnC 

writelnC 

writelnC 

writelnC 

writeln; 

write  C 

read I n( depth); 

end: 


(1)  UINDSPEED 

(2)  SOUND  SPEED  PROOF  I LE 

(3)  BOTTOM  TYPE 

(4)  DEPTH 

Please  input  depth  (meters): 


',wind:4,'  Knots' 
' ,Prof i leinwords 
'.Bottominwords 

', depth :4,'  m' 


>; 
>; 
); 
>; 

>; 


procedure  GetEnvi  ronmental Data ; 

var 

hisanswer   :  char; 

begin 

while  hisanswer  <>  '0'  do  begin 

clrscr; 

writelnC 

writeln<' 

writelnC 

writelnC 

writeln; 

writeln; 

writeln; 

writelnC 

writelnC 

writelnC 

writelnC 

writelnC 

writeln; 

writelnC 

write  C 

hisanswer 


This  is  the  current  environmental  data: 


(0)  BACK  TO  MAIN 

(1)  UINDSPEED 

(2)  SOUND  SPEED  PROOFILE 

(3)  BOTTOM  TYPE 

(4)  DEPTH 


Please  input  the  numer  of  the 
parameters  you  want  to  change 
readkey; 
case  hisanswer  of 


'1' 
'2' 
'3' 

'4' 

end; 
end; 
end: 


Getwindspeed; 
GetSoundProf i le; 
GetBottomtype; 
GetDepth; 


')* 

',wind:4,'  Knots' 

', Prof i leinwords 

', Bottom! nwords 

', depth :4,'  m' 

'); 

'); 
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begin 


aARGET  INIT> 

sigmarinit 

portmostcourse 

courserange 

minsubspeed 

maxsubspeed 

react range 

react iontack 

motioninodel 


:=DEFSIGMARINIT; 

:=DEFPORTMOSTCOURSE; 

:=DEFCOURSERANGE; 

;=DEFMINSUBSPEED; 

:=DEFMAXSUBSPEED; 

;=DEFREACTRANAGE; 

;=DEFREACTIONTACK; 

:=OMNI; 


{SEARCHER  INIT> 
noisesigina:=DEFNOISESIG; 

{SIMULATION  INIT> 

totalsimtime  :=DEFTOTALSIMTIME; 
WannaContinue:=TRUE; 
scoutsdelay  :=DEFSCCXJTSDELAY; 


{ENVIRONMENT  INIT> 


bottomtype 

svtprof 1 letype 

wind 

depth 

end. 


=DEFBOTTOMTYP; 
=DEFSVTPROFIL; 
=DEFUIND; 
=DEFDEPTH; 


111 


112 


APPENDIX  C.  USER  INTERFACE  SEQUENCES 

1.0   GENERAL  MENU 


Please,  select  one  of  the  following  options: 


0.  QUIT  PROGRAM 

1.  UPDATE  TARGET  DATA 

2.  UPDATE  SEARCHER  DATA 

3.  UPDATE  ENVIRONMENTAL  DATA' 

4.  SET  TIME  VARIABLES 

5.  START  SIMULATION  WITH  CURRENT  DATA 

6.  DISPLAY  INITIAL  TARGET  POSITIONS' 

Enter  your  option: 


2.1    TARGET  PARAMETERS  MAIN  MENU 


Default  values  of  the  parameters  used  to  generate  and  update 
500  simulated  tracks: 


(0)  Back  to  Main 

(1)  DATUM  ERROR  IN  NH  : 

(2)  TARGET  MOTION  MODEL  : 
(speed  &  course  uniformly  distributed) 

(3)  COURSE  LIMITS  (FOR  FAN  MOTION  MODEL  ONLY) 

PORT  LIMIT  (0-360)  : 

RANGE     (Not  starbord  limit)  : 

(4)  SPEED 

MINIMUM  ESTIMATED  TARGET  SPEED  IN  KTS  : 

MAXIMUM  ESTIMATED  TARGET  SPEED  IN  KTS  : 

(5)  TARGET  REACTION  MODEL 

REACTION  RANGE  : 

REACTION  TACK  : 

Enter  the  number  of  the  parameters  you  want  to  change  or  0 

to  go  back  to  main.  Return  is  not  necessary  :  ' 
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2.2    ERROR  INPUT 


Please  input  the  following  data  for  500  targets: 


DATUM  OHNIDIRECTIONAL  ERROR  (integer) 
(Standard  Deviation  in  Tgts  Position,  in  nm 

Current  value  is 
Enter  your  value 


2.3. 

MOTION  MODEL  INPUT 

r-. =rr~                    lARGET  DATA       = : — 

— —  . 

Please  input  the  following  data  for  500  targets: 

TARGET  MOTION  MODEL 

(1)  OMNI  :  Omnidirectional  fleeing 

(2)  FAN  :  Directional  movement 

Current  motion  model  is                  : 

Enter  your  model  (1  OMNI,  2  FAN)            : 

2.4.   TARGET  COURSE  LIMITS  INPUT 


.-...  - --    -  .                                 II 

r—r^n                                                          1 

===:==s::     =====55=========      TARGET  DATA     ==s^==m~~~p                                                      || 

Please  input  the  following  data  for  500  targets: 

TARGET  COURSE  LIMITS 

Please  input  the  courses  port  limit 

The  current  value  is                     : 

Enter  your  value                         : 

1     Please  input  the  course  range 

The  current  value  is                     : 

Enter  your  value                         : 
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2.5.   TARGET  SPEED  ESTIMATE  INPUTS 


-=,  =.-.._.«...^-^.^  ._.=^=-_            1 

!r!r!r!r;::n:::::::.i::::i=--H=^            TARGET  DATA     ====tii====# 

T-T!i:!;!»;:i:!!!;!;ni 

Please  input  the  following  data  for  500  targets: 

TARGET  SPEED  LIMITS 

Please  input  the  Minimum  estimated  Target  speeed 

The  current  value  is                      : 

Enter  your  value                       : 

Please  input  the  Maximum  estimated  Target  speed 

The  current  value  is                     : 

Enter  your  value                         : 

2.6.   TARGET  REACTION  PARAMETERS  INPUTS 


_    _  ^^.=.„„.._._  „  ^_            II 

== — -     -■=— 

s ~ 

TARGET 

DATA 

_^;^fcg^   ■ 

Please  input  the 

following 

data  for  500  targets:                           | 

TARGET  REACTION  PARAMETERS 

Please  input  the 

estimated 

REACTION 

RANGE  in  nautical  miles,                   1 

1     The  current  value  is 

= 

Enter  your  value 

Please  input  the 

estimated  REACTION 

TACK  in  nautical  degrees. 

(Evasion  course  relative  tc 

1  the  bearing  to 

the  closest  scout 

The  current  value  is 

■ 

Enter  your  value 

■ 
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3.1 


SCOUTS  MAIN  MENU 


._.._„ II 

m~M       SCOUTS   DATA       ==§i==^===H5iH!H^ 

■~~^=~iS~~~~-~~~^=~ 

Please  select  one 

of  the  following  options: 

0. 

BACK  TO  MAIN 

1. 

DISPLAY  INFO  RECEIVED  FROM  SCOUTS 

2. 

DISPLAY  AND  OR  CHANGE  DEFAULT  SONAR  PARAMETERS 

3. 

DISPLAY  AND  OR  CHANGE  DEFAULT  ERROR  VALUES 

Enter  your  option: 

3.2.    SCOUTS  DATA  LINK  DISPLAY 


SCOUTS   DATA 


Connected  to  Data  Link,  #  of  scouts  present 

Time  searching  so  far 
Receiving  the  following  initial  conditions 


SCOUT  # 

XPOS   YPOS 

COURSE 

SPEED   SONAR_ON 

SCALE 

Hit  any  key  to  continue. 
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3.3. 


SONAR  PARAMETERS  DISPLAY  AND  INPUT 


Scouts  Sonar  Parameters: 


1        (1)     (2)  1  (3)   I  (4)     (5) 

jsCOUT    SL      DT   I  DI    j   Frq.    PulseL. 
1  #     (dB)     (dB)  j  (dB)  1  (Khz)   (msec) 

(6) 

Beam  U. 
(deg) 

1                       i        i 

You  are  going  to  produce  modifications  to  scout  nimiser: 
(eventual) 


Please 
Please 
Please 
Please 
Please 


nput  new  SCXJRCE  LEVEL 

nput  new  DETECTION  THRESHOLD 

nput  new  DIRECTIVITY  INDEX 

nput  new  FREQUENCY (updates  aborption  coeff) 

nput  new  PULSE  LENGTH(in  ms) 


Please  input  new  BEAH  UIDTH(in  degrees) 


3.4. 


SONAR  EQUATION  ERROR  INPUT 


i==H   SCOUTS   DATA 


ERROR  VARIABLE  DEFAULT  VALUE 


Total  estimated  error  in  the  calculation  of  the  terms  in  the 
calculation  of  the  sonar  equation  and  in  the  uncertainties  in 
the  detection  process  in  dB  re  /yPascal, 


"NOISESIGMA"  current  value  is 
Please  input  the  new  valaue 
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4.1 


ENVIRONMENTAL  MAIN  MENU 


ENVIRONMENTAL   DATA 


This  is  the  current  environmental  data: 


(0)  BACK  TO  MAIN 

(1)  WINDSPEED 

(2)  SOUND  SPEED  PROOF  I LE 

(3)  BOTTOM  TYPE 


4.2.    WIND  SPEED  INPUT 


ENVIRONMENTAL   DATA 


This  is  the  current  environmental  data: 


(1)  WINDSPEED 

(2)  SOUND  SPEED  PROOF  I LE 

(3)  BOTTOM  TYPE 

(4)  DEPTH 


Please  input  wind  speed  (knots): 
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4.3. 

SOUND  SPEED  PROFILE  TYPE  INPUT 

j:r:™__;;~-:::Hn:j:H_r:H::^5:H:-: ___j:H:n:n---:b ==nH:=;:=H^-~SHl;i;^li=:H:n— =HSH--=:::-:r--H;:™S 

j=========nj5=g#==s=i=:=      ENVIRONMENTAL   DATA      :==:aPHSi====|:is.:^ 

This  is  the  current  environmental  data: 

(1)  UINDSPEED           : 

(2)  SOUND  SPEED  PROOF I LE 

(3)  BOTTOM  TYPE 

(4)  DEPTH              : 

You  can  choose  the  sound  speed  profile  according  to  the  following 

code: 

(1)  POSITIVE  OR  ISOSPEED 

(2)  NEGATIVE 

Please  enter  your  option  here  : 

4.4.         BOTTOM  TYPE  INPUT 


— ,.  — ,.,■  _.   II 

======:=n=======       bNVIKONMbl 

.„_ .. ._.... ..___:_™.: u.._..i. 

HnU:i: 

<TAL     DATA          _:=:=:-_::=--.,=- 

:sn::K 

This 

is  the  current  bottom  &  depth  environmental  data: 

(1)  WINDSPEED 

1 

(2)  SOUND  SPEED  PROOF  I LE 

1 

(3)  BOTTOM  TYPE 

(4)  DEPTH              ! 

You  can  input  bottom  type  according  to  the  following  code 

(1)  SANDY, FLAT  BOTTOM, NEGLECT ABLE  ROUGHNEESS 

1 

(2)  INTERMEDIATE  BETWEEN  FLAT  AND  VERY  ROUGH 

(3)  ROCKY,  VERY  ROUGH 

PI 

ease  enter  your  option  here  :   ' 
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4.5.    DEPTH  INPUT 


=~="-TTnhrT:r-t;-;r=-rT" — r-r^--=™="=-="===; =_.::.•::=: 

:n=nn=:=::r:n:::nr:::=n=:::nn:::nn:::s=Ks::::::::KSS;:n:::s=::s:K: 

-■"""■"'"--"•''-^ — ■■  ■ :"""""":";:! 

=^=,,,^:H:™:^:i:::4::H      ENVIRONMENTAL   DATA      ====-=======.===:===== 

This  is  the  current  environmental  data: 

(1)  WINDSPEED 

(2)  SOUND  SPEED  PROOF I LE 

(3)  BOTTOM  TYPE 

(4)  DEPTH 

Please  input  depth  (meters):   ' 

5.0    TIME  VARIABLES  INPUT 


:_.=^.,.,  --,,—==      SIMULATION    TIME 

_..=: r — 

1     The  simulation  will  be  carried  out  from  time  0  to  TotalSimtime                || 

in  updating  steps  of  5  minutes.  Totalsimtime  should 

be  smaller                 || 

than  the  last  updating  time  to  the  data- link  device 

file. 

Total  Simulation  Time  (TotalSimtime,  a  multiple  of 

5  min) 

No  info  from  scouts  available  after         : 

Current  value  is                        : 

Enter  your  value                       : 

Hit  any  key  to  continue... 

Delay  Time  of  scouts  after  DATUM  time, (DATUM  "age" 

in  min) 

Current  value  is                        : 

Enter  your  value                       : 
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APPENDIX  D.   TARGETMANAGER  UNIT  CODE 


unit  TargtManager; 

i  nterf ace 

uses 

GRAPH,  CRT; 
var 

function  normaUvar  u:double):real; 

procedure  GenTargets; 

procedure  TargetReaction(ti,ss:integer;distst:real); 

procedure  UpdateTargetPosition(ti : integer); 

procedure  UpdateUeightsUithBayes; 

implementation 

uses  GraphManager.InterManager; 

function  normal (var  u: double): real; 

(converts  uniform  random  nutrbers  to  standard  norma l> 

{by  Alan  Uashburn> 

const 

a=2.515517;b=0.802853;c=0.010328;d=1.432788;e=0.189269;f=0. 001308; 
var 

sign: integer; 

x,y:real; 
begin 

sign:=-1; 

if  u>0.5  then  begin 
sign:=1; 
u:=1.0-u; 

end; 

y:=sqrt(-ln(u*u)); 

x:=(a+y*(btc*y))/(1+y*(d+y*(e+f*y))); 
normal :=sign*(y-x); 
end; {normal} 


function  bearing(xl,yl.x2.y2: integer): integer; 

CFROM  1  TO  2> 

var 

x,y   : integer; 

begin 

x:=x2-x1; 

y:=y2-y1; 

if  x>0  then  begin 

if  y  >  0  then  bearing  :=  round(180  -  arctan(x/y)*180/pi); 

if  y  =  0  then  bearing  :=  90; 

if  y  <  0  then  bearing  :=  -round(arctan(x/y)*180/pi); 
end; 
if  x<0  then  begin 

if  y  <  0  then  bearing  :=  round(360  -  arctan(x/y)*180/pi ); 

if  y  =  0  then  bearing  :=  270; 
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if 

y  >  0  then  beari 

ing  := 

roun 

end; 

if  x=0  then  begin 

if 

y  >=0  then  beari 

ng  := 

180; 

if 

y  <  0  then  beari 

ng  := 

0; 

end; 

end; 

procedure  GenTargets: 

var 

r 

:double; 

i 

:word; 

datunerror  : integer; 

range 

:real; 

round(180  -  arctan<x/y)*180/pi); 


begin 

Randomi  ze; 

for  i:=1   to  NUMBEROFTGTS  do  begin 

r  "random; 

mytargetCi] .X  :=  CENTER  +  round(norinal(r)  *  sigmarinit); 

r:=randoin; 

mytarget[i] .y  :=  CENTER  +  round( normal (r)  *  sigmarinit); 

mytargetEil .datimxrsmytargetEi] .x; 

roytarget[i] .datumy:=mytarget[i] .y; 

mytarget [i] .course:=portmostcourse+random(courserange); 

mytarget [i] .speed:=minsubspeedi-random*(maxsubspeed-minsubspeed); 

mytarget [i] .x:=mytarget[i] .X  +  round(mytarget[i] .speed  *  scoutsdelay  * 

sin(mytarget[i] •course*pi/180)); 

mytarget [i] .y:=mytarget[i] .y  -  round(mytarget[i]  .speed  *  scoutsdelay  * 

cos(mytarget[i] .course*pi/180)); 

mytarget [i] .eel Ix  :=(mytarget[i] .x  div  GRIDSTEP)+1; 

mytarget [i] .eel ly  :=(mytarget[il .y  div  GRIDSTEP)+1; 

mytarget [i] .prob:=0.002; 

mytarget [i] .threatdis:=BIGNUMBER; 
end; 
end; 


procedure  UpdateTargetPosition(ti  rinteger); 

var 

r  : double; 

begin 

mytarget [ti]  .x:=mytarget[ti] .X  +  round(mytarget[ti]  .speed  * 

sin(mytarget[ti] .course*pi/180))  ; 

mytarget [ti] .y:=mytarget[ti] .y  -  round(mytarget[ti]  .speed  * 

cos(mytarget[ti] .course*pi/180)); 

mytarget [ti] .eel Ix  :=(mytarget[ti] .x  div  GRIDSTEP)+1; 

mytarget [ti] .eel  ly  :=(mytarget[ti] .y  div  GRIDSTEP)+1; 
end: 


procedure  TargetReaction(ti .ss:integer:distst:real); 

var 

filehandle  :text; 

r  :double; 

bearts  : integer; 

begin 

if  distst  <  reaetrange  then  begin 

bearts:=bearing(mytarget[ti] . x, mytarget [ti] .y,myscout[ss] .x,myseout [ss] .y); 

if  bearts  >  360  then  bearts:=bearts-360; 

if  bearts  <  0       then  bearts:=360-bearts; 
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if  sin((niyscout[ss] .course-bearts)*pi/180)  >  0  then  begin 

mytargetrti] .course:=bearts  -   react iontack; 
end  else  begin 

mytargetCti) .course:=bearts  +  reactiontack; 
end; 
end  else  begin 

if  motionmodel  =  'OMNI'   then  begin 

mytargetCti] .course:=bearing(mytarget[ti] .datumx.mytarget [ti] .datumy, 

mytargetCti] .x,mytargettti] .y); 
end  else  begin 

mytarget [ti] .course:=portmostcourse«-random(courserange); 
end; 
end; 
if  mytarget [ti] .course>360  then  begin 

mytarget [ti] .course:=mytarget [ti] .course-360; 
end; 
if  mytarget [ti] .course<0       then  begin 

mytargetCti] . course: =360+mytargetCti] .course; 
end; 

if  mytargetCti] .threatdis  >  distst  then  begin 
mytargetCti] .threatdis:=distst; 
mytarget [ti] .evasioncourse:=mytarget[ti] .course; 
end; 

mytarget [ti] .course:=mytarget[ti] .evasioncourse; 
if  ss=numberofscouts  then  mytargetCti] .threatdis:=BIGNUMBER; 
end ; {TargetReacti  on} 


procedure  UpdateWei  ghtsWi  thBayes ; 

var 

i  : integer; 

begin 

for  i:=1  to  NUMBEROFTGTS  do  mytarget Ci] .prob:=mytarget[i]  .prob/totprob; 

end: 


begin 
end. 
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APPENDIX  E.   SCOUTMANAGER  UNIT 

unit  ScoutManager : 
i  nterf ace 
uses 

GRAPH,  CRT; 

procedure  Im'tScouts; 
procedure  CloseScoutsDataLink; 
procedure  ConnectWithScoutsThriJ)ataLink; 
procedure  UpdatePriorUithScoutsInfo; 
procedure  ReadScoutsDataLink; 
procedure  ResetUndetectedMass; 
procedure  UpdateUndetectecMass; 

implementation 

uses  GraphManagerJargtManager. Intermanager.DetecManager; 

procedure  Ini ti  a1 i  zeScoutsVarsWi  thDef aul ts ; 

var 

k       : integer  ; 

begin 

for  k:=1  to  nunberof scouts  do  begin 

myscout [k] .abscoef f  :=DEFabscoef f ; 

myscout [k] .SL       :=DEFSL; 

myscout [k].NL  :=DEFNL; 

myscout [k] .DT  :=DEFDT; 

myscout Ckl.DI       :=DEFDI; 

myscout [k].SS       :=DEFSCATTSTRG; 

myscout  tk] . IS_HOT    :=FALSE; 

myscout [k] .frequency  : =DEF FREQUENCY; 

myscout [k] .beamwidth  :=DEFBEAMUIDTH; 

myscout [k] .pulselength:=DEFPULSELENG; 

myscout [k] .reverbarea:=DEFREVRBAREA; 
end; 
end; 


procedure  Determi  neNumberOfScoutsAndTi  me ; 

var 

1  : integer; 

begin 

readln(datal ink, numberof scouts, dlinktotaltime); 

for  i:=1  to  numberof scouts  do  begin 

readlnCdatal ink, myscout [i] .id); 
end; 
end; 


procedure  InitScouts: 

begin 

assign(datal ink, 'c:\lnk1.dat'); 

reset(datalink); 

Determi  neNumberOf ScoutsAndT  i  me; 

close(datalink); 

I ni t  i  a  I i  zeScoutsVarsWi  thDef aul ts; 

end; 
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procedure  CloseScoutsDataLink; 

begin 

close(datalink); 

end; 


procedure  ConnectWi  thScoutsThruDataLi  nk : 

begin 

reset(datalink); 

DetermineNumiaerOfScoutsAtxJTiine; 

ReadScoutsDatal ink; 

end; 


procedure  ReadScoutsDataLink; 

var 

sonar info  :byte; 

contact  i  nf o: byte; 

k        :byte; 

begin 

for  k:=1  to  numberofscouts  do  begin 
read(datalink,myscout [k] .x); 
read(datalink); 
read(datalink,myscout[k] .y); 
read( data  I i  nk , myscout [k] . course) ; 
read(datalink,niyscout  [k]  .speed); 
read(datalink,sonarinfo); 
read(datal ink, myscout [k].maxscale); 
read(data I i  nk, contact  i  nf o) ; 
if  sonarinfo  =  1  then  begin 

myscout [k] .SONAR_ON:=TRUE; 
end; 
if  contact  info  =  1  then  begin 

if  not  myscout [k] . I S_HOT  then  begin 

myscout [k] .initialhot:=time*5; 
end; 

myscout [k] . I S_HOT  : =TRUE ; 
sound<1500);delay(40); 
sound(1600);delay(10); 
sound(800);delay(200); 
sound(1600);delay(60); 
nosound; 

readC data link, myscout [k] .subx); 
readCdatal ink, myscout [k] .suby); 
read(datal ink, myscout [k] .subcourse); 
readCdatal ink, myscout [k] .subspeed); 
readCdatal ink, myscout [k] .cred); 
end  else  begin 

myscout [k] . I S_HOT      : =FALSE ; 
end; 

readlnCdatalink); 
end; 
end: 


procedure  TempUpdateNegativeInfo(s:byte) ; 

var 

i       : integer; 

cdp5     : rea I ; 

d       ireal; 

begin 

for  i:=1  to  NUMBEROFTGTS  do  begin 
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d:=sqrt(distancer(mytarget[i]  .x,niytargetCi]  .y,myscout[s]  .x,myscout[s]  .y)); 
if  d  <  myscout[s] .maxscale     then  begin 

cdp5 : =Detect  i  onEf f ort  <  i , s , d, depth ) ; 

mytarget  [i]  .prob:=iiiytarget  [i]  .prob  *  {1-cdp5); 
end; 
if  s=numberof scouts  then  begin 

UpdateTargetPosition(i); 

totprob:=totprob  +  mytarget [i] .prob; 
end; 
end; 
end: 

procedure  TenipUpdatePositiveInfo(s:byte) : 

var 

1  : integer; 

cdp5  : rea I ; 

C  : rea I ; 

d  :real; 

radweight:real; 

distoff     :real; 

error         :real; 

begin 

for  i:=1   to  NUMBEROFTGTS  do  begin 

d:=sqrt(distancer(mytarget[i] . x, mytarget [i] .y,myscout[s] .x,myscout[s] .y)); 

cdp5 :=DetectionEff ort (i,s,d, depth); 

C:=myscout[s] .cred; 

distoff :=sqrt(distancer(mytarget[i] . x, mytarget [i] .y,myscout[s] .subx,myscout[s] .suby)); 

error:=  d  *  myscout[s] .beamwidth  *  pi  /180; 

radweight:=exp(-sqr(distoff/error)/2); 

mytarget ti] .prob:=  mytarget [i] .prob  *  (1-C) 

+mytarget[i] .prob  *  cdp5  *  radweight*C; 

TargetReaction(i,s,d); 

if  s=numberof scouts  then  begin 
UpdateTargetPosi t i on( i ) ; 
totprob:=totprob  +  mytarget [i] .prob; 

end; 
end; 
end: 

procedure  UpdatePriorWithScoutsInfo; 

var 

scout      : integer; 

begin 

totprob:=0. 0000000000001 ; 

for  scout :=1  to  numberof scouts  do  begin 
if  odd(time)  then  begin 

UpdateNL( scout); 
end; 
if  not  myscout [scout] . I S_HOT  then  begin 

TempUpdateNegat  i  vel nfo( scout ) ; 
end  else  begin 

TempUpdatePos i t i vel nfo( scout ) ; 
end; 
end; 
end: 

procedure  ResetUndetectedMass: 
begin 
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undetectedmass:=1 ; 
end; 


procedure  UpdateUndetectedMass: 
begin 

undetectedtnass:=undetectedmass*totprob; 
end: 


begin 
end. 
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APPENDIX  F.  DETECMANAGER  UNIT 

unit  DetecManager; 

interface 

uses 

CRT, Intermanager.scoutmanager; 

function  Detect ionEffort( id, sd  : integer;distts, dp: real): real, ■ 

function  AbsorptionCoefficient(f: real): real; 

function  UpdateNL(sc:byte): integer; 

function  ScatteringArea(tau,theta: real): integer; 

function  ScatteringStrength(svtprof i le: integer; freq: real; winspeed, bottom: integer): integer; 

f unct  i  on  di  stancerCxI ,y1 , x2, y2 : rea I ) : rea I ; 

procedure  UpdateSonarParameters; 

implementation 

{  TOOLS} 

{ } 

function  distancer(xl.yl.x2.y2:real):real ; 

var  d1,d2:real; 

begin 

d1:=(x1-x2)*(x1-x2); 

d2:={y1-y2)*(y1-y2); 

di  stancer : =d1+d2; 

end: 

function  power2(  base    :  real: 

exponent  :  real  )  :  real : 
{base  can  not  be  negative> 
begin 

power2:=  exp(  exponent  *  ln(  base  )  ); 
end:  {  function  power  } 

function  nornialcdfwashb(mu,sig,x:real):real: 

Cbased  in  a  program  by  Alan  Uashburn> 

const 

A=0. 319381530;   B=-0, 356563782;   C=1. 781477937;  D=-1 .821255978; 

E=1. 330274429;  G=  0.231641900; 

var 

y,f,p         :real; 

positive  : boo lean; 

begin 

if  sig<=0  then  sig:=1; 

x:=(x-mu)/sig; 

positive:=(x>=0); 

x:=abs(x); 

if  x>=100  then  x:=100; 

y:=1/(1+G*x); 

f:=C+y*(D+y*E); 

p:=1-exp(-0.5*x*x)/sqrt(8*arctan(1))*y*(A+y*(B+y*f)); 

if  positive  then  begin 
nornialcdfwashb:=p; 

end  else  begin 

normalcdfwashb:=1-p; 

end; 
end: 
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function  sumdB(a.b:real):real; 
begin 

sumdB:=(10/ln(10))*ln(   power2<10,a/10>  +   power2(10,b/10)  ); 
end: 


{SCAHERING  METHODS} 

{ 

function  TargetStrength(tgx,tgy,scx.scy,tgc:real) :integer; 

const 

LOWTS=10; 

INCTS=8; 

var 

cospsisqr       :real; 

num,den        :real; 

tstemp         : integer; 

begin 

if  (sex  <>  tgx)  or  (scyo  tgy)  then  begin 

num:=sqr((scx-tgx)*sin(tgc*pi/180)+(scy-tgy)*cos(tgc*pi/180)); 

den:=({scx-tgx>*(scx-tgx)+(scy-tgy)*(scy-tgy>); 

cosps  i  sqr :  =nijn/den; 

tsteiiip:=round(L0WTS+INCTS*(1 -cospsisqr)); 
TargetSt  rength : =tstenip; 
end  else  begin 
TargetStrength:=LOWTS; 
end; 
end: 


function    ScatteringArea(tau. theta: real ): integer; 

const 

C=1500; 

begin 

theta:=theta*pi/180; 

ScatteringArea:=round(  (10/ln(10)  )  *  ln(tau*theta*C/2)  ); 
end; 


function  SurfaceScattering(fr:rea1 ; wind: integer): integer; 

const 

THETA=37*pi/180; 

var 

h       ireal; 

begin 

h:=power2(wind,2.5)*0.026; 

SurfaceScattering:=  round<(10/ln(10))  *  power2(ln(f r*h*sin(THETA)),  0.99)  -  45.3); 
end; 


function  BottomScattering(fr:rea1 ;bot:integer) : integer; 
begin 

case  bot  of 

1  :   BottoinScattering:=round(  0.0120*fr*f r*fr  +  0.0283*fr*fr  +0.0062*fr  -  28.153); 

2  :  BottomScattering:=round(  0.0076*fr*f r*fr  -  0.2337*fr*fr  +2.8700*fr  -  40.410); 
2       :  BottoniScattering:=-18; 

end; 
end; 
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f uncti on  Scatteri ngStrength( svtprof i 1 e : i nteger ;  f req : real : wi  nspeed . bottom : i nteger) : i  nteger ; 
<svtprofile=1  (positive  velocity  gradient)> 
<svtprof ile=2  (negative  velocity  gradient)> 
begin 

case  svtprof ile  of 

1  :  ScatteringStrength:=SurfaceScattering(freq,winspeed); 

2  :  ScatteringStrength:=Botto(nScattering(freq, bottom); 
end; 

end; 

procedure  UpdateSonarParatneters; 

var 

i: integer; 

begin 

for  i:=1  to  nimberofscouts  do  begin 
myscout[i] .reverbarea:=ScatteringArea(niyscout[il .pulselength, 

myscout [i] .beamwidth); 
myscoutCi] .SS:=ScatteringStrength(svtprofiletype, 

myscout [i] . frequency, wind, bottomtype); 
end; 
end: 


{NOISE  METHODS} 
{ 


function  AmbientNoiseCf :real :w: integer): real : 
begin 

AmbientNoise:=46. 121+2. 221*w-0.041*w*w-17.012*(  ln(f)/ln(10)  -  3  ); 
end: 


function  Sel fNoi se(v: real :f:real): real : 
begin 

SelfNoise:=  23+2.02*v+(20/ln(10))*ln(25/f ); 
end; 

function  NoiseLeveKscoutvel :real :f req: real :windspeed: integer): integer: 

var 

An  : rea I ; 

Sn         :real; 

begin 

An:=Anibi  entNoi  se(  f  req,  ui  ndspeed) ; 

Sn:=SelfNoise(scoutvel,freq); 

NoiseLevel:=round(sumdB(An,Sn)); 
end: 


function  UpdateNKsc: byte) : integer: 
begin 

with  myscout [sc]  do  begin 

NL:=  NoiseLevel (speed, frequency, wind); 

end; 
end; 
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{PROPAGATION  METHODS} 

{ 

function    AbsorptionCoefficient(f :real):real : 
begin 

AbsorptionCoefficient:=        8e-3  /(  (0.7  /f*f)+1   ) 

+     0.04/((6000/f*f)+1)        +    4e-7  *  f*f  ; 
end; 


function  TransLossesGeom(d,H:real):real : 

var 

TLG     rreal; 

begin 

if  d  <=  0  then  begin 
TLG:=1; 

end  else  begin 

TLG:=(10/ln(10>)*(ln<d)+ln(H/pi)); 
end; 

TransLossesGeoni:=TLG; 
end; 


{DETECTION  MODEL  METHODS} 

{ 

function  DetectionEffort(id.sd:integer:distts. dp: real) rreal ; 

var 

p  :realj 

TL  :realj 

TLGEOM  treat; 

FOM  :real; 

SE  :real; 

TS  :real; 

SErevbr  rreal; 

SEnoise  :real; 

tx  : integer; 

ty  : integer; 

sx  : integer; 

sy  : integer; 

tc  :real; 


begin 

tx:=mytarget[id] .x; 

ty:=mytarget[id] .y; 

sx:=niyscout  [sd]  .x; 

sy:=myscout [sd] .y; 

tc:=mytarget[id] .course; 

distts:=185*distts;<:i85  accounts  for  the  scale  factor> 

TLGeom:=  TransLossesGeom(distts,dp); 

TL:=TLGeom+distts*myscout[sd] .abscoeff; 

TS :=  TargetSt rength( tx, ty, sx, sy, tc) ; 

FOM:=  myscout [sd] .SL-myscout [sd] .NL-myscout [sd] .DT+TS+myscout [sd] .DI ; 

SEnoise:=F0M-2*TL; 

SErevbr :=TS-TLGeom-myscout[sd] . DT -myscout [sd] . reverbarea-myscout [sd] .SS; 

if  SErevbr  <  SEnoise  then  begin 

SE:=  SErevbr; 

end  else  begin 

SE:=  SEnoise; 
end; 
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p:=power2(1-nornialcclfwashb(0,   NOISESIGMA,SE),   NUMOFPULSES); 
Detect ionEf fort: =1-p; 
end; 


begin 
end. 
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APPENDIX  G.  CELLSMANAGER  UNIT 


unit  CellsManager: 

i  nterf ace 

uses 

GRAPH,  CRT  ; 


type 

probcelltype  =  array [1.. 80,1.. 80]  of  single; 


var 

probincell  : probcelltype; 

maxprob  :real; 

numberout  : integer; 

procedure  UpdatelnCellProbabilities; 
procedure  UpdateOutAreaProbabi I i ties; 
procedure  WriteProbinCell; 


implementation 

uses  GraphManager.  Targtmanager.InterManager; 


procedure  Clearance! Is; 

var 

i,j       : integer; 

begin 

for  i:=1  to  MAXGRID  do 

begin 

for  j:=1  to  MAXGRID  do 
begin 
probincel I [i , j] :=0; 

end; 
end; 
end: 


procedure  Updatel nCel 1 Probabi 1 i  ti  es : 
var 

i   : integer; 
begin 

maxprob: =0; 
numberout :=0; 
ClearallCells; 
for  i:=1  to  NUMBEROFTGTS  do 
begin 
if((mytarget[i].x  >  GRIDSTEP)  and  (mytarget [i] .y  >  GRIDSTEP  ))  and 

((mytarget[i].x  <MAXGRID*GRIDSTEP)  and  (mytarget [i] .y  <MAXGRID*GRIDSTEP  )) 

then 

begin 

probincell [mytarget [i] .eel Ix, mytarget [i] .celly] := 

probincell [mytarget [i] .eel Ix, mytarget [i] .cellyl+mytarget [i] .prob; 

if  maxprob  <  probincel I [mytarget [i] .eel Ix, mytarget [i] .celly]  then 
begin 


135 


maxprob:=  probincell [mytarget[i] .cellx,mytarget[i] .celly]; 
end; 
end  else 
begin 

inc(nijnberout); 
end; 
end; 
end; 


procedure  BroadcastProbinCell: 

var 

i,j       : integer; 

begin 

assign(outdatal ink, "c:\lnk2.dat"); 

for  i:=1  to  MAXGRID  do 

for  j:=1  to  MAXGRID  do 

writeln(outdatalink, probincell [i, j]  }; 

end; 
end; 
end; 


begin 
end. 
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APPENDIX   H.GRAPHMANAGER   UNIT 


unit  GraphManager: 

interface 

uses 
GRAPH,   CRT; 

const 

BOXSIZE=100; 

SCALECOLORS=10; 

BGIdir='..\BGI'; 

GRIDCOLOR=32; 

BASEPATTERN=1; 

TEXTC0L0R=8; 

TARGETSCOLOR=30; 

SCOUTSCOLOR=24; 

PALETTE1=11; 

PALETTE2=12; 

MAXGRID=80; 

GRIDSTEP=5; 

CENTER=MAXGRID*GRIDSTEP  div  2+GRIDSTEP; 

STEP        =  5; 

FRAME        =80; 

type 

colortype      =  array [1..SCALEC0L0RS] of  byte; 

celltype       =  array [1..MAXGRID,1..MAXGRID]  of  integer; 

const 

COLORCOOE    :  colortype  =  (8,1,33,12,4,36,52,62,55,63); 


procedure  InitGraphics; 

procedure  SetGrid; 

procedure  paint(x,y,pat,col:word); 

procedure  SetScale; 

procedure  SetTimed: integer); 

procedure  DisplayTargets; 

procedure  DisplayScouts; 

procedure  PleaseWait; 

procedure  SetProbDistrColors; 

procedure  OutputGraph; 

procedure  Start; 

procedure  DisplayData; 

procedure  DisplayContactSuninary(srchr:byte); 


implementation 

uses  TargtManager . ScoutManager , Cel 1 sManager , I ntermanager : 


procedure  InitGraphics; 

var 

GraphD river, 

GraphMode, 

i        :  integer; 
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begin 

GraphD river  :=  GRAPH. Detect; 

GRAPH.InitGraph(GraphDriver,GraphMode,BGIdir); 

for  i:=1  to  SCALECOLORS  do  begin 

SetPalette(i,COLORCODE[i]); 

end; 
encl;{initgraphics} 


procedure  SetGrid: 

var 

i       : integer; 

oldstyle  :TextSettingsType; 

begin 

GetTextSettings(oldstyle); 

SetPa I et te( 1 5 , GR I DCOLOR ) ; 

for  i:=1  to  NAXGRID  do  begin 

line(i*GRIDSTEP,GRIDSTEP,i*GRIDSTEP,GRIDSTEP*MAXGRID>; 

line(GRIDSTEP,i*GRIDSTEP,GRIDSTEP*MAXGRID,i*GRIDSTEP); 
end; 

SetTextStyle(2,0,(STEP+1>  div  2  ); 

line<GRIDSTEP+MAXGRID*GRIDSTEP  div  2,0,GRIDSTEP+MAXGRID*GRIDSTEP  div  2,GRIDSTEP); 
OutTextXY(MAXGRID*GRIDSTEP  div  2,0,chr(48)); 
OutTextXY(0,MAXGRID*GRIDSTEP  div  2,chr(48)); 
OutTextXY<MAXGRID*GRIDSTEP  div  4,0,concat(chr(49),chr{48))); 
OutTextXY(0,MAXGRID*GRIDSTEP  div  4,concat(chr(49),chr(48))); 
OutTextXY(MAXGRID*GRIDSTEP*3  div  4,0,concatCchr(49),chr(48))); 
OutTextXY(0,MAXGRID*GRIDSTEP*3  div  4,concat(chr(49),chr(48)>); 
OutTextXY(MAXGRID*GRIDSTEP-GRIDSTEP,0,concat(chr(50),chr(48))); 
OutTextXY(0,0,concat(chr(50),chr(48))); 

OutTextXY(0,MAXGRID*GRIDSTEP-GRIDSTEP,concat(chr(50),chr(48))); 
SetColordS); 
end: 


procedure  paint(x.y.pat,col :word); 
begin 

SetFillStyle(pat,  col); 
FloodFill(x,y,  GetMaxColor); 
end; 

procedure  SetScale; 

var 

i       : integer; 

outstring:string; 

begin 

line(STEP,FRAME*STEP+STEP*3, 

STEP*2*SCALEC0L0RS+STEP,FRAME*STEP+STEP*3); 
I ine(STEP, FRAME*STEP+STEP*5, 

STEP*2*SCALEC0LORS+STEP , FRAME*STEP+STEP*5 ) ; 
for  i:=1  to  SCALECOLORS+1  do 
begin 

I ine( i*2*STEP-STEP, FRAME*STEP+STEP*3, 
i  *2*STEP - STEP , FRAME*STEP+STEP*5 ) ; 
end; 

for  i:=1  to  SCALECOLORS  do 
begin 

paint( i*2*STEP, FRAME*STEP+STEP*4,BASEPATTERN, i ); 
SetTextStyle(2,0,(STEP+2)  div  2  ); 
OutTextXY( i  *2*STEP- 1 0 , FRAME*STEP+STEP, 
concat(chr(46>,chr(47+i))); 
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end; 
OutTextXY((SCALECOLORS+1)*2*STEP-STEP,FRAME*STEP+STEP,chr(49)); 

end; 


procedure  DisplayTargets; 

var 

i      : integer; 

begin 

SetPalette(PALETTE1,TARGETSC0LOR); 

for  i:=1  to  NUMBEROFTGTS  do 

begin 

if  <(nrytarget[i]  .X  >GRIDSTEP)  and  (mytargetCi] .x 
((mytargetli]  .y  >GRIDSTEP>  and  (mytargetCi] .y 


<MAXGRID*GRIDSTEP})  and 
<MAXGRID*GRIDSTEP))  then 


PutPixel<mytarget[i) .x, mytargetCi] .y.PALETTEl); 


end; 
end; 


procedure  DisplayScouts; 
var 

i      : i  nteger; 
begin 

SetPalette(PALETTE2,SCCXJTSC0L0R); 
for  i:=1  to  nunberof scouts  do 
begin 

if  ((myscoutEil.x  >GRIDSTEP)  and  (inyscout[i) .x  <MAXGRID*GRIDSTEP))  and 
((myscout[i].y  >GRIDSTEP)  and  (myscoutCi] .y  <MAXGRID*GRIDSTEP))  then 
begin 

PutPixeUmyscoutEi]  .x,myscout[i]  .y,PALETTE2); 
SetColor(PALETTE2); 
OutTextXY(nTyscout[i] .X  +  5,  myscoutCi] .y  -  5,  concat(myscoutCi] .id 

,chr(48),chr{48),chr(48+i))); 
Ci re  I e( myscoutCi] .x, myscoutCi] .y,nnyscoutCi] .maxscale); 
if  myscoutCi] .IS_HOT  then  begin 


PutPixeUmyscout 
PutPixeUmyscout 
PutPixeUmyscout 
PutPixeUmyscout 
PutPixeUmyscout 
PutPixeUmyscout 
PutP  i  xe I (myscout 
PutPixeUmyscout 
PutPixeUmyscout 
PutPixeUmyscout 
PutPixeUmyscout 
PutPixeUmyscout 
PutPixeUmyscout 


] , subx+5 ,  myscout  C  i ] . suby+5 , 8 ) 
] .  subx-*-6 ,  myscout  C  i  ] .  suby+4 , 8 ) 
] . subx+7 ,  myscout  C  i ] . suby+3 , 8 ) 
] . subx+8 ,  myscout  C  i ] . suby+2 , 8 ) 
] . subx+9 ,  myscout  C  i ] . suby+ 1,8) 
] .subx+10,  myscout  Ci] .suby+0,8); 
].subx+11,  myscoutCi] -suby-l, 8); 
].subx-i-4,  myscoutCi]  .suby+4, 8) 
] .subx+3,  myscout  Ci] .suby+3, 8) 
] .subx+2,  myscout  Ci] .suby+2, 8) 
] .subx+1 ,  myscout  Ci] .suby+1 ,8) 
].subx+0,  myscoutCi] .suby+0,8) 
].subx-1,  myscoutCi] .suby- 1,8) 


SetColor(14); 

OutTextXY(myscoutCi] .subx  +5,  myscoutCi] .suby- 10,  concat(chr{48+i))); 

SetColordS); 


end; 
end; 


end; 
end: 


procedure  SetProbDistrColors: 

var 

i,j,c   : integer; 

begin 
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if  tnaxprob  >  0  then  begin 

for  i:=1  to  MAXGRID-1  do  begin 
for  j:=1  to  MAXGRID-1  do  begin 

for  c:=0  to  SCALECOLORS-1  do  begin 
if  (100*probinceU [i, j]/maxprob>c*SCALECOLORS) 

and  (100*probincell[i,j]/inaxprob<=(c*SCALEC0L0RS+100/SCALECOLORS))  then  begin 
paint( i*GRIDSTEP+1 , j*GRIDSTEP+1 ,BASEPATTERN,c+1 ); 

end; 
end; 
end; 
end; 
end  else  begin 

OutTextXY(40,40,'No  target  in  the  area, it  is  gone'); 
end; 
end; 


procedure  PleaseWait: 

begin 

InitGraphics; 

outtextxy(MAXGRID*GRIDSTEP+10,  GRIDSTEP, 

end: 


'Please  wait..'); 


procedure  DisplayLogo: 
van 

i  : integer; 
MAXx  : integer; 
MAXy  : integer; 
begin 

SetPalette(15,GRIDC0L0R}; 
MAXx:=GetMaxX; 
MAXy:=GetMaxY; 

I i ne< round( (MAXx/2 ) - BOXS I ZE ) , round< (MAXy/2 ) - BOXS I ZE ) , 
round( (MAXx/2 )+BOXS I ZE ) , round< ( MAXy/2 > - BOXS I ZE ) ) ; 
I i ne( round( (MAXx/2 ) - BOXS I ZE ) , round( ( MAXy/2 )+BOXS I ZE ) , 
round( CMAXx/2 )+BOXS I ZE ) , round( ( MAXy/2 )+BOXS I ZE ) ) ; 
line(round((MAXx/2)-BOXSIZE),round((MAXy/2)-BOXSIZE), 
round( ( MAXx/2 )- BOXS I ZE ), round( ( MAXy/2 )+BOXS I ZE )) ; 
line(round((MAXx/2)+BOXSIZE),round((MAXy/2)-BOXSIZE), 
round(<MAXx/2)+BOXSIZE),round((MAXy/2)+BOXSIZE)); 
line(round((MAXx/2>-BOXSIZE),round((MAXy/2)+BOXSIZE+5), 

round((MAXx/2>+BOXSIZE),round((MAXy/2)+BOXSIZE+5)); 
line(TOund((MAXx/2)-BOXSIZE),round((MAXy/2)+BOXSIZE+25), 
round((MAXx/2)+BOXSIZE),round(<MAXy/2)+BOXSIZE+25)); 
for  i:=0  to  SCALECOLORS  do 
begin 
line(round((MAXx/2)-BOXSIZE)+i*round(BOXSIZE/5),round((MAXy/2)+BOXSIZE+5), 

round((MAXx/2)-BOXSIZE)+i*round(BOXSIZE/5),round((MAXy/2)+BOXSIZE+25)); 
end; 
for  i:=1  to  SCALECOLORS-1  do 
begin 
paint(round((MAXx/2)-BOXSIZE)+1+i*round(BOXSIZE/5),round((MAXy/2)+BOXSIZE+6), 

BASEPATTERN+8,i); 
delay(160); 
end; 
SetPalette(1,8); 

paint(round(MAXx/2),round<MAXy/2),BASEPATTERN,1); 
SettextStyle(1,0,1); 

OuttextXY(round(MAXx/2)-50,round(MAXy/2)-50,'  NRTCAS  '); 
SettextStyle(2,0,A); 
OuttextXY(round(MAXx/2)-100,round(MAXy/2)-25,   '       Reactive  Target'); 
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OuttextXY(round(MAXx/2)-100,round(MAXy/2)-15   ,'  Computer  Assisted  Passive  Search'); 

SettextStyle(2,0,4); 

OuttextXY(round(MAXx/2)-100,round(MAXy/2),     '         Release  1.1  '); 

OuttextXY(round(MAXx/2)-100,round(MAXy/2)+35,      '  NPS  '); 

SettextStyle(2,0,3); 

OuttextXY(round(MAXx/2)-100,round(MAXy/2)+70,     '  C.RECALDE  '); 

end: 


procedure  Di  spl ayl ntro ; 

var 

MAXx     : integer, ■ 

MAXy    : integer; 

begi  n 

cleardevice; 

SetPalette(15,GRIDC0L0R); 

MAXx:=GetMaxX; 

MAXy:=GetMaxY; 

line(GRIDSTEP,GRIDSTEP,MAXx-GRIDSTEP,GRIDSTEP); 

line(GRIDSTEP,GRIDSTEP,GRIDSTEP,MAXy-GRIDSTEP); 

line(HAXx-GRIDSTEP,MAXy-GRIDSTEP,GRIDSTEP,MAXy-GRIDSTEP); 

line(HAXx-GRIDSTEP,MAXy-GRIDSTEP,MAXx-GRIDSTEP,GRIDSTEP); 

SetColor(TEXTC0L0R-2); 

SetTextStyle<1,0,1); 

SetText Just ify( 1,1); 

outtextxy(round(MAXx/2),GRIDSTEP+30, 

'       Non  Reactive  Target  Connputer  Assisted  Search    '); 

SetTextStyle(2,0,4); 

outtextxy(round(MAXx/2),GRIDSTEP+90, 

'    This  program  will  perform  a  Montecarlo  Simulation  to  assist  in  ' 

outtextxy(round(MAXx/2),GRIDSTEP+100, 

'    conduction  of  a  Lost  Contact  Search  of  a  Non  Reactive  Target  by  ' 

outtextxy(round(MAXx/2),GRIDSTEP+110, 

'    Passive  Searchers'); 

outtextxy(round(MAXx/2),GRIDSTEP+130, 

'    The  user  must  provide  parameters  regarding:  ' 

outtextxy(round(MAXx/2),GRIDSTEP+150, 

'  THE  TARGET  ' 

outtextxy(round(MAXx/2),GRIDSTEP+160, 

'  THE  SEARCHER  ' 

outtextxy(round(MAXx/2),GRIDSTEP+170, 

'  THE  SIMULATION  ' 

outtextxy(round(MAXx/2),GRIDSTEP+210, 

'    The  program  will  output  a  Probability  Distribution  Map  in  Color' 

out textxy( round( MAXx/2 ) , GR I DSTEP+220 , 

'    Codes  relative  to  the  maximum  value  of  probability  encountered' 

outtextxy(round(MAXx/2),GRIDSTEP+230, 

'    in  any  of  the  square  cells  in  wich  the  Local  Area  is  divided.  ' 

outtextxy(round(MAXx/2),GRIDSTEP+300, 

'    Please,  press  any  key  to  continue'); 

readkey; 

end: 


procedure  Start: 

begin 

clrscr; 

InitGraphics; 

Di spl ay Logo; 

<:delay(4000);> 

Displaylntro; 
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Closegraph; 
end; 


procedure  SetTi  me ( t : i  nteger ) ; 

var 

fstdigit     :byte; 

snddigit     :byte; 

trddigit     :byte; 

begin 

t:=t*5; 

SetTextStyle(2,0,5); 

OutTextXY<STEP,FRAME*STEP+STEP*5,  'time  :  '); 

SetColor(TEXTC0L0R-2); 

fstdigit :=round(t  div  100); 

snddigit :=round((t  mod  100)  div  10)  ; 

trddigit:=round((t  mod  100)  mod  10)  ; 

OutTextXY(STEP+STEP*1 1 , FRAME*STEP+STEP*5, 

concat(chr(48+fstdigit),chr(48+snddigit),chr(48+trddigit))); 
end; 


procedure  OutputResults 

var 

MAXx 

: integer; 

MAXy 

: integer; 

fstdigit 

:byte; 

snddigit 

:byte; 

trddigit 

:byte; 

fthdigit 

:byte; 

f ihdigit 

:byte; 

P 

:real; 

count 

:byte; 

holder 

: integer; 

stringholder  rstring; 

begin 

MAXx:=GetMaxX; 

MAXy:=GetMaxY; 

SetPalette(15,GRI DCOLOR ) ; 

Setcolor(15); 

line(FRAME*STEP+5,STEP,MAXx-10,STEP); 

I i  ne( FRAME*STEP+5 , STEP, FRAME*STEP+5,MAXy- 10); 

line(MAXx-10,MAXy-10,STEP*FRAME+5,MAXy-10); 

line(MAXx-10,MAXy-10,MAXx-10,STEP); 

SetColor(TEXTC0L0R-2); 

OutTextXY(FRAME*STEP+15,STEP*3,'  RTCAS  STATUS   '); 

SetColordS) 

SetTextStyle<2,0,4); 

OutTextXY(FRAME*STEP+15,STEP*28,'UNDET  PROB  MASS   '); 

p: =undetecte<iiiass*1 0000; 

fstdigit:=round<p)  div  10000; 

snddigit:=(round(p)  mod  10000)div  1000; 

trddigit:=((round(p)  mod  10000)mod  1000)div  100; 

fthdigit:=(((round(p)  mod  10000)mod  1000)mod  100)div  10; 

fihdigit:=((<round(p)  mod  10000)mod  1000)mod  100)mod  10; 

SetColor(TEXTCOLOR-2); 

OutTextXY(FRAME*STEP+STEP*30,STEP*28, 

concat(chrC48+fstdigit),chr(46),chr(48+snddigit), 
chr(48+trddigit),chr(48+fthdigit),chr(48+f ihdigit))); 

SetColordS)  ; 

SetTextStyle(2,0,4); 
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OutTextXY(FRAME*STEP+15,STEP*32, 'Number  of  TGTS  OFF  '); 
fstcligit:=round(numberout  div  100); 
snddigit:=round((numberout  mod  100)  div  10)  ; 
trddigit:=round((nunberout  mod  100)  mod  10)  ; 
SetColor(TEXTC0L0R-2); 
OutTextXY(FRAME*STEP+STEP*30,STEP*32, 

concat(chr(48+fstdigit),chr(48+snddigit),chr(48+trddigit))); 

SetPalette(15,GRIDC0L0R); 

setcolor(15); 

line(FRAME*STEP+5,STEP*40,MAXx-10,STEP*40); 

SetTextStyle(2,0,4); 

SetColor(TEXTC0L0R-2); 

OutTextXY(FRAME*STEP*15,STEP*42,'      MC  SIMULATION  INFO'); 

SetColor(15); 

OutTextXY ( FRAME*STEP+1 5 , STEP*45 , ' TARGET  DATA ' ) ; 
OutTextXY<FRAME*STEP+15,STEP*48, 'DATUM  ERROR  (rai)  '); 
fstdigit:=round(sigmarinit  div  100); 
snddigit:=round((sigmarinit  mod  100)  div  10)  ; 
trddigit:=round((sigmarinit  mod  100)  mod  10)  ; 
SetColor(TEXTC0L0R-2); 
OutTextXY(FRAME*STEP+STEP*30,STEP*48, 

concat(chr(48+fstdigit),chr(48+snddigit), 
chr(48+trddigit))); 

SetColordS); 

OutTextXY(FRAME*STEP+15,STEP*50, 'REACTION  RANGE  (rmi)  '); 

fstdigit:=round(reactrange  div  100); 

snddigit:=round((reactrange  mod  100)  div  10)  ; 

trddigit:=routxl(( react range  mod  100)  mod  10)  ; 

SetColor(TEXTC0L0R-2); 

OutTextXY(FRAME*STEP+STEP*30,STEP*50, 

concat(chr(48+fstdigit),chr(48+snddigit), 
chr(48+trddigit))); 


SetColordS); 

OutTextXY(FRAME*STEP+15,STEP*52, 'REACTION  TACK  '); 

fstdigit:=round(reactiontack  div  100); 

snddigit:=round((reactiontack  mod  100)  div  10)  ; 

trddigit:=round((reactiontack  mod  100)  mod  10)  ; 

SetColor(TEXTC0L0R-2); 

OutTextXY(FRAME*STEP+STEP*30,STEP*52, 

concat(chr(48+fstdigit),chr(48+snddigit), 
chr(48+trddigit))); 


SetColordS); 

OutTextXY<FRAME*STEP+1S,STEP*S4,'MIN  SPEED         '); 

fstdigit:=round(minsubspeed  div  100); 

snddigit:=round((minsubspeed  mod  100)  div  10)  ; 

trddigit:=round((minsubspeed  mod  100)  mod  10)  ; 

SetColor(TEXTC0L0R-2); 

OutTextXY (FRAME*STEP+STEP*30,STEP*S4, 

concat(chr(48+fstdigit),chr(48+snddigit), 
chr(48+trddigit))); 

SetColordS); 

OutTextXY(FRAME*STEP+1S,STEP*S6,'MAX  SPEED  '); 
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fstdigit:=round(maxsubspeed  div  100); 
snddigit:=round((inaxsubspeed  mod  100)  div  10)   ; 
trddigit:=round((maxsubspeed  mod  100)  mod  10)  ; 
SetCo I or<  TEXTCOLOR- 2 ) ; 
OutTextXY(FRAME*STEP+STEP*30,STEP*56, 

concat(chr(48+fstdigit),chr<48+snddigit), 
chr(48+trddigit))); 

SetColordS); 

line(FRAME*STEP+5,STEP*60,MAXx-10,STEP*60); 

SetColor(TEXTC0L0R-2); 

OutTextXY(FRAME*STEP+15,STEP*62,'     REAL  WORLD  SEARCH  INFO  '); 

SetColordS); 

OutTextXY(FRAME*STEP+15,STEP*65, 'SEARCHERS  DATA'); 

0utTextXY(FRAME*STEP+15,STEP*69,'SC0UT  COURSE  SPEED  SONAR  HOT  '); 

for  count :=1  to  8  do  begin 

if  count  <=  nimberof scouts  then  begin 
SetColor<TEXTC0L0R-2); 
OutTextXY(FRAME*STEP+15,STEP*70+count*2*STEP,   concat( 

myscout [count] .id,chr(48),chr(48),chr(48+count),'  ')); 
fstdigit:=round(myscout [count] .course  div  100); 
snddigit:=round((myscout [count] .course  mod  100)  div  10)  ; 
trddigit:=round((m/scout [count] .course  mod  100)  mod  10)  ; 
OutTextXY ( FRAME*STEP+55 , STEP*70+count*2*STEP , conca t ( 

chr(48+fstdigit),chrC48+snddigit),chr(48+trddigit))); 
holder:=round(myscout [count] .speed); 
fstdigit:=round(holder  div  100); 
snddigit:=round((holder  mod  100)  div  10)  ; 
trddigit:=round((holder  mod  100)  mod  10)  ; 
OutTextXY ( FRAME*STEP+95 , STEP*70+count*2*STEP , concat ( 

chr{48+fstdigit),chr(48+snddigit),chr(48+trddigit))); 
if  myscout [count] .SONAR_ON  then  begin 

stringholder:='  ON  '; 
end  else  begin 

stringholder:='  OFF'; 
end; 

OutTextXY(FRAME*STEP+135,STEP*70+count*2*STEP,stringholder); 
<***>  myscout[2].IS_H0T:=TRUE; 

if  myscout [count] . I S_H0T  then  begin 
stringholder:='  YES'; 
SetColor(14); 

OutTextXY(FRAME*STEP+175,STEP*70+count*2*STEP,stringholder); 
end; 
if  myscout [count] .initialhot  <>  0  then  begin 

D  i  spl ayContactSimnaryC  count ) ; 
end; 
end; 
end; 

SetTextStyle(2,0,4); 
SetColordS); 

OutTextXY(FRAME*STEP+15,MAXy-30,'Hit  any  key  to  go  back  to  main  menu'); 
end; {procedure} 


procedure  Displa>ContactSunanary(srchr:byte) : 

var 

holder: integer; 

fstdi git: integer; 

snddi  gi  t : i nteger; 

trddigit:integer; 
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begin 

SetColor(TEXTC0L0R-2); 

OutTextXY(30*STEP,FRAME*STEP+STEP,' INITIAL  CONTACT  REPORTS   '); 

SetColordS); 

OutTextXY(30*STEP,FRAME*STEP+STEP*3,'SEAR     TIME     COURSE     SPEED    '); 

SetColor(14); 

OutTextXY(30*STEP,FRAME*STEP+STEP*3+STEP*srchr*2,concat( 

myscouttsrchr] .id,chr(48),chr(48),chr(48+srchr>)); 
fstdigit:=round(tnyscout[srchr3 .initialhot  div  100); 
snddigit:=round((myscout[srchr] .initialhot  mod  100)  div  10)  ; 
trddigit:=round((myscout[srchr] .initialhot  n»d  100)  mod  10)  ; 
OutTextXY(30*STEP+35,FRAME*STEP+STEP*3+STEP*srchr*2,concat( 

chr(48+fstdigit),chr(48+snddigit),chr(48+trddigit))); 
fstdigit:=round(myscout[srchr] .subcourse  div  100); 
snddigit:=round((myscout[srchr] .subcourse  mod  100)  div  10)  ; 
trddigit:=round((myscout[srchr] .subcourse  mod  100)  mod  10)  ; 
OutTextXY(30*STEP+80,FRAME*STEP+STEP*3+STEP*srchr*2,concat( 

chr(48+fstdigit),chr(48+snddigit),chr(48+trddigit))); 
holder:=round(myscout[srchr] .subspeed); 
fstdigit:=rowxi(holder  div  100); 
sr>ddi9it:=round((holder  mod  100)  div  10)  ; 
trddigit:=round((holder  mod  100)  mod  10)  ; 
OutTextXY(30*STEP+125,FRAME*STEP+STEP*3+STEP*srchr*2,concat( 

chr<48+fstdigit),chr(48+snddigit),chr(48+trddigit))); 
end; 


procedure  DisplayData; 

begin 

ClearDevice; 

SetGrid; 

SetScale; 

DisplayTargets; 

DisplayScouts; 

undetectedmass:=1 ; 

OutputResults; 

readkey; 

cleardevice; 

CloseGraph; 

end: 


procedure  OutputGraph; 

begin 

ClearDevice; 

SetGrid; 

SetScale; 

SetTime(totalsimtime); 

<DispLayTargets;> 

SetProbOistrColors; 

DisplayScouts; 

OutputResults; 

readkey; 

ClearDevice; 

CloseGraph; 

end; 

begin 
end. 
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